songjunann

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 16
散金: 4
帖子: 164
在線: 94.5小時(shí)
蟲號(hào): 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ù)此樓

» 猜你喜歡

311求調(diào)劑已經(jīng)有4人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)
354求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有3人回復(fù)
274求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有15人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有8人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有10人回復(fù)
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有8人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有10人回復(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的回帖

snoopyzhao

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

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲號(hào): 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的回帖

songjunann

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 16
散金: 4
帖子: 164
在線: 94.5小時(shí)
蟲號(hào): 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的語法，不會(huì)沖突吧？

贊一下

回復(fù)此樓

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

已閱回復(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í)
蟲號(hào): 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的語法，不會(huì)沖突吧？

應(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í)
蟲號(hào): 1198651
注冊: 2011-02-02
性別: MM
專業(yè): 高分子物理與高分子物理化

引用回帖:

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

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

贊一下

回復(fù)此樓

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

已閱回復(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í)
蟲號(hào): 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的大小是動(dòng)態(tài)的嗎？

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

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 271材料工程求調(diào)劑 +4	.6lL 2026-03-18	4/200	2026-03-18 20:05 by 楤哥
[考研] 314求調(diào)劑 +8	無懈可擊的巨人 2026-03-12	8/400	2026-03-18 14:50 by haxia
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 311求調(diào)劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 0703化學(xué)調(diào)劑 +4	pupcoco 2026-03-17	7/350	2026-03-18 12:14 by djl2006
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao