home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #18 / NN_1992_18.iso / spool / comp / programm / 2397 < prev    next >
Encoding:
Text File  |  1992-08-20  |  2.3 KB  |  76 lines
  1. Newsgroups: comp.programming
  2. Path: sparky!uunet!newsgate.watson.ibm.com!yktnews!admin!curt
  3. From: curt@watson.ibm.com (Curt McDowell)
  4. Subject: Re: Any simple Random Number generators?
  5. Sender: news@watson.ibm.com (NNTP News Poster)
  6. Message-ID: <1992Aug20.143244.6305@watson.ibm.com>
  7. Date: Thu, 20 Aug 1992 14:32:44 GMT
  8. Lines: 61
  9. Disclaimer: This posting represents the poster's views, not necessarily those of IBM
  10. References: <13AUG199218224504@judy.uh.edu> <1992Aug14.104121.29374@corax.udac.uu.se> <3283@creatures.cs.vt.edu> <1992Aug18.143439@fdr.jsc.nasa.gov>
  11. Nntp-Posting-Host: gorby.watson.ibm.com
  12. Organization: IBM T.J. Watson Research Center
  13. Keywords: random, shift
  14.  
  15. In article <1992Aug18.143439@fdr.jsc.nasa.gov>, shirley@fdr.jsc.nasa.gov (Bill Shirley [CSC]) writes:
  16. > 
  17. > If you just want it to reseed differently
  18. > 
  19. > random(time())  
  20. > 
  21. > will do it
  22. > -- 
  23. > Bill Shirley        |  ``Computer Science is not about computers any 
  24. > shirley@fdr.jsc.nasa.gov|    more than astronomy is about telescopes.''
  25.  
  26. No, that won't quite do it.  Firstly, the call is srandom.  Secondly,
  27. time() requires an argument or it will core dump.  What you wanted to
  28. say was srandom(time(0)).
  29.  
  30. A slightly better way is srandom(time(0) + getpid()), since if your
  31. program runs quick it might be run more than once a second (e.g. from a
  32. shell script).
  33.  
  34. Here are my favorite convenience #defines when my program is not in
  35. a hurry:
  36.  
  37. #define rander() (((unsigned long) random() << 5) ^ \
  38.            (unsigned long) random()       ^ \
  39.           ((unsigned long) random() >> 5))
  40.  
  41. #define srnd()   srandom(time(0) + getpid())
  42. #define rnd(n)   (rander() % (n))
  43.  
  44. Here is my favorite way to generate vaguely random numbers when my
  45. program is in a tremendous hurry (e.g. as for test data).  It's a
  46. pseudo-random shift register of the type detailed in a previous note
  47. from Andy Collins, cycling through every possible bit combination
  48. except 0 before recycling.  The numbers are far from ideally random,
  49. but it's very fast.  You may be able to find a way to make them
  50. appear more random.
  51.  
  52. /*
  53.  * randomize_buffer
  54.  *
  55.  *   Fill a buffer of 'len' unsigned longs with semi-random data.
  56.  */
  57.  
  58. void randomize_buffer(buf, len)
  59.     unsigned long  *buf;
  60.     int        len;
  61. {
  62.     unsigned long    b;
  63.  
  64.     while ((b = rander()) == 0)
  65.         ;
  66.  
  67.     while (len--) {
  68.         *buf++ = b;
  69.                 b = b << 1 | (b >> 30 ^ b >> 27) & 1;
  70.     }
  71. }
  72.  
  73. --
  74. Curt McDowell
  75. IBM T. J. Watson Research Center
  76.