home *** CD-ROM | disk | FTP | other *** search
/ The CDPD Public Domain Collection for CDTV 3 / CDPDIII.bin / pd / programming / gnusmalltalk / random.st < prev    next >
Text File  |  1992-02-15  |  3KB  |  120 lines
  1. "======================================================================
  2. |
  3. |   Random number Method Definitions
  4. |
  5.  ======================================================================"
  6.  
  7.  
  8. "======================================================================
  9. |
  10. | Copyright (C) 1990, 1991, 1992 Free Software Foundation, Inc.
  11. | Written by Steve Byrne.
  12. |
  13. | This file is part of GNU Smalltalk.
  14. |
  15. | GNU Smalltalk is free software; you can redistribute it and/or modify it
  16. | under the terms of the GNU General Public License as published by the Free
  17. | Software Foundation; either version 1, or (at your option) any later version.
  18. | 
  19. | GNU Smalltalk is distributed in the hope that it will be useful, but WITHOUT
  20. | ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  21. | FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  22. | details.
  23. | 
  24. | You should have received a copy of the GNU General Public License along with
  25. | GNU Smalltalk; see the file COPYING.  If not, write to the Free Software
  26. | Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.  
  27. |
  28.  ======================================================================"
  29.  
  30.  
  31. "
  32. |     Change Log
  33. | ============================================================================
  34. | Author       Date       Change 
  35. | sbb         16 Mar 91      Class creation now separate statement.
  36. |
  37. | sbyrne     19 Sep 89      Converted to use real method categories.
  38. |
  39. | sbyrne      3 Jul 89      created.
  40. |
  41. "
  42.  
  43. Stream subclass: #Random
  44.        instanceVariableNames: 'seed'
  45.        classVariableNames: ''
  46.        poolDictionaries: ''
  47.        category: nil
  48. !
  49.  
  50. Random comment: "Here's a random comment :-)"
  51. 'My instances are generator streams that produce random numbers, which are 
  52. floating point values between 0 and 1.'!
  53.  
  54. !Random class methodsFor: 'instance creation'!
  55.  
  56. new
  57.     ^self basicNew setSeed
  58. !!
  59.  
  60.  
  61.  
  62. !Random methodsFor: 'testing'!
  63.  
  64. chiSquare
  65.     "returns on Sun3 93.40000000000009"
  66.     ^self chiSquare: 1000 range: 100
  67. !
  68.  
  69. chiSquare: n range: r
  70.     | f t seed |
  71.     seed _ 1234567.
  72.     f _ Array new: r + 1.
  73.     1 to: r + 1 do: [ :i | f at: i put: 0 ].
  74.     n timesRepeat:
  75.     [ seed _ (seed * 31415821) + 1 bitAnd: 16r3FFFFFFF.
  76.           t _ seed \\ r.
  77.       f at: t + 1 put: (f at: t + 1) + 1 ].
  78.     t _ 0.
  79.     1 to: r do: [ :i | t _ t + (f at: i) squared ].
  80.     ^r asFloat * t / n - n
  81.  
  82. !!
  83.  
  84.  
  85. !Random methodsFor: 'basic'!
  86.  
  87. atEnd
  88.     ^false
  89. !
  90.  
  91. next
  92.     | value |
  93.     "From Sedgewick's 'Algorithms', based on Lehmer's method"
  94.     seed _ (seed * 31415821) + 1 bitAnd: 16r3FFFFFFF.
  95.     ^seed / 16r3FFFFFFF.0
  96. !
  97.  
  98. nextPut: value
  99.     self shouldNotImplement
  100. !
  101.  
  102. next: anInteger
  103.     | collection |
  104.     collection _ OrderedCollection new.
  105.     anInteger timesRepeat: [ collection add: self next ]. 
  106.     ^collection
  107. !
  108.  
  109. nextMatchFor: aNumber
  110.     ^self next = aNumber
  111. !!
  112.  
  113.  
  114.  
  115. !Random methodsFor: 'private'!
  116.  
  117. setSeed
  118.     seed _ Time secondClock
  119. !!
  120.