 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1115 | 回復: 2 | |||
mplanck至尊木蟲 (職業(yè)作家)
之上又是啥子蟲蟲？
|
[交流]
【求助】fortran 程序windows 下編譯正確，而linux下出錯 已有2人參與
|
|
如題，一個cpc上下的程序，windows 下編譯得到結果正確， 放在linux下就出問題： 用f77, fort77 編譯，在兩個sub模塊結尾提示 Error on line 1679: Declaration error for rvv: adjustable dimension on non- argument Error on line 1679: wr_ardecls: nonconstant array size 但我的數組沒問題呀 再用gfortran，f95編譯，編譯通過，但結果不對 再用ifort編譯，編譯也通過，又出來一個不對的新結果 真搞不懂為什么這樣？ 麻煩蟲蟲門指點下，謝謝！ 附出錯部分代碼： 其中報錯的一個sub模塊代碼如下，1679行對應于最后一行，即‘end’所在行；再就是 整個代碼中沒有找到wr_ardecls這個東西，不知道咋搞的。非常感謝！ c------------------------------------------------------------------------ subroutine dcludcmp(a,n,np,indx,d,ierror) c this is the double complex version of the homonyme sbr. from Numerical c Recipes implicit double complex(a-h,o-q,s-z) implicit double precision (r) parameter (tiny=(1.0d-30,1.0d-30)) dimension a(np,np),indx(n),rvv(n) d=dcmplx(1.d0,0.d0) ierror=0 do 12 i=1,n raamax=0.d0 do 11 j=1,n if (cdabs(a(i,j)).gt.raamax) raamax=cdabs(a(i,j)) 11 continue if (raamax.eq.0.d0)then write (*,*) 'singular matrix in dcludcmp' ierror=1 return endif rvv(i)=1.d0/raamax 12 continue do 19 j=1,n if (j.gt.1) then do 14 i=1,j-1 sum=a(i,j) if (i.gt.1)then do 13 k=1,i-1 sum=sum-a(i,k)*a(k,j) 13 continue a(i,j)=sum endif 14 continue endif raamax=0.d0 do 16 i=j,n sum=a(i,j) if (j.gt.1)then do 15 k=1,j-1 sum=sum-a(i,k)*a(k,j) 15 continue a(i,j)=sum endif rdum=rvv(i)*cdabs(sum) if (rdum.ge.raamax) then imax=i raamax=rdum endif 16 continue if (j.ne.imax)then do 17 k=1,n dum=a(imax,k) a(imax,k)=a(j,k) a(j,k)=dum 17 continue d=-d rvv(imax)=rvv(j) endif indx(j)=imax if(j.ne.n)then if(cdabs(a(j,j)).eq.0.d0) a(j,j)=tiny dum=1.d0/a(j,j) do 18 i=j+1,n a(i,j)=a(i,j)*dum 18 continue endif 19 continue if(cdabs(a(n,n)).eq.0.d0) a(n,n)=tiny return end c--------------------------------------------------------------------------- |
回復此樓» 猜你喜歡
求調劑
已經有3人回復
-
265求調劑
已經有4人回復
-
085700資源與環(huán)境308求調劑
已經有6人回復
-
一志愿吉林大學材料學碩321求調劑
已經有12人回復
-
286分人工智能專業(yè)請求調劑愿意跨考！
已經有3人回復
-
329求調劑
已經有5人回復
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
已經有4人回復
-
材料學碩318求調劑
已經有5人回復
-
一志愿中國海洋大學，生物學，301分，求調劑
已經有6人回復
-
081700化工學碩調劑
已經有3人回復
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強帖: 16
- 應助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
余澤成(金幣+2):辛苦了！ 2010-07-17 10:46:06
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
余澤成(金幣+2):辛苦了！ 2010-07-17 10:46:06
|
不要用 f2c 編譯你的程序，用 gfortran 或其它正宗的 Fortran 編譯器，呵呵 用第二個錯誤提示 google 一下，你會發(fā)現，你不是遇到這個錯誤的第一人，呵呵…… f2c does not support automatic arrays; g77 does, as do most other f77 compilers (and all f90+). 至于結果不同，最大的可能是一些 local variable 沒有賦值，這在 Fortran，尤其是一些老的程序中是非常常見的。用 gfortran 編譯時加上 -fno-automatic -finit-local-zero 之類的參數試試…… |
2樓2010-07-17 08:56:26
mplanck
至尊木蟲 (職業(yè)作家)
之上又是啥子蟲蟲？
- 應助: 0 (幼兒園)
- 金幣: 22899.2
- 散金: 5216
- 紅花: 19
- 帖子: 3980
- 在線: 760.8小時
- 蟲號: 468030
- 注冊: 2007-11-27
- 性別: GG
- 專業(yè): 原子和分子物理
3樓2010-07-17 11:27:26
相關版塊跳轉
我要訂閱樓主 mplanck 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 0703化學調劑 ，六級已過，有科研經歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 287求調劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 材料工程專碩調劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 08工科 320總分 求調劑 +5 | 梨花珞晚風 2026-03-17 | 5/250 |
|
|
|
[考研]
 0854，計算機類招收調劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 274求調劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 0854控制工程 359求調劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 289求調劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 330求調劑 +3 | ?醬給調劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 304求調劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 311求調劑 +3 | 冬十三 2026-03-13 | 3/150 |
|