5 * Revision 1.1.1.1 1996/02/15 17:50:19 mclareni
11 C CERN PROGLIB# V104 RNDM .VERSION KERNLNX 1.00 930507
12 C ORIG. 22-MAR-88 from IBM version of G.Marsaglia et al., Montreal
13 C MODIF 24/06/92, adapted N.Karpenko, Dubna
15 C !!! Version valid only for 32 bit machines !!!
16 C !!! with IEEE floating point representation !!!
18 C- Uniform Random Number Generator,
19 C- giving the same sequence as the IBM and VAX version
21 EQUIVALENCE (AMAN,MANT)
22 PARAMETER (MSK1 = 201326592)
23 PARAMETER (MSK2 = 855638016)
24 C PARAMETER (MSK1 = '0C000000'X)
25 C PARAMETER (MSK2 = '33000000'X)
31 MANT = ISHFT (MCGN,-8)
32 IF (MANT.EQ.0) GO TO 14
34 C- AMAN in the range 1 to 2**24-1
36 C- multiply by 2.**(-24)
40 C-- for zero set RNDM = 2.**(-25)