*-- Author : F. James, modified by Mike Seymour
*- Split in 3 files by M. Kirsanov. Initial seeds ISEED set in HWUDAT
C-----------------------------------------------------------------------
- FUNCTION HWRGEN(I)
+
+C FUNCTION HWRGEN(I)
C-----------------------------------------------------------------------
C MAIN RANDOM NUMBER GENERATOR
C USES METHOD OF l'Ecuyer, (VIA F.JAMES, COMP PHYS COMM 60(1990)329)
C-----------------------------------------------------------------------
- IMPLICIT NONE
- DOUBLE PRECISION HWRGEN
- COMMON/HWSEED/ISEED(2)
- INTEGER ISEED
- INTEGER I,K,IZ
-C
- K=ISEED(1)/53668
- ISEED(1)=40014*(ISEED(1)-K*53668)-K*12211
- IF (ISEED(1).LT.0) ISEED(1)=ISEED(1)+2147483563
- K=ISEED(2)/52774
- ISEED(2)=40692*(ISEED(2)-K*52774)-K*3791
- IF (ISEED(2).LT.0) ISEED(2)=ISEED(2)+2147483399
- IZ=ISEED(1)-ISEED(2)
- IF (IZ.LT.1) IZ=IZ+2147483562
- HWRGEN=DBLE(IZ)*4.656613001013252D-10
+C IMPLICIT NONE
+C DOUBLE PRECISION HWRGEN
+C COMMON/HWSEED/ISEED(2)
+C INTEGER ISEED
+C INTEGER I,K,IZ
+C
+C K=ISEED(1)/53668
+C ISEED(1)=40014*(ISEED(1)-K*53668)-K*12211
+C IF (ISEED(1).LT.0) ISEED(1)=ISEED(1)+2147483563
+C K=ISEED(2)/52774
+C ISEED(2)=40692*(ISEED(2)-K*52774)-K*3791
+C IF (ISEED(2).LT.0) ISEED(2)=ISEED(2)+2147483399
+C IZ=ISEED(1)-ISEED(2)
+C IF (IZ.LT.1) IZ=IZ+2147483562
+C HWRGEN=DBLE(IZ)*4.656613001013252D-10
C---> (4.656613001013252D-10 = 1.D0/2147483589)
- END
+c END
CDECK ID>, HWRSET.
*CMZ :- -26/04/91 12.42.30 by Federico Carminati
*-- Author : F. James, modified by Mike Seymour