d9/df6/setant_8f_source.html

       SUBROUTINE setant(qvalue)

 C**********************************************************************

 C

 C      SUBROUTINE SETANT(QVALUE)

 C               SET ANTithetic

 C

 C     Sets whether the current generator produces antithetic values.  If

 C     X   is  the value  normally returned  from  a uniform [0,1] random

 C     number generator then 1  - X is the antithetic  value. If X is the

 C     value  normally  returned  from a   uniform  [0,N]  random  number

 C     generator then N - 1 - X is the antithetic value.

 C

 C     All generators are initialized to NOT generate antithetic values.

 C

 C     This is a transcription from Pascal to Fortran of routine

 C     Set_Antithetic from the paper

 C

 C     L'Ecuyer, P. and Cote, S. "Implementing a Random Number Package

 C     with Splitting Facilities." ACM Transactions on Mathematical

 C     Software, 17:98-111 (1991)

 C

 C

 C                              Arguments

 C

 C

 C     QVALUE -> .TRUE. if generator G is to generating antithetic

 C                    values, otherwise .FALSE.

 C                                   LOGICAL QVALUE

 C

 C**********************************************************************

 C     .. Parameters ..

       INTEGER numg

       parameter(numg=32)

 C     ..

 C     .. Scalar Arguments ..

       LOGICAL qvalue

 C     ..

 C     .. Scalars in Common ..

       INTEGER a1,a1vw,a1w,a2,a2vw,a2w,m1,m2

 C     ..

 C     .. Arrays in Common ..

       INTEGER cg1(numg),cg2(numg),ig1(numg),ig2(numg),lg1(numg),

      +        lg2(numg)

       LOGICAL qanti(numg)

 C     ..

 C     .. Local Scalars ..

       INTEGER g

 C     ..

 C     .. External Functions ..

       LOGICAL qrgnin

       EXTERNAL qrgnin

 C     ..

 C     .. External Subroutines ..

       EXTERNAL getcgn

 C     ..

 C     .. Common blocks ..

       COMMON /globe/m1,m2,a1,a2,a1w,a2w,a1vw,a2vw,ig1,ig2,lg1,lg2,cg1,

      +       cg2,qanti

 C     ..

 C     .. Save statement ..

       SAVE /globe/

 C     ..

 C     .. Executable Statements ..

 C     Abort unless random number generator initialized

       IF (qrgnin()) go to 10

       WRITE (*,*) ' SETANT called before random number generator ',

      +  ' initialized -- abort!'

       CALL xstopx

      + (' SETANT called before random number generator initialized')


    10 CALL getcgn(g)

       qanti(g) = qvalue

       RETURN


       END

getcgn
subroutine getcgn(g)
Definition: getcgn.f:1

to
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
Definition: Quad-opts.cc:233

setant
subroutine setant(qvalue)
Definition: setant.f:1