Last change
on this file since 2403 was 2403, checked in by cmv, 22 years ago |
Creation du module de code source de MINUIT (CERNLIB) extrait par CMV
cmv 11/06/2003
|
File size:
912 bytes
|
Rev | Line | |
---|
[2403] | 1 | *
|
---|
| 2 | * $Id: mnrn15.F,v 1.1.1.1 2003-06-11 14:18:29 cmv Exp $
|
---|
| 3 | *
|
---|
| 4 | * $Log: not supported by cvs2svn $
|
---|
| 5 | * Revision 1.1.1.1 1996/03/07 14:31:31 mclareni
|
---|
| 6 | * Minuit
|
---|
| 7 | *
|
---|
| 8 | *
|
---|
| 9 | #include "minuit/pilot.h"
|
---|
| 10 | SUBROUTINE MNRN15(VAL,INSEED)
|
---|
| 11 | #include "minuit/d506dp.inc"
|
---|
| 12 | C This is a super-portable random number generator.
|
---|
| 13 | C It should not overflow on any 32-bit machine.
|
---|
| 14 | C The cycle is only ~10**9, so use with care!
|
---|
| 15 | C Note especially that VAL must not be undefined on input.
|
---|
| 16 | C Set Default Starting Seed
|
---|
| 17 | PARAMETER (THREE=3.0)
|
---|
| 18 | DATA ISEED/12345/
|
---|
| 19 | IF (VAL .EQ. THREE) GO TO 100
|
---|
| 20 | C
|
---|
| 21 | INSEED = ISEED
|
---|
| 22 | K = ISEED/53668
|
---|
| 23 | ISEED = 40014*(ISEED-K*53668) - K*12211
|
---|
| 24 | IF (ISEED .LT. 0) ISEED = ISEED + 2147483563
|
---|
| 25 | VAL = REAL(ISEED) * 4.656613E-10
|
---|
| 26 | RETURN
|
---|
| 27 | C "entry" to set seed, flag is VAL=3.
|
---|
| 28 | 100 ISEED = INSEED
|
---|
| 29 | RETURN
|
---|
| 30 | END
|
---|
Note:
See
TracBrowser
for help on using the repository browser.