irbit1.f90

来自「FORTRANvisualfortran常用数值算法集及源码」· F90 代码 · 共 16 行

FUNCTION irbit1(iseed)
INTEGER irbit1,iseed,IB1,IB2,IB5,IB18
PARAMETER (IB1=1,IB2=2,IB5=16,IB18=131072)
LOGICAL newbit
newbit=iand(iseed,IB18)/=0
if(iand(iseed,IB5)/=0) newbit=.not.newbit
if(iand(iseed,IB2)/=0) newbit=.not.newbit
if(iand(iseed,IB1)/=0) newbit=.not.newbit
irbit1=0
iseed=iand(ishft(iseed,1),not(IB1))
if(newbit) then
  irbit1=1
  iseed=ior(iseed,IB1)
endif
END FUNCTION irbit1