]>
Commit | Line | Data |
---|---|---|
df152cc8 UD |
1 | #include <float.h> |
2 | #include <stdio.h> | |
3 | #include <stdlib.h> | |
4 | #include <string.h> | |
5 | ||
6 | #ifndef DECIMAL_DIG | |
7 | # define DECIMAL_DIG 21 | |
8 | #endif | |
9 | ||
10 | ||
29955b5d AS |
11 | static int |
12 | do_test (void) | |
df152cc8 UD |
13 | { |
14 | unsigned short int xs[3] = { 0x0001, 0x0012, 0x0123 }; | |
15 | unsigned short int lxs[7]; | |
16 | unsigned short int *xsp; | |
17 | int result = 0; | |
18 | long int l; | |
19 | double d; | |
20 | double e; | |
21 | ||
22 | /* Test srand48. */ | |
23 | srand48 (0x98765432); | |
24 | /* Get the values of the internal Xi array. */ | |
25 | xsp = seed48 (xs); | |
26 | if (xsp[0] != 0x330e || xsp[1] != 0x5432 || xsp[2] != 0x9876) | |
27 | { | |
28 | puts ("srand48(0x98765432) didn't set correct value"); | |
29 | printf (" expected: { %04hx, %04hx, %04hx }\n", 0x330e, 0x5432, 0x9876); | |
30 | printf (" seen: { %04hx, %04hx, %04hx }\n", xsp[0], xsp[1], xsp[2]); | |
31 | result = 1; | |
32 | } | |
33 | /* Put the values back. */ | |
34 | memcpy (xs, xsp, sizeof (xs)); | |
35 | (void) seed48 (xs); | |
36 | ||
37 | /* See whether the correct values are installed. */ | |
38 | l = lrand48 (); | |
39 | if (l != 0x2fed1413l) | |
40 | { | |
41 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
42 | __LINE__ - 4, 0x2fed1413l, l); | |
43 | result = 1; | |
44 | } | |
45 | ||
46 | l = mrand48 (); | |
5bcba0ed | 47 | if (l != -0x5d73effdl) |
df152cc8 UD |
48 | { |
49 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 50 | __LINE__ - 4, -0x5d73effdl, l); |
df152cc8 UD |
51 | result = 1; |
52 | } | |
53 | ||
54 | l = lrand48 (); | |
55 | if (l != 0x585fcfb7l) | |
56 | { | |
57 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
58 | __LINE__ - 4, 0x585fcfb7l, l); | |
59 | result = 1; | |
60 | } | |
61 | ||
62 | l = mrand48 (); | |
5bcba0ed | 63 | if (l != -0x61770b8cl) |
df152cc8 UD |
64 | { |
65 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 66 | __LINE__ - 4, -0x61770b8cl, l); |
df152cc8 UD |
67 | result = 1; |
68 | } | |
69 | ||
70 | /* Test seed48. The previous call should have install the values in | |
71 | the initialization of `xs' above. */ | |
72 | xs[0] = 0x1234; | |
73 | xs[1] = 0x5678; | |
74 | xs[2] = 0x9012; | |
75 | xsp = seed48 (xs); | |
76 | if (xsp[0] != 0x62f2 || xsp[1] != 0xf474 || xsp[2] != 0x9e88) | |
77 | { | |
78 | puts ("seed48() did not install the values correctly"); | |
79 | printf (" expected: { %04hx, %04hx, %04hx }\n", 0x62f2, 0xf474, 0x9e88); | |
80 | printf (" seen: { %04hx, %04hx, %04hx }\n", xsp[0], xsp[1], xsp[2]); | |
81 | result = 1; | |
82 | } | |
83 | ||
84 | /* Test lrand48 and mrand48. We continue from the seed established | |
85 | above. */ | |
86 | l = lrand48 (); | |
87 | if (l != 0x017e48b5l) | |
88 | { | |
89 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
90 | __LINE__ - 4, 0x017e48b5l, l); | |
91 | result = 1; | |
92 | } | |
93 | ||
94 | l = mrand48 (); | |
5bcba0ed | 95 | if (l != -0x1485e05dl) |
df152cc8 UD |
96 | { |
97 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 98 | __LINE__ - 4, -0x1485e05dl, l); |
df152cc8 UD |
99 | result = 1; |
100 | } | |
101 | ||
102 | l = lrand48 (); | |
103 | if (l != 0x6b6a3f95l) | |
104 | { | |
105 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
106 | __LINE__ - 4, 0x6b6a3f95l, l); | |
107 | result = 1; | |
108 | } | |
109 | ||
110 | l = mrand48 (); | |
111 | if (l != 0x175c0d6fl) | |
112 | { | |
113 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
114 | __LINE__ - 4, 0x175c0d6fl, l); | |
115 | result = 1; | |
116 | } | |
117 | ||
118 | /* Test lcong48. */ | |
119 | lxs[0] = 0x4567; | |
120 | lxs[1] = 0x6789; | |
121 | lxs[2] = 0x8901; | |
122 | lxs[3] = 0x0123; | |
123 | lxs[4] = 0x2345; | |
124 | lxs[5] = 0x1111; | |
125 | lxs[6] = 0x2222; | |
126 | lcong48 (lxs); | |
127 | ||
128 | /* See whether the correct values are installed. */ | |
129 | l = lrand48 (); | |
130 | if (l != 0x6df63d66l) | |
131 | { | |
132 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
133 | __LINE__ - 4, 0x6df63d66l, l); | |
134 | result = 1; | |
135 | } | |
136 | ||
137 | l = mrand48 (); | |
138 | if (l != 0x2f92c8e1l) | |
139 | { | |
140 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
141 | __LINE__ - 4, 0x2f92c8e1l, l); | |
142 | result = 1; | |
143 | } | |
144 | ||
145 | l = lrand48 (); | |
146 | if (l != 0x3b4869ffl) | |
147 | { | |
148 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
149 | __LINE__ - 4, 0x3b4869ffl, l); | |
150 | result = 1; | |
151 | } | |
152 | ||
153 | l = mrand48 (); | |
154 | if (l != 0x5cd4cc3el) | |
155 | { | |
156 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
157 | __LINE__ - 4, 0x5cd4cc3el, l); | |
158 | result = 1; | |
159 | } | |
160 | ||
161 | /* Check whether srand48() restores the A and C parameters. */ | |
162 | srand48 (0x98765432); | |
163 | ||
164 | /* See whether the correct values are installed. */ | |
165 | l = lrand48 (); | |
166 | if (l != 0x2fed1413l) | |
167 | { | |
168 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
169 | __LINE__ - 4, 0x2fed1413l, l); | |
170 | result = 1; | |
171 | } | |
172 | ||
173 | l = mrand48 (); | |
5bcba0ed | 174 | if (l != -0x5d73effdl) |
df152cc8 UD |
175 | { |
176 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 177 | __LINE__ - 4, -0x5d73effdl, l); |
df152cc8 UD |
178 | result = 1; |
179 | } | |
180 | ||
181 | l = lrand48 (); | |
182 | if (l != 0x585fcfb7l) | |
183 | { | |
184 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
185 | __LINE__ - 4, 0x585fcfb7l, l); | |
186 | result = 1; | |
187 | } | |
188 | ||
189 | l = mrand48 (); | |
5bcba0ed | 190 | if (l != -0x61770b8cl) |
df152cc8 UD |
191 | { |
192 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 193 | __LINE__ - 4, -0x61770b8cl, l); |
df152cc8 UD |
194 | result = 1; |
195 | } | |
196 | ||
197 | /* And again to see whether seed48() does the same. */ | |
198 | lcong48 (lxs); | |
199 | ||
200 | /* See whether lxs wasn't modified. */ | |
201 | l = lrand48 (); | |
202 | if (l != 0x6df63d66l) | |
203 | { | |
204 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
205 | __LINE__ - 4, 0x6df63d66l, l); | |
206 | result = 1; | |
207 | } | |
208 | ||
209 | /* Test seed48. The previous call should have install the values in | |
210 | the initialization of `xs' above. */ | |
211 | xs[0] = 0x1234; | |
212 | xs[1] = 0x5678; | |
213 | xs[2] = 0x9012; | |
214 | xsp = seed48 (xs); | |
215 | if (xsp[0] != 0x0637 || xsp[1] != 0x7acd || xsp[2] != 0xdbec) | |
216 | { | |
217 | puts ("seed48() did not install the values correctly"); | |
218 | printf (" expected: { %04hx, %04hx, %04hx }\n", 0x0637, 0x7acd, 0xdbec); | |
219 | printf (" seen: { %04hx, %04hx, %04hx }\n", xsp[0], xsp[1], xsp[2]); | |
220 | result = 1; | |
221 | } | |
222 | ||
223 | /* Test lrand48 and mrand48. We continue from the seed established | |
224 | above. */ | |
225 | l = lrand48 (); | |
226 | if (l != 0x017e48b5l) | |
227 | { | |
228 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
229 | __LINE__ - 4, 0x017e48b5l, l); | |
230 | result = 1; | |
231 | } | |
232 | ||
233 | l = mrand48 (); | |
5bcba0ed | 234 | if (l != -0x1485e05dl) |
df152cc8 UD |
235 | { |
236 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 237 | __LINE__ - 4, -0x1485e05dl, l); |
df152cc8 UD |
238 | result = 1; |
239 | } | |
240 | ||
241 | l = lrand48 (); | |
242 | if (l != 0x6b6a3f95l) | |
243 | { | |
244 | printf ("lrand48() in line %d failed: expected %lx, seen %lx\n", | |
245 | __LINE__ - 4, 0x6b6a3f95l, l); | |
246 | result = 1; | |
247 | } | |
248 | ||
249 | l = mrand48 (); | |
250 | if (l != 0x175c0d6fl) | |
251 | { | |
252 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", | |
253 | __LINE__ - 4, 0x175c0d6fl, l); | |
254 | result = 1; | |
255 | } | |
256 | ||
257 | /* Test drand48. */ | |
258 | d = drand48 (); | |
259 | if (d != 0.0908832261858485424) | |
260 | { | |
261 | printf ("drand48() in line %d failed: expected %.*g, seen %.*g\n", | |
262 | __LINE__ - 4, DECIMAL_DIG, 0.0908832261858485424, | |
263 | DECIMAL_DIG, d); | |
264 | result = 1; | |
265 | } | |
266 | ||
267 | d = drand48 (); | |
268 | if (d != 0.943149381730059133133) | |
269 | { | |
270 | printf ("drand48() in line %d failed: expected %.*g, seen %.*g\n", | |
271 | __LINE__ - 4, DECIMAL_DIG, 0.943149381730059133133, | |
272 | DECIMAL_DIG, d); | |
273 | result = 1; | |
274 | } | |
275 | ||
276 | /* Now the functions which get the Xis passed. */ | |
277 | xs[0] = 0x3849; | |
278 | xs[1] = 0x5061; | |
279 | xs[2] = 0x7283; | |
280 | ||
281 | l = nrand48 (xs); | |
282 | if (l != 0x1efe61a1l) | |
283 | { | |
284 | printf ("nrand48() in line %d failed: expected %lx, seen %lx\n", | |
285 | __LINE__ - 4, 0x1efe61a1l, l); | |
286 | result = 1; | |
287 | } | |
288 | ||
289 | l = jrand48 (xs); | |
5bcba0ed | 290 | if (l != -0xa973860l) |
df152cc8 UD |
291 | { |
292 | printf ("jrand48() in line %d failed: expected %lx, seen %lx\n", | |
5bcba0ed | 293 | __LINE__ - 4, -0xa973860l, l); |
df152cc8 UD |
294 | result = 1; |
295 | } | |
296 | ||
297 | l = nrand48 (xs); | |
298 | if (l != 0x2a5e57fel) | |
299 | { | |
300 | printf ("nrand48() in line %d failed: expected %lx, seen %lx\n", | |
301 | __LINE__ - 4, 0x2a5e57fel, l); | |
302 | result = 1; | |
303 | } | |
304 | ||
305 | l = jrand48 (xs); | |
306 | if (l != 0x71a779a8l) | |
307 | { | |
308 | printf ("jrand48() in line %d failed: expected %lx, seen %lx\n", | |
309 | __LINE__ - 4, 0x71a779a8l, l); | |
310 | result = 1; | |
311 | } | |
312 | ||
313 | /* Test whether the global A and C are used. */ | |
314 | lcong48 (lxs); | |
315 | ||
316 | l = nrand48 (xs); | |
317 | if (l != 0x32beee9fl) | |
318 | { | |
319 | printf ("nrand48() in line %d failed: expected %lx, seen %lx\n", | |
320 | __LINE__ - 4, 0x32beee9fl, l); | |
321 | result = 1; | |
322 | } | |
323 | ||
324 | l = jrand48 (xs); | |
325 | if (l != 0x7bddf3bal) | |
326 | { | |
327 | printf ("jrand48() in line %d failed: expected %lx, seen %lx\n", | |
328 | __LINE__ - 4, 0x7bddf3bal, l); | |
329 | result = 1; | |
330 | } | |
331 | ||
332 | l = nrand48 (xs); | |
333 | if (l != 0x85bdf28l) | |
334 | { | |
335 | printf ("nrand48() in line %d failed: expected %lx, seen %lx\n", | |
336 | __LINE__ - 4, 0x85bdf28l, l); | |
337 | result = 1; | |
338 | } | |
339 | ||
340 | l = jrand48 (xs); | |
341 | if (l != 0x7b433e47l) | |
342 | { | |
343 | printf ("jrand48() in line %d failed: expected %lx, seen %lx\n", | |
344 | __LINE__ - 4, 0x7b433e47l, l); | |
345 | result = 1; | |
346 | } | |
347 | ||
348 | /* Test erand48. Also compare with the drand48 results. */ | |
349 | (void) seed48 (xs); | |
350 | ||
351 | d = drand48 (); | |
352 | e = erand48 (xs); | |
353 | if (d != e) | |
354 | { | |
355 | printf ("\ | |
356 | drand48() and erand48 in lines %d and %d produce different results\n", | |
357 | __LINE__ - 6, __LINE__ - 5); | |
358 | printf (" drand48() = %g, erand48() = %g\n", d, e); | |
359 | result = 1; | |
360 | } | |
361 | else if (e != 0.640650904452755298735) | |
362 | { | |
363 | printf ("erand48() in line %d failed: expected %.*g, seen %.*g\n", | |
364 | __LINE__ - 4, DECIMAL_DIG, 0.640650904452755298735, | |
365 | DECIMAL_DIG, e); | |
366 | result = 1; | |
367 | ||
368 | } | |
369 | ||
370 | d = drand48 (); | |
371 | e = erand48 (xs); | |
372 | if (d != e) | |
373 | { | |
374 | printf ("\ | |
375 | drand48() and erand48 in lines %d and %d produce different results\n", | |
376 | __LINE__ - 6, __LINE__ - 5); | |
377 | printf (" drand48() = %g, erand48() = %g\n", d, e); | |
378 | result = 1; | |
379 | } | |
380 | else if (e != 0.115372323508150742555) | |
381 | { | |
382 | printf ("erand48() in line %d failed: expected %.*g, seen %.*g\n", | |
383 | __LINE__ - 4, DECIMAL_DIG, 0.0115372323508150742555, | |
384 | DECIMAL_DIG, e); | |
385 | result = 1; | |
386 | ||
387 | } | |
388 | ||
389 | return result; | |
390 | } | |
29955b5d AS |
391 | |
392 | #define TEST_FUNCTION do_test () | |
393 | #include "../test-skeleton.c" |