 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1728 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】fortran出錯(cuò) 編譯通過(guò) 運(yùn)行溢出
|
|||
|
subroutine adams(t,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2) external f,rkt dimension y(n),d(n),z(n,m),b(4,n),e(n),s(n),g(n) real*8 kp1,kp2 a=t do 10 i=1,n 10 z(i,1)=y(i) call f(t,y,n,d,kp1,kp2) do 20 i=1,n 20 b(1,i)=d(i) do 50 i=2,4 if (i.le.m)then t=a+(i-1)*h call rkt(t,h,y,n,f,eps,d,e,s,g,kp1,kp2) do 30 k=1,n 30 z(k,i)=y(k) call f(t,y,n,d,kp1,kp2) do 40 k=1,n 40 b(i,k)=d(k) end if 50 continue do 100 i=5,m do 60 j=1,n q=55*b(4,j)-59*b(3,j)+37*b(2,j)-9*b(1,j) y(j)=z(j,i-1)+h*q/24.0 b(1,j)=b(2,j) b(2,j)=b(3,j) b(3,j)=b(4,j) 60 continue t=a+(i-1)*h call f(t,y,n,d,kp1,kp2) do 70 k=1,n 70 b(4,k)=d(k) do 80 j=1,n q=9*b(4,j)+19*b(3,j)-5*b(2,j)+b(1,j) z(j,i)=z(j,i-1)+h*q/24.0 y(j)=z(j,i) 80 continue do 90 k=1,n 90 b(4,k)=d(k) 100 continue t=a return end subroutine rkt(t,h,y,n,f,eps,d,b,c,g,kp1,kp2) dimension y(n),d(n),a(4),b(n),c(n),g(n) real*8 kp1,kp2 h1=h m=1 p=1+eps x=t do 10 i=1,n 10 c(i)=y(i) 20 if (p.ge.eps)then a(1)=h1/2.0 a(2)=a(1) a(3)=h1 a(4)=h1 do 30 i=1,n g(i)=y(i) y(i)=c(i) 30 continue dt=h/m t=x do 80 j=1,m call f(t,y,n,d,kp1,kp2) do 40 i=1,n 40 b(i)=y(i) do 60 k=1,3 do 50 i=1,n y(i)=y(i)+a(k)*d(i) b(i)=b(i)+a(k+1)*d(i)/3.0 50 continue tt=t+a(k) call f(tt,y,n,d,kp1,kp2) 60 continue do 70 i=1,n 70 y(i)=b(i)+h1*d(i)/6.0 t=t+dt 80 continue p=0.0 do 90 i=1,n q=abs(y(i)-g(i)) if(q.gt.p) p=q 90 continue h1=h1/2.0 m=m+m goto 20 endif t=x return end program madams c ----------------------------------------------------------------------- c this rpogram is used to solve a systwem of differential equagtion c by using adams method c ------------------------------------------------------------------- external f dimension y(11),d(11),z(11,500000),b(4,11),e(11),s(11),g(11) real*8 wucha,kp1,kp2,t,tao,variable,Parameters open(26,file='variable.dat') open(28,file='Parameters.dat') open(10,file='wucha.dat') h=0.001d0 y(1)=0.01d0 y(2)=0.02d0 y(3)=0.03d0 y(4)=0.4d0 y(5)=0.2d0 y(6)=0.8d0 y(7)=0.8d0 y(8)=2.5d0 y(9)=0.9d0 y(10)=6d0 y(11)=1d0 m=500000 n=11 eps=0.00001d0 do kp1=1.50d0,1.60d0,0.1d0 do kp2=0.00d0,0.50d0,0.1d0 call adams(a,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2) do i=1,m t=(i-1)*h detae=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+ $ (z(3,i)-z(6,i))**2.0)**0.5 ! write(26,'(2f15.6)')t,detae ! write(28,'(6f15.6)')t,z(7,i),z(8,i),z(9,i),z(10,i),z(11,i) if(detae>50)then write(10,*)kp1,kp2,1 write(*,*)kp1,kp2,1 goto 20 elseif(detae>1.0d-3)then tao=t endif enddo if(tao<=50)then write(10,*)kp1,kp2,0 write(*,*)kp1,kp2,0 else write(10,*)kp1,kp2,1 write(*,*)kp1,kp2,1 endif 20 continue enddo enddo close(2) end subroutine f(t,y,n,d,kp1,kp2) c ---------------------------------------------------------------------- ! real*8 kp1,kp2 dimension y(n),d(n) double precision ex,ey,ez,u1,kp1,kp2 cccccccccccccc ex=y(1)-y(4) ey=y(2)+y(5) ez=y(3)-y(6) c u1=ey-ez+y(8)*(y(1)+y(4))*ex-y(10)*(y(1)+y(4))*ey+4d0*y(11)*ez c u2=0d0 cc u1=ey+y(8)*(y(1)+y(4))*ex-y(10)*(y(1)+y(4))*ey cc u2=-ex+4d0*y(11)*ex u1=ey+y(8)*(y(1)+y(4))*ex-y(10)*(y(1)+y(4))*ey+4d0*y(11)*ez u2=-ex ! do kp1=0.00d0,10.0d0,0.01d0 ! do kp2=0.00d0,10.0d0,0.01d0 ! kp1=2.0d0 ! kp2=0.5d0 cccccccccccccc d(1)=y(2)-1d0*y(1)**3+3d0*y(1)**2-y(3)+3d0 d(2)=1d0-5d0*y(1)**2-y(2) d(3)=0.006d0*(4d0*(y(1)+1.6d0)-y(3)) cccccccccccccccccccccccc d(4)=y(5)-y(7)*y(4)**3+y(8)*y(4)**2-y(6)+3d0+u1 d(5)=y(9)-y(10)*y(4)**2-y(5) d(6)=y(11)*(4d0*(y(4)+1.6d0)-y(6))+u2 cccccccccccccccccccccccccccccccccc d(7)=-kp1*kp2*y(1)**3*ex d(8)=kp1*kp2*y(1)**2*ex d(9)=kp1*kp2*ey d(10)=-kp1*kp2*y(1)**2*ey d(11)=kp1*kp2*(4d0*(y(1)+1.6d0)-y(3))*ez ! enddo ! enddo return end 高手看看 這個(gè)程序 出錯(cuò)了 為什么 溢出了��！  |
回復(fù)此樓» 猜你喜歡
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑!
已經(jīng)有5人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
275求調(diào)劑
已經(jīng)有20人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有9人回復(fù)
-
294分080500材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】求fortra程序
已經(jīng)有4人回復(fù)
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é)
4樓2010-09-14 22:43:47
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é)
2樓2010-09-14 20:22:34
maomao1210
金蟲 (正式寫手)
- 程序強(qiáng)帖: 5
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1432.8
- 散金: 242
- 紅花: 16
- 沙發(fā): 1
- 帖子: 991
- 在線: 442.1小時(shí)
- 蟲號(hào): 253215
- 注冊(cè): 2006-05-20
- 性別: MM
- 專業(yè): 考古理論
★
余澤成(金幣+1):謝謝參與應(yīng)助！ 2010-09-15 13:32:24
余澤成(金幣+1):謝謝參與應(yīng)助！ 2010-09-15 13:32:24
|
do kp1=1.50d0,1.60d0,0.1d0 do kp2=0.00d0,0.50d0,0.1d0 call adams(a,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2) do i=1,m t=(i-1)*h detae=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+ $ (z(3,i)-z(6,i))**2.0)**0.5 ! write(26,'(2f15.6)')t,detae ! write(28,'(6f15.6)')t,z(7,i),z(8,i),z(9,i),z(10,i),z(11,i) if(detae>50)then write(10,*)kp1,kp2,1 write(*,*)kp1,kp2,1 goto 20 elseif(detae>1.0d-3)then tao=t endif enddo if(tao<=50)then write(10,*)kp1,kp2,0 write(*,*)kp1,kp2,0 else write(10,*)kp1,kp2,1 write(*,*)kp1,kp2,1 endif 20 continue enddo enddo 在循環(huán)的時(shí)候，=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+ $ (z(3,i)-z(6,i))**2.0 是一個(gè)無(wú)窮大或者無(wú)窮小的數(shù)字，cvf不識(shí)別，不予計(jì)算，就提示溢出了。 gfortran 編譯計(jì)算的時(shí)候，即使是NaN 也會(huì)繼續(xù)計(jì)算不提示任何錯(cuò)誤。 |
3樓2010-09-14 22:22:30
maomao1210
金蟲 (正式寫手)
- 程序強(qiáng)帖: 5
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1432.8
- 散金: 242
- 紅花: 16
- 沙發(fā): 1
- 帖子: 991
- 在線: 442.1小時(shí)
- 蟲號(hào): 253215
- 注冊(cè): 2006-05-20
- 性別: MM
- 專業(yè): 考古理論
snoopyzhao ，嘿嘿，
5樓2010-09-14 22:57:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +9 | leaflight 2026-03-24 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|