當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

songjunann

銅蟲 (小有名氣)

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

[求助] 幫忙看一下，這兩個(gè)程序（關(guān)于點(diǎ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子程序?yàn)椋?br />    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

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

回復(fù)此樓

» 猜你喜歡

299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有9人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有8人回復(fù)

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

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

songjunann

銅蟲 (小有名氣)

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

引用回帖:

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

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

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

贊一下

回復(fù)此樓

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

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

查看全部 6 個(gè)回答

snoopyzhao

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

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

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

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

snoopyzhao

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

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

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

jjdg: 感謝參與 2011-09-27 00:42:46
songjunann(金幣+20): 謝謝，再問一個(gè)小問題， 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小時(shí)
蟲號: 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的回帖

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

[考研] 299求調(diào)劑 +8	15188958825 2026-03-25	8/400	2026-03-29 01:36 by fmesaito
[考研] 0856求調(diào)劑 +6	楒桉 2026-03-28	6/300	2026-03-29 00:31 by 544594351
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 315求調(diào)劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 083000學(xué)碩274求調(diào)劑 +8	Li李魚 2026-03-26	8/400	2026-03-28 20:33 by 加油向未來啊
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +4	傾____蕭 2026-03-27	5/250	2026-03-27 11:52 by zhshch
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 359求調(diào)劑 +4	王了個(gè)楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚的貓
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 352求調(diào)劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚的貓
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +6	cccchenso 2026-03-23	6/300	2026-03-25 22:25 by 544594351
[考研] 296求調(diào)劑 +4	汪��？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪��？！