Retrieves the current table used in the 64-bit Mersenne Twister generator.

Required Arguments

MTABLE — Integer(kind=8) array of length 313 containing the table used in the 64-bit Mersenne Twister  generator. (Output)

FORTRAN 90 Interface

Generic:                              CALL RNGE64 (MTABLE)

Specific:                             The specific interface name is RNGE64

FORTRAN 77 Interface

Single:                                CALL RNGE64 (MTABLE)

Description

The values in the table contain the state of the 64-bit Mersenne Twister random number generator. The table can be used by RNSE64 to set the generator back to this state.

Example

In this example, four simulation streams are generated. The first series is generated with the seed used for initialization. The second series is generated using an array for initialization. The third series is obtained by resetting the generator back to the state it had at the beginning of the second stream. Therefore, the second and third streams are identical. The fourth stream is obtained by resetting the generator back to its original, uninitialized state, and having it reinitialize using the seed. The first and fourth streams are therefore the same.

 

      USE RNIN64_INT

      USE RNGE64_INT

      USE RNSET_INT

      USE UMACH_INT

      USE RNUN_INT

      IMPLICIT   NONE

      INTEGER    I, ISEED, NOUT

      INTEGER(KIND=8) INIT(4)

      DATA INIT/291,564,837,1110/

      DATA ISEED/123457/

      INTEGER NR

      REAL R(5)

      INTEGER(KIND=8) MTABLE(313)

      CHARACTER CLABEL(5)*5, FMT*8, RLABEL(3)*5

      RLABEL(1)='NONE'

      CLABEL(1)='NONE'

      DATA FMT/'(W10.4)'/

      NR=5

      CALL UMACH (2, NOUT)

      ISEED = 123457

      CALL RNOPT(9)

      CALL RNSET(ISEED)

      CALL RNUN(R)

      CALL WRRRL('FIRST STREAM OUTPUT',1,5,R,1,0, &

                  FMT, RLABEL, CLABEL)

!     REINITIALIZE MERSENNE TWISTER SERIES WITH AN ARRAY

      CALL RNIN64(INIT)

!     SAVE THE STATE OF THE SERIES

      CALL RNGE64(MTABLE)

      CALL RNUN(R)

      CALL WRRRL('SECOND STREAM OUTPUT',1,5,R,1,0, &

                   FMT, RLABEL, CLABEL)

!     RESTORE THE STATE OF THE TABLE

      CALL RNSE64(MTABLE)

      CALL RNUN(R)

      CALL WRRRL('THIRD STREAM OUTPUT',1,5,R,1,0, &

                   FMT, RLABEL, CLABEL)

!     RESET THE SERIES - IT WILL REINITIALIZE FROM THE SEED

      MTABLE(1)=1000

      CALL RNSE64(MTABLE)

      CALL RNUN(R)

      CALL WRRRL('FOURTH STREAM OUTPUT',1,5,R,1,0, &

                   FMT, RLABEL, CLABEL)

      END

Output

 

                     First stream output

     0.5799      0.9401      0.7102      0.1640      0.5457

                    Second stream output

     0.4894      0.7397      0.5725      0.0863      0.7588

                     Third stream output

     0.4894      0.7397      0.5725      0.0863      0.7588

                    Fourth stream output

     0.5799      0.9401      0.7102      0.1640      0.5457

Visual Numerics, Inc. Visual Numerics - Developers of IMSL and PV-WAVE http://www.vni.com/ PHONE: 713.784.3131 FAX:713.781.9260