home *** CD-ROM | disk | FTP | other *** search
/ DP Tool Club 8 / CDASC08.ISO / NEWS / 554 / JUIN / RANDOM.PAS < prev    next >
Pascal/Delphi Source File  |  1993-10-07  |  2KB  |  60 lines
  1. {─ Fido Pascal Conference ────────────────────────────────────────────── PASCAL ─
  2. Msg  : 333 of 389                                                               
  3. From : Erik Johnson                        1:104/28.0           15 Jun 93  18:43 
  4. To   : Ben Glazer                                                                
  5. Subj : RANDOM NUMBER GENERATOR                                                
  6. ────────────────────────────────────────────────────────────────────────────────
  7. > EJ> The random number generator in Pascal is not very 
  8. > EJ> good.  I have one that
  9. > EJ> takes 3 seeds and will not repeat for some trillion or so numbers.
  10. > EJ> can supply you with the code if you are interested.
  11. > 
  12. >I'd like to see this one, Erik!  Please, post away...
  13.  
  14. Here is the one I started with (as copied out of some long forgotten
  15. programming journal).  Unfortunately it is uncommented, and we are
  16. advised by the original columnist that modying the results can cripple
  17. the 'randomness' unless one knows they are doing.  This produces
  18. something like 'billions' of random numbers w/o repeating.  I have since
  19. modified (object oriented, longints...) but this version is smaller and
  20. simpler.}
  21.  
  22. unit random;
  23. {$D-}
  24.  
  25. interface
  26.  
  27. type
  28.     seedtype=record
  29.       seed1,seed2,seed3:word
  30.     end;
  31.  
  32. function ran(range:word):integer;
  33. procedure getseed(var s1,s2,s3:word);
  34. procedure setseed(s1,s2,s3:word);
  35. procedure incseed(num:longint);
  36.  
  37. var randomseed:seedtype;
  38.  
  39.  
  40. implementation
  41.  
  42. function ran(range:word):integer;
  43. begin
  44.   with randomseed do begin
  45.     seed1:=seed1*179 mod 32771;
  46.     seed2:=seed2*183 mod 32779;
  47.     seed3:=seed3*182 mod 32783
  48.   end;
  49.   ran:=0;
  50.   if range>1 then
  51. ran:=(randomseed.seed1+randomseed.seed2+randomseed.seed3) mod range
  52. end;
  53.  
  54. ** The rest got garbled, but they aren't crucial.  Procs for setting and
  55. getting the 3 current seed values, and incrementing the current seed
  56. values without getting a random number back.
  57.  
  58. I'm not sure how the numerical constants in the RAN function are
  59. derived, but I don't recommend changing them, unless you know how
  60. they're derived.