home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #19 / NN_1992_19.iso / spool / comp / lang / c / 12930 < prev    next >
Encoding:
Text File  |  1992-08-29  |  2.2 KB  |  71 lines
  1. Path: sparky!uunet!gatech!purdue!yuma!longs.LANCE.ColoState.EDU!tp718390
  2. From: tp718390@longs.LANCE.ColoState.EDU (Thomas Proell)
  3. Newsgroups: comp.lang.c
  4. Subject: Re: random number generation
  5. Message-ID: <Aug28.205915.65889@yuma.ACNS.ColoState.EDU>
  6. Date: 28 Aug 92 20:59:15 GMT
  7. References: <1992Aug27.133805.21292@u.washington.edu>
  8. Sender: news@yuma.ACNS.ColoState.EDU (News Account)
  9. Reply-To: tp718390@longs.LANCE.ColoState.EDU (Thomas Proell)
  10. Organization: Engineering College, Colorado State University
  11. Lines: 58
  12.  
  13. In article <1992Aug27.133805.21292@u.washington.edu>,
  14. scipio@milton.u.washington.edu (jim lee) writes:
  15. > Can someone help me?  I need some coding that will generate random
  16. > numbers.
  17. > 
  18. > Thanks for any help.
  19. > 
  20. > 
  21. > Jim Lee (scipio@u.washington.edu)
  22.  
  23. A very elegant way to generate many different streams of uniformly
  24. distributed pseudo-random numbers is to use the ?rand48() library.
  25. The ? stands for either one of the following letters [d, e, l, n,
  26. m, j]. This assumes that you have access to a good C-compiler. 
  27.  
  28. drand48(), lrand48(), and mrand48() require to be initialized, i.e.
  29. srand48(time(NULL)); and are called without an argument. 
  30.  
  31. erand48(x sub i), nrand48(x sub i), and jrand48(x sub i) do not need
  32. to be initialized. x stands for an >unsigned short< integer array. Using
  33. different arrays as arguments for calling the subroutine gives you the 
  34. ability to generate different streams of i.e. floating point values; i.e: 
  35.  
  36. main(int argc, char **argv)
  37. {
  38.  
  39.   unsigned short *x, *y, ....
  40.  
  41.   float          *a, *b, ....
  42.  
  43.   x = calloc(MAX, sizeof(unsigned short)); 
  44.   y = .......
  45.   a = calloc(MAX, sizeof(float)); 
  46.   .....
  47.  
  48. /* Seed the first element. All subsequent values are generated 
  49.    by the function.                                               */
  50.  
  51.   x[0] = 0; 
  52.   y[0] = 10; 
  53.  
  54.   START LOOP
  55.       a[loop] = (float) nrand48(x); 
  56.       b[loop] = (float) nrand48(y); 
  57.   .
  58.   .
  59.  
  60. lrand48() and nrand48() gives a non-negative long integer in [0.0....2^31)
  61. mrand48() and jrand48() gives signed long integer in [-2^31....2^31]
  62. drand48() and erand48() gives a non-negative double precision in [0.0...1.0]
  63. {Had troubles with the last one!}
  64.  
  65.  
  66. Hope I could help, 
  67. Thomas Proell
  68. Department of Chemical Engineering
  69. Colorado State University
  70. Fort Collins, CO 80525  
  71.