songjunann

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 16
散金: 4
帖子: 164
在線: 94.5小時
蟲號: 1198651
注冊: 2011-02-02
性別: MM
專業(yè): 高分子物理與高分子物理化

[求助] 幫忙看一下，這兩個程序（關(guān)于點積的）講的是什么，為什么要寫的這么繁瑣

dot.f程序：

  function dot (x,y,n)
cc-begin
c
c-----------------------------------------------------------------------
c
c    dot (x,y,n)
c
c
c    purpose
c    -------
c
c    dot calculates the dot product of the two n-vectors
c    x and y.
c
cc-doc
c    input arguments
c    ---------------
c
c    x    = vector of length n.
c    y    = vector of length n.
c    n    = dimension of x and y.
c
c    output arguments
c    ----------------
c
c    dot    = value of the function, i.e. dot product of x and y.
c
c    notes
c    -----
c
c    1. dot contains a test for n = 0.
c
c-----------------------------------------------------------------------
cc-end
c
   implicit real*8 (a-h,o-z)
   dimension x(*),y(*)
c
   dot = 0.
   if (n .le. 0) return
c    do 100 i = 1,n
c       dot = dot + x(i)*y(i)
c  100 continue
   dot = ddot(n,x,1,y,1)
   end
%%%%%%%%%%%%%%%%%%%%%%%%%
其中ddot.f子程序為：
   DOUBLE PRECISION FUNCTION DDOT(N,DX,INCX,DY,INCY)
*    .. Scalar Arguments ..
   INTEGER INCX,INCY,N
*    ..
*    .. Array Arguments ..
   DOUBLE PRECISION DX(*),DY(*)
*    ..
*
*  Purpose
*  =======
*
*    DDOT forms the dot product of two vectors.
*    uses unrolled loops for increments equal to one.
*
*  Further Details
*  ===============
*
*    jack dongarra, linpack, 3/11/78.
*    modified 12/3/93, array(1) declarations changed to array(*)
*
*  =====================================================================
*
*    .. Local Scalars ..
   DOUBLE PRECISION DTEMP
   INTEGER I,IX,IY,M,MP1
*    ..
*    .. Intrinsic Functions ..
   INTRINSIC MOD
*    ..
   DDOT = 0.0d0
   DTEMP = 0.0d0
   IF (N.LE.0) RETURN
   IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
*
*       code for both increments equal to 1
*
*
*       clean-up loop
*
      M = MOD(N,5)
      IF (M.NE.0) THEN
         DO I = 1,M
            DTEMP = DTEMP + DX(I)*DY(I)
         END DO
         IF (N.LT.5) THEN
            DDOT=DTEMP
         RETURN
         END IF
      END IF
      MP1 = M + 1
      DO I = MP1,N,5
      DTEMP = DTEMP + DX(I)*DY(I) + DX(I+1)*DY(I+1) +
   $          DX(I+2)*DY(I+2) + DX(I+3)*DY(I+3) + DX(I+4)*DY(I+4)
      END DO
   ELSE
*
*       code for unequal increments or equal increments
*       not equal to 1
*
      IX = 1
      IY = 1
      IF (INCX.LT.0) IX = (-N+1)*INCX + 1
      IF (INCY.LT.0) IY = (-N+1)*INCY + 1
      DO I = 1,N
         DTEMP = DTEMP + DX(IX)*DY(IY)
         IX = IX + INCX
         IY = IY + INCY
      END DO
   END IF
   DDOT = DTEMP
   RETURN
   END

個人感覺不就是兩個向量的點積，為什么兩個程序要寫的這么復(fù)雜呢，有什么特別的地方嗎，請大家多多指教，謝謝！

回復(fù)此樓

» 猜你喜歡

【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有4人回復(fù)
317求調(diào)劑已經(jīng)有5人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有10人回復(fù)
一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。已經(jīng)有3人回復(fù)
一志愿福大288有機化學(xué)，求調(diào)劑已經(jīng)有3人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有10人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有11人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

請幫忙看一下教授的回復(fù)是什么意思？已經(jīng)有5人回復(fù)
請問對于生產(chǎn)制劑的原料藥產(chǎn)地變更的申請，其變更說明該怎么寫已經(jīng)有8人回復(fù)
求高手幫忙鑒定一下植物已經(jīng)有14人回復(fù)
求助高斯腳本同時提交兩個任務(wù)，一個結(jié)束，另一個也結(jié)束了已經(jīng)有8人回復(fù)
幫忙看一下這兩個粒徑分布圖的區(qū)別已經(jīng)有3人回復(fù)
【急】請大家?guī)兔匆幌逻@個offer是否合乎要求，新手散光所有金幣求助~~ 已經(jīng)有12人回復(fù)
鎂基復(fù)合材料的XRD圖高手幫忙看一下啊急�。�！已經(jīng)有12人回復(fù)
幫忙看一段程序已經(jīng)有3人回復(fù)
大家?guī)兔匆幌? 這個反應(yīng)為什么不能進行已經(jīng)有6人回復(fù)
請大俠幫忙看一下這個包裝和內(nèi)含的是什么材料已經(jīng)有7人回復(fù)
有這方面經(jīng)驗的人幫忙給介紹一下吧已經(jīng)有6人回復(fù)
初寫論文菜鳥，大家?guī)兔匆幌逻@幾張圖是用什么軟件畫的或生成的？已經(jīng)有10人回復(fù)
文書寫作經(jīng)驗交流月——分享幾個PS寫作指導(dǎo)網(wǎng)站已經(jīng)有8人回復(fù)
【求助/交流】大家?guī)兔匆幌逻@是什么菌已經(jīng)有55人回復(fù)

1樓 2011-09-25 14:58:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★
xzhdty(金幣+1): 歡迎常來程序語言看看 2011-09-25 19:19:42

ddot.f 應(yīng)該是 blas 包中的函數(shù)，之所以寫成這樣，呵呵，或許是受那個時代計算機的計算能力的限制吧，里面有一些手工的向量化操作，呵呵……

在 Fortran 90 后，可以用內(nèi)置函數(shù)dot_product來做這個工作了……

贊一下(1人)

回復(fù)此樓

2樓2011-09-25 16:43:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

songjunann

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 16
散金: 4
帖子: 164
在線: 94.5小時
蟲號: 1198651
注冊: 2011-02-02
性別: MM
專業(yè): 高分子物理與高分子物理化

引用回帖:

2樓: Originally posted by snoopyzhao at 2011-09-25 16:43:23:
ddot.f 應(yīng)該是 blas 包中的函數(shù)，之所以寫成這樣，呵呵，或許是受那個時代計算機的計算能力的限制吧，里面有一些手工的向量化操作，呵呵……

在 Fortran 90 后，可以用內(nèi)置函數(shù)dot_product來做這個工作了……

那在Visual Compaq Fortran里面，混合用F77和F90的語法，不會沖突吧？

贊一下

回復(fù)此樓

3樓2011-09-25 23:09:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

jjdg: 感謝參與 2011-09-27 00:42:46
songjunann(金幣+20): 謝謝，再問一個小問題， 2011-09-27 01:00:08

引用回帖:

3樓: Originally posted by songjunann at 2011-09-25 23:09:13:
那在Visual Compaq Fortran里面，混合用F77和F90的語法，不會沖突吧？

應(yīng)該沒有什么問題，但從編程規(guī)范的角度來說，不建議這樣做……

贊一下(1人)

回復(fù)此樓

4樓2011-09-25 23:49:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

songjunann

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 16
散金: 4
帖子: 164
在線: 94.5小時
蟲號: 1198651
注冊: 2011-02-02
性別: MM
專業(yè): 高分子物理與高分子物理化

引用回帖:

3樓: Originally posted by songjunann at 2011-09-25 23:09:13:
那在Visual Compaq Fortran里面，混合用F77和F90的語法，不會沖突吧？

F77里面：
implicit real*8 (a-h,o-z)
dimension a(*),b(*),mdiag(*)
這樣的定義用在哪里呢，a,b的大小是動態(tài)的嗎？

贊一下

回復(fù)此樓

5樓2011-09-27 01:01:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
余澤成(金幣+1): 謝謝參與應(yīng)助！ 2011-09-28 10:31:43

引用回帖:

5樓: Originally posted by songjunann at 2011-09-27 01:01:29:
F77里面：
implicit real*8 (a-h,o-z)
dimension a(*),b(*),mdiag(*)
這樣的定義用在哪里呢，a,b的大小是動態(tài)的嗎？

這個不叫動態(tài)的，而叫 assumed size array，一般用在子程序中，用以表示 dummy variable 的維度是由實際傳遞的參數(shù)決定的。

贊一下(1人)

回復(fù)此樓

6樓2011-09-27 01:28:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +4	吃吃吃才有意義 2026-03-19	4/200	2026-03-19 13:48 by houyaoxu
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	6/300	2026-03-19 13:24 by 楓橋ZL
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	10/500	2026-03-19 13:23 by houyaoxu
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3	Ncdx123456 2026-03-19	3/150	2026-03-19 13:18 by houyaoxu
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 求材料調(diào)劑 +10	隔壁陳先生 2026-03-12	10/500	2026-03-19 09:15 by 無際的草原
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	9/450	2026-03-18 15:31 by vgtyfty
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz