home *** CD-ROM | disk | FTP | other *** search
/ Microsoft Programmer's Library 1.3 / Microsoft-Programers-Library-v1.3.iso / sampcode / fortran / demoran.for < prev    next >
Encoding:
Text File  |  1988-08-11  |  2.9 KB  |  109 lines
  1.     PROGRAM DEMORAN
  2. c
  3. c  This program demonstrates a uniform pseudo-random number
  4. c  generator. 
  5. c
  6. c  Portions of this program are based on ideas presented in the
  7. c  the book "Numerical Recipes - The Art of Scientific Computing"
  8. c  by William H. Press, Brian P. Flannery, Saul A. Teukolsky, and 
  9. c  William T. Vetterling, Cambridge University Press 1986
  10. c
  11. c  If 1000 numbers and 10 bins are requested, each bin should
  12. c  (ideally) be filled with 100 numbers. The percentage error
  13. c  is printed for each bin.
  14. c
  15. c  The following routines are provided:
  16. c
  17. c  REAL FUNCTION RAND () 
  18. c    returns a real number in the range 0. to 1.0
  19. c 
  20. c  INTEGER FUNCTION RANDLIM (ILO,IHI)
  21. c    returns a random integer in the range ILO to IHI
  22. c
  23. c  REAL FUNCTION SRAND (SEED)
  24. c    initializes either generator (seed = 0. to 259199.)
  25. c
  26. c  SUBROUTINE SECOND (TX)
  27. c    returns the number of seconds and hundreths of seconds elapsed
  28. c    since midnight
  29. c
  30. c
  31. c  NOTE -- Both generators should produce identical results
  32. c
  33.     INTEGER BINS(0:999), RANDLIM
  34.     ERR(I) = (I-FLOAT(NREP/NBINS))/(NREP/NBINS)*100.
  35.     WRITE (*,*) 'You will be asked to provide the following:'
  36.     WRITE (*,*) 'how many random numbers to generate'
  37.     WRITE (*,*) 'how many bins to use (1-1000)'
  38.     WRITE (*,*) 'which generator to use (1 or 2)'
  39.     WRITE (*,*)
  40. 10    CONTINUE
  41.     WRITE (*,*) 'Input three numbers separated by blanks or commas'
  42.     WRITE (*,*) 'or CTRL-Z to end'
  43.     READ (*,*,END=999) NREP,NBINS,IGEN
  44.     SEED = SRAND(1.0)
  45.     CALL SECOND (T1)
  46.     DO 100 I=1,NREP
  47.       IF (IGEN .EQ. 1) THEN
  48.             IX = RANDLIM(0,NBINS-1)
  49.       ELSE
  50.         IX = NBINS*RAND()
  51.       ENDIF
  52.       BINS(IX) = BINS(IX)+1
  53. 100    CONTINUE
  54.     CALL SECOND (T2)
  55.     WRITE (*,*) 'Time elapsed=',t2-t1
  56.     WRITE (*,*) 'Numbers generated per second=',nrep/(t2-t1)
  57.     WRITE (*,*)
  58.     WRITE (*,*) 'Bin    Count  % Error'
  59.     WRITE (*,*) '----  ------- -------'
  60.     DO 200 I=0,NBINS-1
  61.       WRITE (*,'(1x,i4,i9,f7.1,''%'')') i+1,bins(i),err(bins(i))
  62.           BINS(I) = 0
  63. 200    CONTINUE
  64.     GO TO 10
  65. 999    CONTINUE
  66.     END
  67.     FUNCTION RANDOM ()
  68. c
  69. c  If called, RANDOM just returns 0.0
  70. c
  71.     INTEGER RANDLIM
  72.     PARAMETER (IA=7141, IC=54773, IM=259200)
  73.     RANDOM = 0.0
  74.     RETURN
  75. c
  76. c  REAL FUNCTION SRAND (SEED)
  77. c    initializes either generator (seed = 0. to 259199.)
  78. c
  79.     ENTRY SRAND (X)
  80.     SRAND = X
  81.         SEED = X
  82.     RETURN
  83. c 
  84. c  INTEGER FUNCTION RANDLIM (ILO,IHI)
  85. c    returns a random integer in the range ILO to IHI
  86. c
  87.     ENTRY RANDLIM (ILO,IHI)
  88.     SEED = MOD (INT(SEED)*IA+IC,IM)
  89.     RANDLIM = ILO+(IHI-ILO+1)*SEED/IM
  90.     RETURN
  91. c
  92. c  REAL FUNCTION RAND () 
  93. c    returns a real number in the range 0. to 1.0
  94. c 
  95.     ENTRY RAND ()
  96.     SEED = MOD (INT(SEED)*IA+IC,IM)
  97.     RAND = SEED/IM
  98.     END
  99.     SUBROUTINE SECOND (TX)
  100. c
  101. c  SUBROUTINE SECOND (TX)
  102. c    returns the number of seconds and hundredths of seconds elapsed
  103. c    since midnight
  104. c
  105.     INTEGER*2 IH,IM,IS,IHU
  106.     CALL GETTIM (IH,IM,IS,IHU)
  107.     TX = IH*3600.+IM*60+IS+IHU/100.
  108.     END
  109.