| fe4da5cc |
1 | * |
| 2 | * $Id$ |
| 3 | * |
| 4 | * $Log$ |
| 5 | * Revision 1.1.1.1 1996/02/15 17:50:24 mclareni |
| 6 | * Kernlib |
| 7 | * |
| 8 | * |
| 9 | FUNCTION RNDM (ISEED) |
| 10 | C |
| 11 | C CERN PROGLIB# V104 RNDM .VERSION KERNMPW 1.00 891208 |
| 12 | C ORIG. 22-MAR-88 from IBM version of G.Marsaglia et al., Montreal |
| 13 | C |
| 14 | C- Uniform Random Number Generator, |
| 15 | C- giving the same sequence as the IBM and VAX version |
| 16 | |
| 17 | EQUIVALENCE (AMAN,MANT) |
| 18 | PARAMETER (MSK1 = 201326592) |
| 19 | PARAMETER (MSK2 = 855638016) |
| 20 | C PARAMETER (MSK1 = '0C000000'X) |
| 21 | C PARAMETER (MSK2 = '33000000'X) |
| 22 | COMMON / V104RN / MCGN |
| 23 | DATA MCGN /12345/ |
| 24 | |
| 25 | |
| 26 | MCGN = MCGN * 69069 |
| 27 | MANT = ISHFT (MCGN,-8) |
| 28 | IF (MANT.EQ.0) GO TO 14 |
| 29 | AMAN = MANT |
| 30 | C- AMAN in the range 1 to 2**24-1 |
| 31 | MANT = MANT - MSK1 |
| 32 | C- multiply by 2.**(-24) |
| 33 | RNDM = AMAN |
| 34 | RETURN |
| 35 | |
| 36 | C-- for zero set RNDM = 2.**(-25) |
| 37 | 14 MANT = MSK2 |
| 38 | RNDM = AMAN |
| 39 | RETURN |
| 40 | |
| 41 | C-- Integer random number |
| 42 | ENTRY IRNDM (ISEED) |
| 43 | MCGN = MCGN * 69069 |
| 44 | IRNDM = ISHFT (MCGN,-1) |
| 45 | END |
| 46 | |
| 47 | SUBROUTINE RDMIN (ISEED) |
| 48 | C-- Set the seed |
| 49 | COMMON / V104RN / MCGN |
| 50 | MCGN = ISEED |
| 51 | RETURN |
| 52 | |
| 53 | C-- Get the seed |
| 54 | ENTRY RDMOUT (ISEED) |
| 55 | ISEED = MCGN |
| 56 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame