一个关于dlog算不出来的问题

woshixiaobai

求路过的大神看一下为什么这样算不出来abc1，是定义的不够合理吗

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program main

02	implicit real(kind=8)(a-h,o-z)

03

04	pi=3.1416926d0

05	vm_ini=2.0d0

06	abc1=dsqrt(-dlog(rf(0)))

07	abc2=2.0d0*pi*rf(0)

08	p=abc1*dsin(abc2)*vm_ini

09

end

10

11	function rf(idum)

12	implicit real(kind=8)(a-h,o-z)

13	save ma,inext,inextp

14	parameter (mbig=1000000000,mseed=161803398,mz=0,fac=1.e-9)

15	dimension ma(55)

16	data iff/0/

17	if(idum<0.or.iff==0) then

18

iff=1

19	mj=mseed-iabs(idum)

20	mj=mod(mj,mbig)

21

ma(55)=mj

22

mk=1

23	do 50 i=1,54

24	ii=mod(21*i,55)

25

ma(ii)=mk

26

mk=mj-mk

27	if(mk<mz) mk=mk+mbig

28

mj=ma(ii)

29	50      continue

30	do 100 k=1,4

31	do 60 i=1,55

32	ma(i)=ma(i)-ma(1+mod(i+30,55))

33	if(ma(i)<mz) ma(i)=ma(i)+mbig

34	60          continue

35	100     continue

36

inext=0

37

inextp=31

38

endif

39	200 inext=inext+1

40	if(inext==56) inext=1

41	inextp=inextp+1

42	if(inextp==56) inextp=1

43	mj=ma(inext)-ma(inextp)

44	if(mj<mz) mj=mj+mbig

45	ma(inext)=mj

46

rf=mj*fac

47	if(rf>1.e-8.and.rf<0.999999999) return

48

goto 200

49

pause

50

end

woshixiaobai

吉大渣渣-固体 发表于 2018-4-17 01:07
计算不出来是因为您没有写输出语句，无法看到执行结果。

直接加上对应的输出语句即可。But重点来了，您的 ...

万分感谢！

吉大渣渣-固体

计算不出来是因为您没有写输出语句，无法看到执行结果。

直接加上对应的输出语句即可。But重点来了，您的代码还是有问题的，主要在声明和数据类型上。
1.“ implicit real(kind=8)(a-h,o-z)”这个懒人神句危害很大的，没有声明的以a-h,o-z开头的变量数组都可以在程序中畅通无阻执行，这样写错了个字母很难排查出来且编译器不会报错，其初始化也是问题，有重大安全隐患；强烈建议改implicit none,强迫症保平安；
2.请重视数据类型。
针对以上问题不妨做以下修改：
test.f90 (1.42 KB, 下载次数: 0)

2018-4-17 01:05 上传

点击文件名下载附件

执行结果：

至于二者为什么执行结果不同，因时间关系，建议您检查初始化和数据类型。祝好！