|
65
|
1 |
"
|
|
|
2 |
======================================================================
|
|
0
|
3 |
|
|
|
|
4 |
| Copyright (C) 1988, 1989 Free Software Foundation, Inc.
|
|
|
5 |
| Written by Steve Byrne.
|
|
|
6 |
|
|
|
|
7 |
| This file is part of GNU Smalltalk.
|
|
|
8 |
|
|
|
|
9 |
| GNU Smalltalk is free software; you can redistribute it and/or modify it
|
|
|
10 |
| under the terms of the GNU General Public License as published by the Free
|
|
|
11 |
| Software Foundation; either version 1, or (at your option) any later version.
|
|
|
12 |
|
|
|
|
13 |
| GNU Smalltalk is distributed in the hope that it will be useful, but WITHOUT
|
|
|
14 |
| ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
|
|
|
15 |
| FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
|
|
|
16 |
| details.
|
|
|
17 |
|
|
|
|
18 |
| You should have received a copy of the GNU General Public License along with
|
|
|
19 |
| GNU Smalltalk; see the file LICENSE. If not, write to the Free Software
|
|
|
20 |
| Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
|
|
|
21 |
|
|
|
65
|
22 |
======================================================================
|
|
0
|
23 |
|
|
65
|
24 |
see notice in (Random>>documentation)
|
|
0
|
25 |
"
|
|
|
26 |
|
|
|
27 |
Stream subclass:#Random
|
|
65
|
28 |
instanceVariableNames:'seed'
|
|
|
29 |
classVariableNames:''
|
|
|
30 |
poolDictionaries:''
|
|
|
31 |
category:'Magnitude-Numbers'
|
|
0
|
32 |
!
|
|
|
33 |
|
|
65
|
34 |
!Random class methodsFor:'documentation'!
|
|
0
|
35 |
|
|
31
|
36 |
copyright
|
|
|
37 |
"
|
|
|
38 |
======================================================================
|
|
|
39 |
|
|
|
|
40 |
| Copyright (C) 1988, 1989 Free Software Foundation, Inc.
|
|
|
41 |
| Written by Steve Byrne.
|
|
|
42 |
|
|
|
|
43 |
| This file is part of GNU Smalltalk.
|
|
|
44 |
|
|
|
|
45 |
| GNU Smalltalk is free software; you can redistribute it and/or modify it
|
|
|
46 |
| under the terms of the GNU General Public License as published by the Free
|
|
|
47 |
| Software Foundation; either version 1, or (at your option) any later version.
|
|
|
48 |
|
|
|
|
49 |
| GNU Smalltalk is distributed in the hope that it will be useful, but WITHOUT
|
|
|
50 |
| ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
|
|
|
51 |
| FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
|
|
|
52 |
| details.
|
|
|
53 |
|
|
|
|
54 |
| You should have received a copy of the GNU General Public License along with
|
|
|
55 |
| GNU Smalltalk; see the file LICENSE. If not, write to the Free Software
|
|
|
56 |
| Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
|
|
|
57 |
|
|
|
|
58 |
======================================================================
|
|
|
59 |
|
|
|
60 |
see notice in (Random>>documentation)
|
|
|
61 |
"
|
|
|
62 |
!
|
|
|
63 |
|
|
|
64 |
documentation
|
|
|
65 |
"
|
|
|
66 |
random numbers - thanks to Steves GNU Smalltalk
|
|
|
67 |
|
|
|
68 |
Notice: although being included here,
|
|
75
|
69 |
this file is NOT covered by the ST/X license, but by
|
|
|
70 |
the FSF copyLeft (see copyright method).
|
|
31
|
71 |
|
|
75
|
72 |
You can redistribute it under the terms stated there ...
|
|
|
73 |
Also, the price you pay for ST/X does not include a charge for
|
|
|
74 |
this file - it has to be considered as a separate piece of
|
|
|
75 |
software, which can be copied and given away without any
|
|
31
|
76 |
restriction from my (CG) side.
|
|
75
|
77 |
|
|
|
78 |
claus: the algorithm may need a rewrite for better numbers
|
|
31
|
79 |
"
|
|
|
80 |
! !
|
|
0
|
81 |
|
|
65
|
82 |
!Random class methodsFor:'instance creation'!
|
|
0
|
83 |
|
|
|
84 |
new
|
|
17
|
85 |
"return a new random generator"
|
|
|
86 |
|
|
0
|
87 |
^self basicNew setSeed
|
|
|
88 |
! !
|
|
|
89 |
|
|
65
|
90 |
!Random methodsFor:'accessing-reading'!
|
|
0
|
91 |
|
|
142
|
92 |
next
|
|
|
93 |
"return the next random number in the range 0..1"
|
|
17
|
94 |
|
|
|
95 |
seed := (seed times:31415821) + 1 bitAnd: 16r3FFFFFFF.
|
|
142
|
96 |
^ self nextInteger / 16r3FFFFFFF asFloat
|
|
17
|
97 |
|
|
|
98 |
"|r|
|
|
|
99 |
r := Random new.
|
|
142
|
100 |
Transcript showCr:r next.
|
|
|
101 |
Transcript showCr:r next.
|
|
|
102 |
Transcript showCr:r next.
|
|
|
103 |
Transcript showCr:r next.
|
|
17
|
104 |
"
|
|
0
|
105 |
!
|
|
|
106 |
|
|
65
|
107 |
nextBetween:start and:stop
|
|
|
108 |
"return a random between start and stop.
|
|
|
109 |
claus: the original GNU version has a bug in returning values
|
|
|
110 |
from the interval [start .. stop+1]"
|
|
|
111 |
|
|
|
112 |
|rnd|
|
|
|
113 |
|
|
|
114 |
rnd := self next.
|
|
|
115 |
rnd := rnd * (stop asFloat - start asFloat).
|
|
|
116 |
rnd := rnd + start asFloat.
|
|
|
117 |
^ rnd
|
|
|
118 |
|
|
|
119 |
"|r|
|
|
|
120 |
r := Random new.
|
|
|
121 |
Transcript showCr:(r nextBetween:1 and:10).
|
|
|
122 |
Transcript showCr:(r nextBetween:1 and:10).
|
|
|
123 |
Transcript showCr:(r nextBetween:1 and:10).
|
|
|
124 |
Transcript showCr:(r nextBetween:1 and:10).
|
|
|
125 |
"
|
|
|
126 |
!
|
|
|
127 |
|
|
0
|
128 |
nextBoolean
|
|
17
|
129 |
"return true or false by random"
|
|
|
130 |
|
|
0
|
131 |
"thanks to Peter Deutsch ..."
|
|
|
132 |
seed := (seed times:31415821) + 1 bitAnd: 16r3FFFFFFF.
|
|
|
133 |
^ seed < 16r20000000
|
|
17
|
134 |
|
|
|
135 |
"|r|
|
|
|
136 |
r := Random new.
|
|
|
137 |
Transcript showCr:r nextBoolean.
|
|
|
138 |
Transcript showCr:r nextBoolean.
|
|
|
139 |
Transcript showCr:r nextBoolean.
|
|
|
140 |
Transcript showCr:r nextBoolean.
|
|
|
141 |
"
|
|
0
|
142 |
!
|
|
|
143 |
|
|
142
|
144 |
nextInteger
|
|
|
145 |
"return the next integral random number,
|
|
|
146 |
in the range 0 .. 16r3FFFFFFF.
|
|
|
147 |
From Sedgewick's 'Algorithms', based on Lehmer's method"
|
|
|
148 |
|
|
|
149 |
"the times: is a kludge - times does not convert to LargeInteger on overflow"
|
|
|
150 |
seed := (seed times:31415821) + 1 bitAnd: 16r3FFFFFFF.
|
|
|
151 |
^ seed
|
|
|
152 |
|
|
|
153 |
"|r|
|
|
|
154 |
r := Random new.
|
|
|
155 |
Transcript showCr:r nextInteger.
|
|
|
156 |
Transcript showCr:r nextInteger.
|
|
|
157 |
Transcript showCr:r nextInteger.
|
|
|
158 |
Transcript showCr:r nextInteger.
|
|
|
159 |
"
|
|
|
160 |
!
|
|
|
161 |
|
|
17
|
162 |
nextIntegerBetween:start and:stop
|
|
|
163 |
"return an integral random between start and stop"
|
|
|
164 |
|
|
|
165 |
|rnd|
|
|
|
166 |
|
|
|
167 |
rnd := self next.
|
|
|
168 |
rnd := rnd * (stop asFloat - start asFloat + 1.0).
|
|
|
169 |
^ (rnd + start) truncated.
|
|
|
170 |
|
|
|
171 |
"|r|
|
|
|
172 |
r := Random new.
|
|
|
173 |
Transcript showCr:(r nextIntegerBetween:1 and:10).
|
|
|
174 |
Transcript showCr:(r nextIntegerBetween:1 and:10).
|
|
|
175 |
Transcript showCr:(r nextIntegerBetween:1 and:10).
|
|
|
176 |
Transcript showCr:(r nextIntegerBetween:1 and:10).
|
|
|
177 |
"
|
|
0
|
178 |
!
|
|
|
179 |
|
|
|
180 |
nextMatchFor: aNumber
|
|
17
|
181 |
"generate the next random, return true iff it has the same
|
|
|
182 |
value as aNumber. Redefined to avoid endless reading."
|
|
|
183 |
|
|
0
|
184 |
^self next = aNumber
|
|
|
185 |
! !
|
|
|
186 |
|
|
65
|
187 |
!Random methodsFor:'blocked methods'!
|
|
0
|
188 |
|
|
65
|
189 |
contents
|
|
|
190 |
"blocked from use - contents makes no sense for random generators"
|
|
|
191 |
|
|
|
192 |
self shouldNotImplement
|
|
|
193 |
!
|
|
|
194 |
|
|
|
195 |
nextPut: value
|
|
|
196 |
"blocked from use - it makes no sense for randoms"
|
|
|
197 |
|
|
|
198 |
self shouldNotImplement
|
|
|
199 |
! !
|
|
|
200 |
|
|
|
201 |
!Random methodsFor:'private'!
|
|
0
|
202 |
|
|
|
203 |
setSeed
|
|
17
|
204 |
"set the initial seed value based on the current time"
|
|
|
205 |
|
|
0
|
206 |
seed := Time secondClock bitAnd: 16r3FFFFFFF
|
|
|
207 |
! !
|
|
65
|
208 |
|
|
142
|
209 |
!Random methodsFor:'testing'!
|
|
|
210 |
|
|
|
211 |
atEnd
|
|
|
212 |
"instances of Random can always give more numbers"
|
|
|
213 |
|
|
|
214 |
^ false
|
|
|
215 |
!
|
|
|
216 |
|
|
|
217 |
chiSquare
|
|
|
218 |
"perform a chiSquare-test on the receiver"
|
|
|
219 |
|
|
|
220 |
"returns on Sun3 93.40000000000009"
|
|
|
221 |
^self chiSquare: 1000 range: 100
|
|
|
222 |
|
|
|
223 |
"Random new chiSquare"
|
|
|
224 |
!
|
|
|
225 |
|
|
|
226 |
chiSquare: n range: r
|
|
|
227 |
| f t s |
|
|
|
228 |
s := 1234567.
|
|
|
229 |
f := Array new: r + 1.
|
|
|
230 |
1 to: r + 1 do: [ :i | f at: i put: 0 ].
|
|
|
231 |
n timesRepeat:
|
|
|
232 |
[ s := (seed * 31415821) + 1 bitAnd: 16r3FFFFFFF.
|
|
|
233 |
t := s \\ r.
|
|
|
234 |
f at: t + 1 put: (f at: t + 1) + 1 ].
|
|
|
235 |
t := 0.
|
|
|
236 |
1 to: r do: [ :i | t := t + (f at: i) squared ].
|
|
|
237 |
^r asFloat * t / n - n
|
|
|
238 |
!
|
|
|
239 |
|
|
|
240 |
isReadable
|
|
|
241 |
^ true
|
|
|
242 |
! !
|
|
|
243 |
|
|
|
244 |
!Random class methodsFor:'documentation'!
|
|
|
245 |
|
|
|
246 |
version
|
|
|
247 |
^ '$Header: /cvs/stx/stx/libbasic2/Random.st,v 1.13 1995-12-07 21:36:14 cg Exp $'
|
|
|
248 |
! !
|