| 查看: 450 | 回复: 4 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请教fortran程序出错问题
|
|||
|
问题一：程序如下： program initial_velocity implicit none !we assumed m=1,k=1 real*8,parameter::temp=300.0 integer,parameter::N=10 integer::i real*8,dimension(1:N)::vx,vy,vz write(*,'(1x,a,8x,a,10x,a,10x,a)')'n','vx','vy','vz' call comvel(temp) do i=1,N write(*,'(i2,2x,f10.3,2x,f10.3,2x,f10.3)')i,vx(i),vy(i),vz(i) end do end program initial_velocity subroutine comvel(temp) integer,parameter::N=10 integer::i real*8,dimension(1:N)::vx,vy,vz real*8::temp,Rtemp,sumx,sumy,sumz real*8::gauss,dummy Rtemp=sqrt(temp) do i=1,N vx(i)=Rtemp*gauss(dummy) vy(i)=Rtemp*gauss(dummy) vz(i)=Rtemp*gauss(dummy) end do sumx=0.0 sumy=0.0 sumz=0.0 do i=1,N sumx=sumx+vx(i) sumy=sumy+vy(i) sumz=sumz+vz(i) end do sumx=sumx/real(N) sumy=sumy/real(N) sumz=sumz/real(N) do i=1,N vx(i)=vx(i)-sumx vy(i)=vy(i)-sumy vz(i)=vz(i)-sumz end do return end subroutine comvel function gauss(dummy) real*8,parameter::A1=3.949846138,A3=0.252408784 real*8,parameter::A5=0.076542912,A7=0.008355968 real*8,parameter::A9=0.029899776 real*8::sum,R,R2,dummy integer::i sum=0.0 call random_seed() do i=1,12 call random_number(dummy) sum=sum+dummy end do R=(sum-6.0)/4.0 R2=R*R gauss=((((A9*R2+A7)*R2+A5)*R2+A3)*R2+A1)*R return end function gauss 运行后结果： n vx vy vz 1 0.000 0.000 0.000 2 0.000 0.000 0.000 3 0.000 0.000 0.000 4 0.000 0.000 0.000 5 0.000 0.000 0.000 6 0.000 0.000 0.000 7 0.000 0.000 0.000 8 0.000 0.000 0.000 9 0.000 0.000 0.000 10 0.000 0.000 0.000 Press any key to continue 请问怎样能得到所需的随机数值？？ 问题二：为验证上面程序中部分程序段，我写另一程序如下： program a_1 implicit none real*8::gauss,dummy,Rtemp integer,parameter::N=10 integer::i real,dimension(1:N)::vx,vy,vz real,parameter::temp=300.0 Rtemp=sqrt(temp) do i=1,N vx(i)=Rtemp*gauss(dummy) vy(i)=Rtemp*gauss(dummy) vz(i)=Rtemp*gauss(dummy) print*,vx(i),vy(i),vz(i) end do end program a_1 function gauss(dummy) real*8,parameter::A1=3.949846138,A3=0.252408784 real*8,parameter::A5=0.076542912,A7=0.008355968 real*8,parameter::A9=0.029899776 real*8::sum,R,R2,dummy integer::i sum=0.0 do i=1,12 call random_seed() call random_number(dummy) sum=sum+dummy end do R=(sum-6.0)/4.0 R2=R*R gauss=((((A9*R2+A7)*R2+A5)*R2+A3)*R2+A1)*R return end function gauss 运行结果有时值全部一样，有时部分一样，但对应同一个i 时vx(i),vy(i),vz(i)始终值一样结果如下： 结果（1）： 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 49.20191 Press any key to continue 结果（2）： 82.36326 82.36326 82.36326 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 -72.15305 Press any key to continue 请多多指教！！ |
回复此楼» 猜你喜欢
评审有感
已经有10人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
博士申请
已经有5人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
2楼2009-04-14 18:01:29
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★
jianchaoyv(金币+5,VIP+0): 4-14 21:25
jianchaoyv(金币+5,VIP+0): 4-14 21:25
|
第一个问题：为什么输出全是零呢？因为 call comvel(temp) do i=1,N write(*,'(i2,2x,f10.3,2x,f10.3,2x,f10.3)')i,vx(i),vy(i),vz(i) end do 在这里，你没有把vx,vy,vz传到子程序里，当然输出的都是零。 第二个问题，我建议你运行以下这两个小程序，测试一下 random_seed()和random_number(dummy)的用法。 程序一：（这是你的程序写法） do i=1,12 call random_seed() call random_number(dummy) write(*,*) dummy enddo end 得到什么？呵呵，全都一样的随机数。但是，每运行一次，就有不同的序列。 程序二： 把random_seed()放在外面： call random_seed() do i=1,12 call random_number(dummy) write(*,*) dummy enddo end 得到什么？不同的随机数！ 所以，你应该把 call random_seed()写在循环的外边！ 但是，你把Call random_seed()写在子程序中循环的外边之后，又如何呢？ 又是同样的随机数！！！为什么呢？因为你在主程序中反复调用了Gauss，还是相当于把Call random_seed()写在了循环内部。 最后我们知道，正确的写法应该把Call random_seed()写在主程序的开头！！！呵呵。 注：我只管编程，没有考察算法的正确与否。 |
3楼2009-04-14 19:45:56
4楼2009-04-14 21:46:05
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
5楼2009-04-15 08:32:16
30