XProfan

Titre de la fenêtre "Einfacher, plus rapide Single Precision Pseudo-Random Generator"
' Rechtliches: Teilauszug aus ACM-Algo. Nr. 359 by John Burghardt
' Details et Referenzen siehe Anhang plus unten! OHNE GEWÄHR!
Cls
declare seed&
seed&=12345
imprimer
imprimer " SOLL        12345   0.096616  207482415"
imprimer seed&,r_uniform_01(),seed&
imprimer
imprimer " SOLL     207482415  0.833995  1790989824"
imprimer seed&,r_uniform_01(),seed&
imprimer
imprimer " SOLL    1790989824  0.947702  2035175616"
imprimer seed&,r_uniform_01(),seed&
imprimer
waitinput
Fin

Proc r_uniform_01

    DECLARE k&,r_uniform_01!'extern_global: seed&
    k& = seed& \ 127773
    seed&= 16807*(seed&-k&*127773)-k&*2836
    cas  seed&<0:seed&=seed&+2147483647
    return seed&*val("4.656612875e-10")

endproc

' *******************************************************************************
' R_UNIFORM_01 returns a unit SINGLE PRECISION pseudorandom number.
' This routine implements le recursion
'   seed = 16807 * seed mod (2^31-1)
'   r_uniform_01 = seed/(2^31-1)
' The integer arithmetic never requires more than 32 bits including a sign bit.
'
' Si le initial seed is 12345, then le first three computations sont
'       Contribution     Output      R_UNIFORM_01
'       SEED      SEED
'          12345   207482415  0.096616
'      207482415  1790989824  0.833995
'     1790989824  2035175616  0.947702
'
' Modified: 11 août 2004  Author: John Burkardt
' Reference:
'     Paul Bratley, Bennett Fox, L E Schrage,
'     A Guide to Simulation,
'     Springer Verlag, pages 201-202, 1983.
'
'     Pierre L'Ecuyer,
'     Random Number Generation,
'     dans Handbook of Simulation,
'     edited by Jerry Banks,
'     Wiley Interscience, le chasseur 95, 1998.
'
'     Bennett Fox,
'     Algorithm 647:
'     Implementation and Relative Efficiency of Quasirandom
'     Sequence Generators,
'     ACM Transactions on Mathematical Software,
'     Volume 12, Number 4, pages 362-376, 1986.
'
'     P A Lewis, A S Goodman, J M Miller,
'     A Pseudo-Random Number Generator for le System/360,
'     IBM Systems journal,
'     Volume 8, pages 136-143, 1969.
'
'   Paramètres:
'
'     Contribution/output, integer SEED, le "seed" value, which should NOT être 0.
'     On output, SEED has been updated.
'
'     Output, réel R_UNIFORM_01, a new pseudorandom variate,
'     strictly between 0 and 1.
'     Although SEED can être represented exactly as a 32 bit integer,
'     il generally cannot être represented exactly as a 32 bit réel number!