tstbot.for

Go to the documentation of this file.

      PROGRAM tstbot
C**********************************************************************
C
C     A test program for the bottom level routines
C
C**********************************************************************
C     Set up the random number generator
C     .. Local Scalars ..
      INTEGER ians,iblock,igen,iseed1,iseed2,itmp,ix,ixgen,nbad
C     ..
C     .. Local Arrays ..
      INTEGER answer(10000),genlst(5)
C     ..
C     .. External Functions ..
      INTEGER ignlgi
      EXTERNAL ignlgi
C     ..
C     .. External Subroutines ..
      EXTERNAL getsd,initgn,setall,setcgn
C     ..
C     .. Data statements ..
      DATA genlst/1,5,10,20,32/
C     ..
C     .. Executable Statements ..
      nbad = 0
      WRITE (*,9000)
 
 9000 FORMAT (' For five virual generators of the 32'/
     +       ' This test generates 10000 numbers then resets the block'/
     +       '      and does it again'/
     +       ' Any disagreements are reported -- there should be none'/)
C
C     Set up Generators
C
      CALL setall(12345,54321)
C
C     For a selected set of generators
C
      DO 60,ixgen = 1,5
          igen = genlst(ixgen)
          CALL setcgn(igen)
          WRITE (*,*) ' Testing generator ',igen
C
C     Use 10 blocks
C
          CALL initgn(-1)
          CALL getsd(iseed1,iseed2)
          DO 20,iblock = 1,10
C
C     Generate 1000 numbers
C
              DO 10,ians = 1,1000
                  ix = ians + (iblock-1)*1000
                  answer(ix) = ignlgi()
   10         CONTINUE
              CALL initgn(+1)
   20     CONTINUE
          CALL initgn(-1)
C
C     Do it again and compare answers
C
          CALL getsd(iseed1,iseed2)
C
C     Use 10 blocks
C
          DO 50,iblock = 1,10
C
C     Generate 1000 numbers
C
              DO 40,ians = 1,1000
                  ix = ians + (iblock-1)*1000
C      ANSWER( IX ) = IGNLGI()
                  itmp = ignlgi()
                  IF (.NOT. (itmp.NE.answer(ix))) GO TO 30
                  WRITE (*,9010) iblock,ians,ix,answer(ix),itmp
 
 9010             FORMAT (' Disagreement on regeneration of numbers'/
     +                   ' Block ',i2,' N within Block ',i2,
     +                   ' Index in answer ',i5/
     +                   ' Originally Generated ',i10,' Regenerated ',
     +                   i10)
 
                  nbad = nbad + 1
                  IF (nbad.GT.10) stop ' More than 10 mismatches'
   30             CONTINUE
   40         CONTINUE
              CALL initgn(+1)
   50     CONTINUE
          WRITE (*,*) ' Finished testing generator ',igen
          WRITE (*,*) ' Test completed successfully'
   60 CONTINUE
      stop
 
      END