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

songjunann

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 16
散金: 4
帖子: 164
在線: 94.5小時(shí)
蟲號(hào): 1198651
注冊(cè): 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ù)雜呢，有什么特別的地方嗎，請(qǐng)大家多多指教，謝謝！

回復(fù)此樓

» 猜你喜歡

考研調(diào)劑已經(jīng)有3人回復(fù)
317分一志愿南理工材料工程本科湖工大求調(diào)劑已經(jīng)有9人回復(fù)
297求調(diào)劑已經(jīng)有10人回復(fù)
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。已經(jīng)有10人回復(fù)
343求調(diào)劑已經(jīng)有5人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有7人回復(fù)
求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分已經(jīng)有6人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有5人回復(fù)
340求調(diào)劑已經(jīng)有5人回復(fù)
一志愿北化求調(diào)劑已經(jīng)有4人回復(fù)

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

請(qǐng)幫忙看一下教授的回復(fù)是什么意思？已經(jīng)有5人回復(fù)
請(qǐng)問對(duì)于生產(chǎn)制劑的原料藥產(chǎn)地變更的申請(qǐng)，其變更說明該怎么寫已經(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ù)
【急】請(qǐng)大家?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ù)
請(qǐng)大俠幫忙看一下這個(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
注冊(cè): 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的回帖

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

snoopyzhao

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

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

songjunann

銅蟲 (小有名氣)

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

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

[考研] 0856材料化工調(diào)劑總分330 +14	zhubinhao 2026-03-27	14/700	2026-03-29 10:01 by Sjndkwm
[考研] 332求調(diào)劑 +4	@MZB382400 2026-03-28	4/200	2026-03-28 21:02 by 唐沐兒
[考研] 343求調(diào)劑 +5	愛羈絆 2026-03-28	5/250	2026-03-28 20:53 by 唐沐兒
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3	nhj_ 2026-03-25	3/150	2026-03-28 13:24 by Iveryant
[考研] 291求調(diào)劑 +6	HanBeiNingZC 2026-03-24	6/300	2026-03-28 07:55 by baoball
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4	Higraduate 2026-03-22	8/400	2026-03-26 20:34 by Higraduate
[考研] 生物學(xué) 296 求調(diào)劑 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚的貓
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5	軒小寧—— 2026-03-26	5/250	2026-03-26 12:05 by i_cooler
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 332求調(diào)劑 +6	032500 2026-03-25	6/300	2026-03-25 22:45 by 418490947
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚(yáng)。 m
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin