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

Required Arguments

MTABLE — Integer array of length 625 containing the table used in the 32-bit Mersenne Twister generator. (Output)

FORTRAN 90 Interface

Generic:                              CALL RNGE32 (MTABLE)

Specific:                             The specific interface name is RNGE32

FORTRAN 77 Interface

Single:                                CALL RNGE32 (MTABLE)

Description

The values in the table contain the state of the 32-bit Mersenne Twister random number generator. The table can be used by RNSE32 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 RNIN32_INT

      USE RNGE32_INT

      USE RNSET_INT

      USE UMACH_INT

      USE RNUN_INT

      IMPLICIT   NONE

      INTEGER    I, ISEED, NOUT

      INTEGER INIT(4)

      DATA INIT/291,564,837,1110/

      DATA ISEED/123457/

      INTEGER NR

      REAL R(5)

      INTEGER MTABLE(625)

      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(8)

      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 RNIN32(INIT)

!     SAVE THE STATE OF THE SERIES

      CALL RNGE32(MTABLE)

      CALL RNUN(R)

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

                   FMT, RLABEL, CLABEL)

!     RESTORE THE STATE OF THE TABLE

      CALL RNSE32(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 RNSE32(MTABLE)

      CALL RNUN(R)

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

                   FMT, RLABEL, CLABEL)

      END

Output

 

                     First stream output

     0.4347      0.3522      0.0139      0.2091      0.4956

                    Second stream output

     0.2486      0.2226      0.1111      0.9563      0.9846

                     Third stream output

     0.2486      0.2226      0.1111      0.9563      0.9846

                    Fourth stream output

     0.4347      0.3522      0.0139      0.2091      0.4956

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