PYTHIA/pythia/rkrand.F

   1
   2 *==================================================================
   3
   4       FUNCTION RKRAND(IDUMMY)
   5       IMPLICIT DOUBLE PRECISION(A-H,O-Z)
   6       SAVE
   7       DATA INIT/0/
   8       IF(INIT.EQ.0) THEN
   9         INIT=1
  10         X=DMOD(DSQRT(2D0),1D0)
  11         Y=DMOD(DSQRT(3D0),1D0)
  12         Z=DMOD(DSQRT(5D0),1D0)
  13       ELSE
  14         X=DMOD(X+Y+Z,1D0)
  15         Y=DMOD(X+Y+Z,1D0)
  16         Z=DMOD(X+Y+Z,1D0)
  17       ENDIF
  18       RKRAND=X
  19       END