lh | 9ed821d | 2023-04-07 01:36:19 -0700 | [diff] [blame^] | 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 | |
| 11 | static int |
| 12 | do_test (void) |
| 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 (); |
| 47 | if (l != -0x5d73effdl) |
| 48 | { |
| 49 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", |
| 50 | __LINE__ - 4, -0x5d73effdl, l); |
| 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 (); |
| 63 | if (l != -0x61770b8cl) |
| 64 | { |
| 65 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", |
| 66 | __LINE__ - 4, -0x61770b8cl, l); |
| 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 (); |
| 95 | if (l != -0x1485e05dl) |
| 96 | { |
| 97 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", |
| 98 | __LINE__ - 4, -0x1485e05dl, l); |
| 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 (); |
| 174 | if (l != -0x5d73effdl) |
| 175 | { |
| 176 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", |
| 177 | __LINE__ - 4, -0x5d73effdl, l); |
| 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 (); |
| 190 | if (l != -0x61770b8cl) |
| 191 | { |
| 192 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", |
| 193 | __LINE__ - 4, -0x61770b8cl, l); |
| 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 (); |
| 234 | if (l != -0x1485e05dl) |
| 235 | { |
| 236 | printf ("mrand48() in line %d failed: expected %lx, seen %lx\n", |
| 237 | __LINE__ - 4, -0x1485e05dl, l); |
| 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); |
| 290 | if (l != -0xa973860l) |
| 291 | { |
| 292 | printf ("jrand48() in line %d failed: expected %lx, seen %lx\n", |
| 293 | __LINE__ - 4, -0xa973860l, l); |
| 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 | } |
| 391 | |
| 392 | #define TEST_FUNCTION do_test () |
| 393 | #include "../test-skeleton.c" |