stdlib/tst-rand48.c

   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 int
  12 main (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 }