1
SUBROUTINE getsd(iseed1,iseed2)
2
C**********************************************************************
4
C SUBROUTINE GETSD(G,ISEED1,ISEED2)
7
C Returns the value of two integer seeds of the current generator
9
C This is a transcription from Pascal to Fortran of routine
10
C Get_State from the paper
12
C L'Ecuyer, P. and Cote, S. "Implementing a Random Number Package
13
C with Splitting Facilities." ACM Transactions on Mathematical
14
C Software, 17:98-111 (1991)
21
C ISEED1 <- First integer seed of generator G
24
C ISEED2 <- Second integer seed of generator G
27
C**********************************************************************
33
C .. Scalar Arguments ..
36
C .. Scalars in Common ..
37
INTEGER a1,a1vw,a1w,a2,a2vw,a2w,m1,m2
39
C .. Arrays in Common ..
40
INTEGER cg1(numg),cg2(numg),ig1(numg),ig2(numg),lg1(numg),
47
C .. External Functions ..
51
C .. External Subroutines ..
55
COMMON /globe/m1,m2,a1,a2,a1w,a2w,a1vw,a2vw,ig1,ig2,lg1,lg2,cg1,
58
C .. Save statement ..
61
C .. Executable Statements ..
62
C Abort unless random number generator initialized
63
IF (qrgnin()) GO TO 10
64
call basout(io,wte,'GETSD called before random number generator ')