 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1727 | 回復(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ù)此樓» 猜你喜歡
面上5B能上會(huì)嗎？
已經(jīng)有4人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有9人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。
已經(jīng)有4人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有14人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有5人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有21人回復(fù)
-
26申博自薦
已經(jīng)有5人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】求fortra程序
已經(jīng)有4人回復(fù)
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
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
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
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ā)表 | |
|---|---|---|---|---|
|
[考研] 總分293求調(diào)劑 +7 | 加一一九 2026-03-25 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 085600 286分 材料求調(diào)劑 +8 | 麻辣魷魚 2026-03-27 | 9/450 |
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4 | 12.27 2026-03-27 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +9 | 15188958825 2026-03-25 | 9/450 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 0703一志愿9，初試成績(jī)：338，四六級(jí)已過(guò)，有科研經(jīng)歷，求調(diào)劑！ +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過(guò) 2026-03-27 | 6/300 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|