3

我正在使用数字配方方案来生成随机数（此 PDF 文件中的ran3第 7 页）。我没有发现任何奇怪的地方，但是这一次，我在“热身”阶段得到了一个大于. 代码看起来好像这不应该发生。我可以通过将 if 语句更改为 while 语句来轻松解决此问题，但我想知道后果是什么。MBIGif(mk.lt.MZ)mk=mk+MBIG

编辑：这是功能

      FUNCTION ran3a(idum) 
      INTEGER idum
      INTEGER MBIG,MSEED,MZ 
C     REAL MBIG,MSEED,MZ
      REAL ran3a,FAC
      PARAMETER (MBIG=1000000000,MSEED=161803398,MZ=0,FAC=1./MBIG)  
C     PARAMETER (MBIG=4000000.,MSEED=1618033.,MZ=0.,FAC=1./MBIG) 
      INTEGER i,iff,ii,inext,inextp,k
      INTEGER mj,mk,ma(55)  
C     REAL mj,mk,ma(55) 
      SAVE iff,inext,inextp,ma
      DATA iff /0/
      if(idum.lt.0.or.iff.eq.0)then
        iff=1
        mj=MSEED-iabs(idum) 
        mj=mod(mj,MBIG) 
        ma(55)=mj
        mk=1
        do 11 i=1,54
          ii=mod(21*i,55) 
          ma(ii)=mk
          mk=mj-mk
          if(mk.lt.MZ)mk=mk+MBIG
          mj=ma(ii)  
11      continue
        do 13 k=1,4
          do 12 i=1,55
            ma(i)=ma(i)-ma(1+mod(i+30,55)) 
            if(ma(i).lt.MZ)ma(i)=ma(i)+MBIG 
12        continue 
13      continue 
        inext=0
        inextp=31
        idum=1
      endif
      inext=inext+1
      if(inext.eq.56)inext=1
      inextp=inextp+1
      if(inextp.eq.56)inextp=1
      mj=ma(inext)-ma(inextp) 
      if(mj.lt.MZ)mj=mj+MBIG
      ma(inext)=mj
      ran3a=mj*FAC
      return
      END
4

2 回答 2

3

我得到了 Seg Faults（使用 gfortran 4.8），因为该函数试图将输入值idum从负数更改为 1。该行没有理由（也没有任何内容iff），所以我将其删除并打印出ma数组几个不同的地方，在数组中没有找到负数。

但是，一种可能性是，如果iabs(idum)大于MSEED，您可能对线有问题mj=MSEED - iabs(idum)。你应该mj=abs(MSEED-abs(idum))像书中写的那样使用来保护它。

于 2013-05-28T11:53:48.107 回答
-3

看了pdf。你需要做的是

1) 播种： value = ran3(-1)

2）使用它：值=ran3（0）

于 2013-05-27T08:08:52.367 回答