home *** CD-ROM | disk | FTP | other *** search
/ PsL Monthly 1993 December / PSL Monthly Shareware CD-ROM (December 1993).iso / prgmming / dos / c / snippet.exe / RAND2.C < prev    next >
Encoding:
C/C++ Source or Header  |  1990-06-25  |  1.2 KB  |  46 lines
  1. /********************************************************************
  2.           The McGill Super-Duper Random Number Generator
  3.              G. Marsaglia, K. Ananthanarayana, N. Paul
  4.   
  5.     Incorporating the Ziggurat method of sampling from decreasing
  6.               or symmetric unimodal density functions.
  7.                       G. Marsaglia, W.W. Tsang
  8.   
  9.                  Rewritten into C by E. Schneider
  10.  *********************************************************************/
  11.   
  12.  static unsigned long mcgn, srgn;
  13.   
  14.  #define MULT 69069L
  15.   
  16.  void rstart (i1, i2)
  17.  long i1, i2;
  18.   
  19.  { mcgn = (i1 == 0) ? 0 : i1 | 1;
  20.    srgn = (i2 == 0) ? 0 : (i2 & 0x7FF) | 1;
  21.   }
  22.   
  23.  long uni()
  24.  { unsigned long r0, r1;
  25.    r0 = (srgn >> 15);
  26.    r1 = srgn ^ r0;
  27.    r0 = (r1 << 17);
  28.    srgn = r0 ^ r1;
  29.    mcgn = MULT * mcgn;
  30.    r1 = mcgn ^ srgn;
  31.    return (r1 >> 1);  }
  32.   
  33.  long vni()
  34.  { unsigned long r0, r1;
  35.    r0 = (srgn >> 15);
  36.    r1 = srgn ^ r0;
  37.    r0 = (r1 << 17);
  38.    srgn = r0 ^ r1;
  39.    mcgn = MULT * mcgn;
  40.    r1 = mcgn ^ srgn;
  41.    return r1;  }
  42. /* 
  43. "Anyone who consider arithmetic means of producing random number is,
  44.  of course, in a state of sin" - John Von Neumann
  45. */
  46.