math/test-misc.c

   1 /* Miscellaneous tests which don't fit anywhere else.
   2    Copyright (C) 2000, 2001 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4
   5    The GNU C Library is free software; you can redistribute it and/or
   6    modify it under the terms of the GNU Lesser General Public
   7    License as published by the Free Software Foundation; either
   8    version 2.1 of the License, or (at your option) any later version.
   9
  10    The GNU C Library is distributed in the hope that it will be useful,
  11    but WITHOUT ANY WARRANTY; without even the implied warranty of
  12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13    Lesser General Public License for more details.
  14
  15    You should have received a copy of the GNU Lesser General Public
  16    License along with the GNU C Library; if not, write to the Free
  17    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
  18    02111-1307 USA.  */
  19
  20 #include <fenv.h>
  21 #include <float.h>
  22 #include <ieee754.h>
  23 #include <math.h>
  24 #include <stdio.h>
  25 #include <string.h>
  26
  27
  28 int
  29 main (void)
  30 {
  31   int result = 0;
  32
  33 #ifndef NO_LONG_DOUBLE
  34   {
  35     long double x = 0x100000001ll + (long double) 0.5;
  36     long double q;
  37     long double r;
  38
  39     r = modfl (x, &q);
  40     if (q != (long double) 0x100000001ll || r != 0.5)
  41       {
  42         printf ("modfl (%Lg, ...) failed\n", x);
  43         result = 1;
  44       }
  45   }
  46
  47 # if __GNUC__ >= 3 || __GNUC_MINOR__ >= 96
  48   {
  49     long double x;
  50     long double m;
  51     long double r;
  52     int e;
  53     int i;
  54
  55 #  if LDBL_MANT_DIG == 64
  56     m = 0xf.fffffffffffffffp-4L;
  57 #  elif LDBL_MANT_DIG == 113
  58     m = 0x1.ffffffffffffffffffffffffffffp-1L;
  59 #  else
  60 #   error "Please adjust"
  61 #  endif
  62
  63     for (i = LDBL_MAX_EXP, x = LDBL_MAX; i >= LDBL_MIN_EXP; --i, x /= 2.0L)
  64       {
  65         printf ("2^%d: ", i);
  66
  67         r = frexpl (x, &e);
  68         if (r != m)
  69           {
  70             printf ("mantissa incorrect: %.20La\n", r);
  71             result = 1;
  72             continue;
  73           }
  74         if (e != i)
  75           {
  76             printf ("exponent wrong %d (%.20Lg)\n", e, x);
  77             result = 1;
  78             continue;
  79           }
  80         puts ("ok");
  81       }
  82
  83     for (i = LDBL_MIN_EXP, x = LDBL_MIN; i >= LDBL_MIN_EXP - LDBL_MANT_DIG + 1;
  84          --i, x /= 2.0L)
  85       {
  86         printf ("2^%d: ", i);
  87
  88         r = frexpl (x, &e);
  89         if (r != 0.5L)
  90           {
  91             printf ("mantissa incorrect: %.20La\n", r);
  92             result = 1;
  93             continue;
  94           }
  95         if (e != i)
  96           {
  97             printf ("exponent wrong %d (%.20Lg)\n", e, x);
  98             result = 1;
  99             continue;
 100           }
 101         puts ("ok");
 102       }
 103
 104   }
 105 # endif
 106
 107 #if 0
 108   {
 109     int e;
 110     long double r = frexpl (LDBL_MIN * LDBL_EPSILON, &e);
 111
 112     if (r != 0.5)
 113       {
 114         printf ("frexpl (LDBL_MIN * LDBL_EPSILON, ...): mantissa wrong: %Lg\n",
 115                 r);
 116         result = 1;
 117       }
 118     else if (e != -16444)
 119       {
 120         printf ("frexpl (LDBL_MIN * LDBL_EPSILON, ...): exponent wrong: %d\n",
 121                 e);
 122         result = 1;
 123       }
 124   }
 125 #endif
 126 #endif
 127
 128   {
 129     double x = 0x100000001ll + (double) 0.5;
 130     double q;
 131     double r;
 132
 133     r = modf (x, &q);
 134     if (q != (double) 0x100000001ll || r != 0.5)
 135       {
 136         printf ("modf (%g, ...) failed\n", x);
 137         result = 1;
 138       }
 139   }
 140
 141   {
 142     union ieee754_float v1;
 143     union ieee754_float v2;
 144     float f;
 145
 146     v1.f = f = FLT_MIN;
 147     if (fpclassify (f) != FP_NORMAL)
 148       {
 149         printf ("fpclassify (FLT_MIN) failed: %d\n", fpclassify (f));
 150         result = 1;
 151       }
 152     f = nextafterf (f, FLT_MIN / 2.0f);
 153     if (fpclassify (f) != FP_SUBNORMAL)
 154       {
 155         printf ("fpclassify (FLT_MIN-epsilon) failed: %d\n", fpclassify (f));
 156         result = 1;
 157       }
 158     v2.f = f = nextafterf (f, FLT_MIN);
 159     if (fpclassify (f) != FP_NORMAL)
 160       {
 161         printf ("fpclassify (FLT_MIN-epsilon+epsilon) failed: %d\n",
 162                 fpclassify (f));
 163         result = 1;
 164       }
 165
 166     if (v1.ieee.mantissa != v2.ieee.mantissa)
 167       {
 168         printf ("FLT_MIN: mantissa differs: %8x vs %8x\n",
 169                 v1.ieee.mantissa, v2.ieee.mantissa);
 170         result = 1;
 171       }
 172     if (v1.ieee.exponent != v2.ieee.exponent)
 173       {
 174         printf ("FLT_MIN: exponent differs: %4x vs %4x\n",
 175                 v1.ieee.exponent, v2.ieee.exponent);
 176         result = 1;
 177       }
 178     if (v1.ieee.negative != v2.ieee.negative)
 179       {
 180         printf ("FLT_MIN: negative differs: %d vs %d\n",
 181                 v1.ieee.negative, v2.ieee.negative);
 182         result = 1;
 183       }
 184
 185     v1.f = f = -FLT_MIN;
 186     if (fpclassify (f) != FP_NORMAL)
 187       {
 188         printf ("fpclassify (-FLT_MIN) failed: %d\n", fpclassify (f));
 189         result = 1;
 190       }
 191     f = nextafterf (f, -FLT_MIN / 2.0f);
 192     if (fpclassify (f) != FP_SUBNORMAL)
 193       {
 194         printf ("fpclassify (-FLT_MIN-epsilon) failed: %d\n", fpclassify (f));
 195         result = 1;
 196       }
 197     v2.f = f = nextafterf (f, -FLT_MIN);
 198     if (fpclassify (f) != FP_NORMAL)
 199       {
 200         printf ("fpclassify (-FLT_MIN-epsilon+epsilon) failed: %d\n",
 201                 fpclassify (f));
 202         result = 1;
 203       }
 204
 205     if (v1.ieee.mantissa != v2.ieee.mantissa)
 206       {
 207         printf ("-FLT_MIN: mantissa differs: %8x vs %8x\n",
 208                 v1.ieee.mantissa, v2.ieee.mantissa);
 209         result = 1;
 210       }
 211     if (v1.ieee.exponent != v2.ieee.exponent)
 212       {
 213         printf ("-FLT_MIN: exponent differs: %4x vs %4x\n",
 214                 v1.ieee.exponent, v2.ieee.exponent);
 215         result = 1;
 216       }
 217     if (v1.ieee.negative != v2.ieee.negative)
 218       {
 219         printf ("-FLT_MIN: negative differs: %d vs %d\n",
 220                 v1.ieee.negative, v2.ieee.negative);
 221         result = 1;
 222       }
 223
 224     f = FLT_MAX;
 225     if (fpclassify (f) != FP_NORMAL)
 226       {
 227         printf ("fpclassify (FLT_MAX) failed: %d\n", fpclassify (f));
 228         result = 1;
 229       }
 230     f = nextafterf (f, INFINITY);
 231     if (fpclassify (f) != FP_INFINITE)
 232       {
 233         printf ("fpclassify (FLT_MAX+epsilon) failed: %d\n", fpclassify (f));
 234         result = 1;
 235       }
 236
 237     f = -FLT_MAX;
 238     if (fpclassify (f) != FP_NORMAL)
 239       {
 240         printf ("fpclassify (-FLT_MAX) failed: %d\n", fpclassify (f));
 241         result = 1;
 242       }
 243     f = nextafterf (f, -INFINITY);
 244     if (fpclassify (f) != FP_INFINITE)
 245       {
 246         printf ("fpclassify (-FLT_MAX-epsilon) failed: %d\n", fpclassify (f));
 247         result = 1;
 248       }
 249
 250     v1.f = f = 0.0625;
 251     f = nextafterf (f, 0.0);
 252     v2.f = f = nextafterf (f, 1.0);
 253
 254     if (v1.ieee.mantissa != v2.ieee.mantissa)
 255       {
 256         printf ("0.0625f down: mantissa differs: %8x vs %8x\n",
 257                 v1.ieee.mantissa, v2.ieee.mantissa);
 258         result = 1;
 259       }
 260     if (v1.ieee.exponent != v2.ieee.exponent)
 261       {
 262         printf ("0.0625f down: exponent differs: %4x vs %4x\n",
 263                 v1.ieee.exponent, v2.ieee.exponent);
 264         result = 1;
 265       }
 266     if (v1.ieee.negative != v2.ieee.negative)
 267       {
 268         printf ("0.0625f down: negative differs: %d vs %d\n",
 269                 v1.ieee.negative, v2.ieee.negative);
 270         result = 1;
 271       }
 272
 273     v1.f = f = 0.0625;
 274     f = nextafterf (f, 1.0);
 275     v2.f = f = nextafterf (f, 0.0);
 276
 277     if (v1.ieee.mantissa != v2.ieee.mantissa)
 278       {
 279         printf ("0.0625f up: mantissa differs: %8x vs %8x\n",
 280                 v1.ieee.mantissa, v2.ieee.mantissa);
 281         result = 1;
 282       }
 283     if (v1.ieee.exponent != v2.ieee.exponent)
 284       {
 285         printf ("0.0625f up: exponent differs: %4x vs %4x\n",
 286                 v1.ieee.exponent, v2.ieee.exponent);
 287         result = 1;
 288       }
 289     if (v1.ieee.negative != v2.ieee.negative)
 290       {
 291         printf ("0.0625f up: negative differs: %d vs %d\n",
 292                 v1.ieee.negative, v2.ieee.negative);
 293         result = 1;
 294       }
 295
 296     v1.f = f = -0.0625;
 297     f = nextafterf (f, 0.0);
 298     v2.f = f = nextafterf (f, -1.0);
 299
 300     if (v1.ieee.mantissa != v2.ieee.mantissa)
 301       {
 302         printf ("-0.0625f up: mantissa differs: %8x vs %8x\n",
 303                 v1.ieee.mantissa, v2.ieee.mantissa);
 304         result = 1;
 305       }
 306     if (v1.ieee.exponent != v2.ieee.exponent)
 307       {
 308         printf ("-0.0625f up: exponent differs: %4x vs %4x\n",
 309                 v1.ieee.exponent, v2.ieee.exponent);
 310         result = 1;
 311       }
 312     if (v1.ieee.negative != v2.ieee.negative)
 313       {
 314         printf ("-0.0625f up: negative differs: %d vs %d\n",
 315                 v1.ieee.negative, v2.ieee.negative);
 316         result = 1;
 317       }
 318
 319     v1.f = f = -0.0625;
 320     f = nextafterf (f, -1.0);
 321     v2.f = f = nextafterf (f, 0.0);
 322
 323     if (v1.ieee.mantissa != v2.ieee.mantissa)
 324       {
 325         printf ("-0.0625f down: mantissa differs: %8x vs %8x\n",
 326                 v1.ieee.mantissa, v2.ieee.mantissa);
 327         result = 1;
 328       }
 329     if (v1.ieee.exponent != v2.ieee.exponent)
 330       {
 331         printf ("-0.0625f down: exponent differs: %4x vs %4x\n",
 332                 v1.ieee.exponent, v2.ieee.exponent);
 333         result = 1;
 334       }
 335     if (v1.ieee.negative != v2.ieee.negative)
 336       {
 337         printf ("-0.0625f down: negative differs: %d vs %d\n",
 338                 v1.ieee.negative, v2.ieee.negative);
 339         result = 1;
 340       }
 341
 342     v1.f = f = 0.0f;
 343     f = nextafterf (f, 1.0);
 344     v2.f = nextafterf (f, -1.0);
 345
 346     if (v1.ieee.mantissa != v2.ieee.mantissa)
 347       {
 348         printf ("0.0f up: mantissa differs: %8x vs %8x\n",
 349                 v1.ieee.mantissa, v2.ieee.mantissa);
 350         result = 1;
 351       }
 352     if (v1.ieee.exponent != v2.ieee.exponent)
 353       {
 354         printf ("0.0f up: exponent differs: %4x vs %4x\n",
 355                 v1.ieee.exponent, v2.ieee.exponent);
 356         result = 1;
 357       }
 358     if (0 != v2.ieee.negative)
 359       {
 360         printf ("0.0f up: negative differs: 0 vs %d\n",
 361                 v2.ieee.negative);
 362         result = 1;
 363       }
 364
 365     v1.f = f = 0.0f;
 366     f = nextafterf (f, -1.0);
 367     v2.f = nextafterf (f, 1.0);
 368
 369     if (v1.ieee.mantissa != v2.ieee.mantissa)
 370       {
 371         printf ("0.0f down: mantissa differs: %8x vs %8x\n",
 372                 v1.ieee.mantissa, v2.ieee.mantissa);
 373         result = 1;
 374       }
 375     if (v1.ieee.exponent != v2.ieee.exponent)
 376       {
 377         printf ("0.0f down: exponent differs: %4x vs %4x\n",
 378                 v1.ieee.exponent, v2.ieee.exponent);
 379         result = 1;
 380       }
 381     if (1 != v2.ieee.negative)
 382       {
 383         printf ("0.0f down: negative differs: 1 vs %d\n",
 384                 v2.ieee.negative);
 385         result = 1;
 386       }
 387
 388     if (nextafterf (0.0f, INFINITY) != nextafterf (0.0f, 1.0f)
 389         || nextafterf (-0.0f, INFINITY) != nextafterf (-0.0f, 1.0f)
 390         || nextafterf (0.0f, -INFINITY) != nextafterf (0.0f, -1.0f)
 391         || nextafterf (-0.0f, -INFINITY) != nextafterf (-0.0f, -1.0f))
 392       {
 393         printf ("nextafterf (+-0, +-Inf) != nextafterf (+-0, +-1)\n");
 394         result = 1;
 395       }
 396
 397     if (nexttowardf (0.0f, INFINITY) != nexttowardf (0.0f, 1.0f)
 398         || nexttowardf (-0.0f, INFINITY) != nexttowardf (-0.0f, 1.0f)
 399         || nexttowardf (0.0f, -INFINITY) != nexttowardf (0.0f, -1.0f)
 400         || nexttowardf (-0.0f, -INFINITY) != nexttowardf (-0.0f, -1.0f))
 401       {
 402         printf ("nexttowardf (+-0, +-Inf) != nexttowardf (+-0, +-1)\n");
 403         result = 1;
 404       }
 405   }
 406
 407   {
 408     union ieee754_double v1;
 409     union ieee754_double v2;
 410     double d;
 411
 412     v1.d = d = DBL_MIN;
 413     if (fpclassify (d) != FP_NORMAL)
 414       {
 415         printf ("fpclassify (DBL_MIN) failed: %d\n", fpclassify (d));
 416         result = 1;
 417       }
 418     d = nextafter (d, DBL_MIN / 2.0);
 419     if (fpclassify (d) != FP_SUBNORMAL)
 420       {
 421         printf ("fpclassify (DBL_MIN-epsilon) failed: %d\n", fpclassify (d));
 422         result = 1;
 423       }
 424     v2.d = d = nextafter (d, DBL_MIN);
 425     if (fpclassify (d) != FP_NORMAL)
 426       {
 427         printf ("fpclassify (DBL_MIN-epsilon+epsilon) failed: %d\n",
 428                 fpclassify (d));
 429         result = 1;
 430       }
 431
 432     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 433       {
 434         printf ("DBL_MIN: mantissa0 differs: %8x vs %8x\n",
 435                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 436         result = 1;
 437       }
 438     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 439       {
 440         printf ("DBL_MIN: mantissa1 differs: %8x vs %8x\n",
 441                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 442         result = 1;
 443       }
 444     if (v1.ieee.exponent != v2.ieee.exponent)
 445       {
 446         printf ("DBL_MIN: exponent differs: %4x vs %4x\n",
 447                 v1.ieee.exponent, v2.ieee.exponent);
 448         result = 1;
 449       }
 450     if (v1.ieee.negative != v2.ieee.negative)
 451       {
 452         printf ("DBL_MIN: negative differs: %d vs %d\n",
 453                 v1.ieee.negative, v2.ieee.negative);
 454         result = 1;
 455       }
 456
 457     v1.d = d = -DBL_MIN;
 458     if (fpclassify (d) != FP_NORMAL)
 459       {
 460         printf ("fpclassify (-DBL_MIN) failed: %d\n", fpclassify (d));
 461         result = 1;
 462       }
 463     d = nextafter (d, -DBL_MIN / 2.0);
 464     if (fpclassify (d) != FP_SUBNORMAL)
 465       {
 466         printf ("fpclassify (-DBL_MIN-epsilon) failed: %d\n", fpclassify (d));
 467         result = 1;
 468       }
 469     v2.d = d = nextafter (d, -DBL_MIN);
 470     if (fpclassify (d) != FP_NORMAL)
 471       {
 472         printf ("fpclassify (-DBL_MIN-epsilon+epsilon) failed: %d\n",
 473                 fpclassify (d));
 474         result = 1;
 475       }
 476
 477     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 478       {
 479         printf ("-DBL_MIN: mantissa0 differs: %8x vs %8x\n",
 480                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 481         result = 1;
 482       }
 483     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 484       {
 485         printf ("-DBL_MIN: mantissa1 differs: %8x vs %8x\n",
 486                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 487         result = 1;
 488       }
 489     if (v1.ieee.exponent != v2.ieee.exponent)
 490       {
 491         printf ("-DBL_MIN: exponent differs: %4x vs %4x\n",
 492                 v1.ieee.exponent, v2.ieee.exponent);
 493         result = 1;
 494       }
 495     if (v1.ieee.negative != v2.ieee.negative)
 496       {
 497         printf ("-DBL_MIN: negative differs: %d vs %d\n",
 498                 v1.ieee.negative, v2.ieee.negative);
 499         result = 1;
 500       }
 501
 502     d = DBL_MAX;
 503     if (fpclassify (d) != FP_NORMAL)
 504       {
 505         printf ("fpclassify (DBL_MAX) failed: %d\n", fpclassify (d));
 506         result = 1;
 507       }
 508     d = nextafter (d, INFINITY);
 509     if (fpclassify (d) != FP_INFINITE)
 510       {
 511         printf ("fpclassify (DBL_MAX+epsilon) failed: %d\n", fpclassify (d));
 512         result = 1;
 513       }
 514
 515     d = -DBL_MAX;
 516     if (fpclassify (d) != FP_NORMAL)
 517       {
 518         printf ("fpclassify (-DBL_MAX) failed: %d\n", fpclassify (d));
 519         result = 1;
 520       }
 521     d = nextafter (d, -INFINITY);
 522     if (fpclassify (d) != FP_INFINITE)
 523       {
 524         printf ("fpclassify (-DBL_MAX-epsilon) failed: %d\n", fpclassify (d));
 525         result = 1;
 526       }
 527
 528     v1.d = d = 0.0625;
 529     d = nextafter (d, 0.0);
 530     v2.d = d = nextafter (d, 1.0);
 531
 532     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 533       {
 534         printf ("0.0625 down: mantissa0 differs: %8x vs %8x\n",
 535                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 536         result = 1;
 537       }
 538     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 539       {
 540         printf ("0.0625 down: mantissa1 differs: %8x vs %8x\n",
 541                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 542         result = 1;
 543       }
 544     if (v1.ieee.exponent != v2.ieee.exponent)
 545       {
 546         printf ("0.0625 down: exponent differs: %4x vs %4x\n",
 547                 v1.ieee.exponent, v2.ieee.exponent);
 548         result = 1;
 549       }
 550     if (v1.ieee.negative != v2.ieee.negative)
 551       {
 552         printf ("0.0625 down: negative differs: %d vs %d\n",
 553                 v1.ieee.negative, v2.ieee.negative);
 554         result = 1;
 555       }
 556
 557     v1.d = d = 0.0625;
 558     d = nextafter (d, 1.0);
 559     v2.d = d = nextafter (d, 0.0);
 560
 561     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 562       {
 563         printf ("0.0625 up: mantissa0 differs: %8x vs %8x\n",
 564                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 565         result = 1;
 566       }
 567     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 568       {
 569         printf ("0.0625 up: mantissa1 differs: %8x vs %8x\n",
 570                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 571         result = 1;
 572       }
 573     if (v1.ieee.exponent != v2.ieee.exponent)
 574       {
 575         printf ("0.0625 up: exponent differs: %4x vs %4x\n",
 576                 v1.ieee.exponent, v2.ieee.exponent);
 577         result = 1;
 578       }
 579     if (v1.ieee.negative != v2.ieee.negative)
 580       {
 581         printf ("0.0625 up: negative differs: %d vs %d\n",
 582                 v1.ieee.negative, v2.ieee.negative);
 583         result = 1;
 584       }
 585
 586     v1.d = d = -0.0625;
 587     d = nextafter (d, 0.0);
 588     v2.d = d = nextafter (d, -1.0);
 589
 590     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 591       {
 592         printf ("-0.0625 up: mantissa0 differs: %8x vs %8x\n",
 593                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 594         result = 1;
 595       }
 596     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 597       {
 598         printf ("-0.0625 up: mantissa1 differs: %8x vs %8x\n",
 599                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 600         result = 1;
 601       }
 602     if (v1.ieee.exponent != v2.ieee.exponent)
 603       {
 604         printf ("-0.0625 up: exponent differs: %4x vs %4x\n",
 605                 v1.ieee.exponent, v2.ieee.exponent);
 606         result = 1;
 607       }
 608     if (v1.ieee.negative != v2.ieee.negative)
 609       {
 610         printf ("-0.0625 up: negative differs: %d vs %d\n",
 611                 v1.ieee.negative, v2.ieee.negative);
 612         result = 1;
 613       }
 614
 615     v1.d = d = -0.0625;
 616     d = nextafter (d, -1.0);
 617     v2.d = d = nextafter (d, 0.0);
 618
 619     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 620       {
 621         printf ("-0.0625 down: mantissa0 differs: %8x vs %8x\n",
 622                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 623         result = 1;
 624       }
 625     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 626       {
 627         printf ("-0.0625 down: mantissa1 differs: %8x vs %8x\n",
 628                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 629         result = 1;
 630       }
 631     if (v1.ieee.exponent != v2.ieee.exponent)
 632       {
 633         printf ("-0.0625 down: exponent differs: %4x vs %4x\n",
 634                 v1.ieee.exponent, v2.ieee.exponent);
 635         result = 1;
 636       }
 637     if (v1.ieee.negative != v2.ieee.negative)
 638       {
 639         printf ("-0.0625 down: negative differs: %d vs %d\n",
 640                 v1.ieee.negative, v2.ieee.negative);
 641         result = 1;
 642       }
 643
 644     v1.d = d = 0.0;
 645     d = nextafter (d, 1.0);
 646     v2.d = nextafter (d, -1.0);
 647
 648     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 649       {
 650         printf ("0.0 up: mantissa0 differs: %8x vs %8x\n",
 651                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 652         result = 1;
 653       }
 654     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 655       {
 656         printf ("0.0 up: mantissa1 differs: %8x vs %8x\n",
 657                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 658         result = 1;
 659       }
 660     if (v1.ieee.exponent != v2.ieee.exponent)
 661       {
 662         printf ("0.0 up: exponent differs: %4x vs %4x\n",
 663                 v1.ieee.exponent, v2.ieee.exponent);
 664         result = 1;
 665       }
 666     if (0 != v2.ieee.negative)
 667       {
 668         printf ("0.0 up: negative differs: 0 vs %d\n",
 669                 v2.ieee.negative);
 670         result = 1;
 671       }
 672
 673     v1.d = d = 0.0;
 674     d = nextafter (d, -1.0);
 675     v2.d = nextafter (d, 1.0);
 676
 677     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 678       {
 679         printf ("0.0 down: mantissa0 differs: %8x vs %8x\n",
 680                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 681         result = 1;
 682       }
 683     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 684       {
 685         printf ("0.0 down: mantissa1 differs: %8x vs %8x\n",
 686                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 687         result = 1;
 688       }
 689     if (v1.ieee.exponent != v2.ieee.exponent)
 690       {
 691         printf ("0.0 down: exponent differs: %4x vs %4x\n",
 692                 v1.ieee.exponent, v2.ieee.exponent);
 693         result = 1;
 694       }
 695     if (1 != v2.ieee.negative)
 696       {
 697         printf ("0.0 down: negative differs: 1 vs %d\n",
 698                 v2.ieee.negative);
 699         result = 1;
 700       }
 701
 702     if (nextafter (0.0, INFINITY) != nextafter (0.0, 1.0)
 703         || nextafter (-0.0, INFINITY) != nextafter (-0.0, 1.0)
 704         || nextafter (0.0, -INFINITY) != nextafter (0.0, -1.0)
 705         || nextafter (-0.0, -INFINITY) != nextafter (-0.0, -1.0))
 706       {
 707         printf ("nextafter (+-0, +-Inf) != nextafter (+-0, +-1)\n");
 708         result = 1;
 709       }
 710
 711     if (nexttoward (0.0, INFINITY) != nexttoward (0.0, 1.0)
 712         || nexttoward (-0.0, INFINITY) != nexttoward (-0.0, 1.0)
 713         || nexttoward (0.0, -INFINITY) != nexttoward (0.0, -1.0)
 714         || nexttoward (-0.0, -INFINITY) != nexttoward (-0.0, -1.0))
 715       {
 716         printf ("nexttoward (+-0, +-Inf) != nexttoward (+-0, +-1)\n");
 717         result = 1;
 718       }
 719   }
 720
 721 #ifndef NO_LONG_DOUBLE
 722   {
 723     union ieee854_long_double v1;
 724     union ieee854_long_double v2;
 725     long double ld;
 726
 727     v1.d = ld = LDBL_MIN;
 728     if (fpclassify (ld) != FP_NORMAL)
 729       {
 730         printf ("fpclassify (LDBL_MIN) failed: %d\n", fpclassify (ld));
 731         result = 1;
 732       }
 733     ld = nextafterl (ld, LDBL_MIN / 2.0);
 734     if (fpclassify (ld) != FP_SUBNORMAL)
 735       {
 736         printf ("fpclassify (LDBL_MIN-epsilon) failed: %d (%La)\n",
 737                 fpclassify (ld), ld);
 738         result = 1;
 739       }
 740     v2.d = ld = nextafterl (ld, LDBL_MIN);
 741     if (fpclassify (ld) != FP_NORMAL)
 742       {
 743         printf ("fpclassify (LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
 744                 fpclassify (ld), ld);
 745         result = 1;
 746       }
 747
 748     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 749       {
 750         printf ("LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
 751                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 752         result = 1;
 753       }
 754     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 755       {
 756         printf ("LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
 757                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 758         result = 1;
 759       }
 760     if (v1.ieee.exponent != v2.ieee.exponent)
 761       {
 762         printf ("LDBL_MIN: exponent differs: %4x vs %4x\n",
 763                 v1.ieee.exponent, v2.ieee.exponent);
 764         result = 1;
 765       }
 766     if (v1.ieee.negative != v2.ieee.negative)
 767       {
 768         printf ("LDBL_MIN: negative differs: %d vs %d\n",
 769                 v1.ieee.negative, v2.ieee.negative);
 770         result = 1;
 771       }
 772
 773     v1.d = ld = -LDBL_MIN;
 774     if (fpclassify (ld) != FP_NORMAL)
 775       {
 776         printf ("fpclassify (-LDBL_MIN) failed: %d\n", fpclassify (ld));
 777         result = 1;
 778       }
 779     ld = nextafterl (ld, -LDBL_MIN / 2.0);
 780     if (fpclassify (ld) != FP_SUBNORMAL)
 781       {
 782         printf ("fpclassify (-LDBL_MIN-epsilon) failed: %d (%La)\n",
 783                 fpclassify (ld), ld);
 784         result = 1;
 785       }
 786     v2.d = ld = nextafterl (ld, -LDBL_MIN);
 787     if (fpclassify (ld) != FP_NORMAL)
 788       {
 789         printf ("fpclassify (-LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
 790                 fpclassify (ld), ld);
 791         result = 1;
 792       }
 793
 794     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 795       {
 796         printf ("-LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
 797                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 798         result = 1;
 799       }
 800     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 801       {
 802         printf ("-LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
 803                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 804         result = 1;
 805       }
 806     if (v1.ieee.exponent != v2.ieee.exponent)
 807       {
 808         printf ("-LDBL_MIN: exponent differs: %4x vs %4x\n",
 809                 v1.ieee.exponent, v2.ieee.exponent);
 810         result = 1;
 811       }
 812     if (v1.ieee.negative != v2.ieee.negative)
 813       {
 814         printf ("-LDBL_MIN: negative differs: %d vs %d\n",
 815                 v1.ieee.negative, v2.ieee.negative);
 816         result = 1;
 817       }
 818
 819     ld = LDBL_MAX;
 820     if (fpclassify (ld) != FP_NORMAL)
 821       {
 822         printf ("fpclassify (LDBL_MAX) failed: %d\n", fpclassify (ld));
 823         result = 1;
 824       }
 825     ld = nextafterl (ld, INFINITY);
 826     if (fpclassify (ld) != FP_INFINITE)
 827       {
 828         printf ("fpclassify (LDBL_MAX+epsilon) failed: %d\n", fpclassify (ld));
 829         result = 1;
 830       }
 831
 832     ld = -LDBL_MAX;
 833     if (fpclassify (ld) != FP_NORMAL)
 834       {
 835         printf ("fpclassify (-LDBL_MAX) failed: %d\n", fpclassify (ld));
 836         result = 1;
 837       }
 838     ld = nextafterl (ld, -INFINITY);
 839     if (fpclassify (ld) != FP_INFINITE)
 840       {
 841         printf ("fpclassify (-LDBL_MAX-epsilon) failed: %d\n",
 842                 fpclassify (ld));
 843         result = 1;
 844       }
 845
 846     v1.d = ld = 0.0625;
 847     ld = nextafterl (ld, 0.0);
 848     v2.d = ld = nextafterl (ld, 1.0);
 849
 850     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 851       {
 852         printf ("0.0625L down: mantissa0 differs: %8x vs %8x\n",
 853                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 854         result = 1;
 855       }
 856     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 857       {
 858         printf ("0.0625L down: mantissa1 differs: %8x vs %8x\n",
 859                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 860         result = 1;
 861       }
 862     if (v1.ieee.exponent != v2.ieee.exponent)
 863       {
 864         printf ("0.0625L down: exponent differs: %4x vs %4x\n",
 865                 v1.ieee.exponent, v2.ieee.exponent);
 866         result = 1;
 867       }
 868     if (v1.ieee.negative != v2.ieee.negative)
 869       {
 870         printf ("0.0625L down: negative differs: %d vs %d\n",
 871                 v1.ieee.negative, v2.ieee.negative);
 872         result = 1;
 873       }
 874
 875     v1.d = ld = 0.0625;
 876     ld = nextafterl (ld, 1.0);
 877     v2.d = ld = nextafterl (ld, 0.0);
 878
 879     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 880       {
 881         printf ("0.0625L up: mantissa0 differs: %8x vs %8x\n",
 882                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 883         result = 1;
 884       }
 885     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 886       {
 887         printf ("0.0625L up: mantissa1 differs: %8x vs %8x\n",
 888                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 889         result = 1;
 890       }
 891     if (v1.ieee.exponent != v2.ieee.exponent)
 892       {
 893         printf ("0.0625L up: exponent differs: %4x vs %4x\n",
 894                 v1.ieee.exponent, v2.ieee.exponent);
 895         result = 1;
 896       }
 897     if (v1.ieee.negative != v2.ieee.negative)
 898       {
 899         printf ("0.0625L up: negative differs: %d vs %d\n",
 900                 v1.ieee.negative, v2.ieee.negative);
 901         result = 1;
 902       }
 903
 904     v1.d = ld = -0.0625;
 905     ld = nextafterl (ld, 0.0);
 906     v2.d = ld = nextafterl (ld, -1.0);
 907
 908     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 909       {
 910         printf ("-0.0625L up: mantissa0 differs: %8x vs %8x\n",
 911                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 912         result = 1;
 913       }
 914     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 915       {
 916         printf ("-0.0625L up: mantissa1 differs: %8x vs %8x\n",
 917                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 918         result = 1;
 919       }
 920     if (v1.ieee.exponent != v2.ieee.exponent)
 921       {
 922         printf ("-0.0625L up: exponent differs: %4x vs %4x\n",
 923                 v1.ieee.exponent, v2.ieee.exponent);
 924         result = 1;
 925       }
 926     if (v1.ieee.negative != v2.ieee.negative)
 927       {
 928         printf ("-0.0625L up: negative differs: %d vs %d\n",
 929                 v1.ieee.negative, v2.ieee.negative);
 930         result = 1;
 931       }
 932
 933     v1.d = ld = -0.0625;
 934     ld = nextafterl (ld, -1.0);
 935     v2.d = ld = nextafterl (ld, 0.0);
 936
 937     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 938       {
 939         printf ("-0.0625L down: mantissa0 differs: %8x vs %8x\n",
 940                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 941         result = 1;
 942       }
 943     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 944       {
 945         printf ("-0.0625L down: mantissa1 differs: %8x vs %8x\n",
 946                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 947         result = 1;
 948       }
 949     if (v1.ieee.exponent != v2.ieee.exponent)
 950       {
 951         printf ("-0.0625L down: exponent differs: %4x vs %4x\n",
 952                 v1.ieee.exponent, v2.ieee.exponent);
 953         result = 1;
 954       }
 955     if (v1.ieee.negative != v2.ieee.negative)
 956       {
 957         printf ("-0.0625L down: negative differs: %d vs %d\n",
 958                 v1.ieee.negative, v2.ieee.negative);
 959         result = 1;
 960       }
 961
 962     v1.d = ld = 0.0;
 963     ld = nextafterl (ld, 1.0);
 964     v2.d = nextafterl (ld, -1.0);
 965
 966     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 967       {
 968         printf ("0.0L up: mantissa0 differs: %8x vs %8x\n",
 969                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 970         result = 1;
 971       }
 972     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 973       {
 974         printf ("0.0L up: mantissa1 differs: %8x vs %8x\n",
 975                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 976         result = 1;
 977       }
 978     if (v1.ieee.exponent != v2.ieee.exponent)
 979       {
 980         printf ("0.0L up: exponent differs: %4x vs %4x\n",
 981                 v1.ieee.exponent, v2.ieee.exponent);
 982         result = 1;
 983       }
 984     if (0 != v2.ieee.negative)
 985       {
 986         printf ("0.0L up: negative differs: 0 vs %d\n",
 987                 v2.ieee.negative);
 988         result = 1;
 989       }
 990
 991     v1.d = ld = 0.0;
 992     ld = nextafterl (ld, -1.0);
 993     v2.d = nextafterl (ld, 1.0);
 994
 995     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 996       {
 997         printf ("0.0L down: mantissa0 differs: %8x vs %8x\n",
 998                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 999         result = 1;
1000       }
1001     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
1002       {
1003         printf ("0.0L down: mantissa1 differs: %8x vs %8x\n",
1004                 v1.ieee.mantissa1, v2.ieee.mantissa1);
1005         result = 1;
1006       }
1007     if (v1.ieee.exponent != v2.ieee.exponent)
1008       {
1009         printf ("0.0L down: exponent differs: %4x vs %4x\n",
1010                 v1.ieee.exponent, v2.ieee.exponent);
1011         result = 1;
1012       }
1013     if (1 != v2.ieee.negative)
1014       {
1015         printf ("0.0L down: negative differs: 1 vs %d\n",
1016                 v2.ieee.negative);
1017         result = 1;
1018       }
1019
1020     if (nextafterl (0.0, INFINITY) != nextafterl (0.0, 1.0)
1021         || nextafterl (-0.0, INFINITY) != nextafterl (-0.0, 1.0)
1022         || nextafterl (0.0, -INFINITY) != nextafterl (0.0, -1.0)
1023         || nextafterl (-0.0, -INFINITY) != nextafterl (-0.0, -1.0))
1024       {
1025         printf ("nextafterl (+-0, +-Inf) != nextafterl (+-0, +-1)\n");
1026         result = 1;
1027       }
1028
1029     if (nexttowardl (0.0L, INFINITY) != nexttowardl (0.0L, 1.0L)
1030         || nexttowardl (-0.0L, INFINITY) != nexttowardl (-0.0L, 1.0L)
1031         || nexttowardl (0.0L, -INFINITY) != nexttowardl (0.0L, -1.0L)
1032         || nexttowardl (-0.0L, -INFINITY) != nexttowardl (-0.0L, -1.0L))
1033       {
1034         printf ("nexttowardl (+-0, +-Inf) != nexttowardl (+-0, +-1)\n");
1035         result = 1;
1036       }
1037   }
1038 #endif
1039
1040   if (! isnormal (FLT_MIN))
1041     {
1042       puts ("isnormal (FLT_MIN) failed");
1043       result = 1;
1044     }
1045   if (! isnormal (DBL_MIN))
1046     {
1047       puts ("isnormal (DBL_MIN) failed");
1048       result = 1;
1049     }
1050 #ifndef NO_LONG_DOUBLE
1051   if (! isnormal (LDBL_MIN))
1052     {
1053       puts ("isnormal (LDBL_MIN) failed");
1054       result = 1;
1055     }
1056 #endif
1057
1058 #ifdef __i386__
1059   /* This is a test for the strange long doubles in x86 FPUs.  */
1060   {
1061     union
1062     {
1063       char b[10];
1064       long double d;
1065     } u =
1066       { .b = { 0, 0, 0, 0, 0, 0, 0, 0x80, 0, 0 } };
1067
1068     if (fpclassify (u.d) != FP_NORMAL)
1069       {
1070         printf ("fpclassify (0x00008000000000000000) failed: %d (%Lg)\n",
1071                 fpclassify (u.d), u.d);
1072         result = 1;
1073       }
1074   }
1075
1076   /* Special NaNs in x86 long double.  Test for scalbl.  */
1077   {
1078     union
1079     {
1080       char b[10];
1081       long double d;
1082     } u =
1083       { .b = { 0, 1, 0, 0, 0, 0, 0, 0xc0, 0xff, 0x7f } };
1084     long double r;
1085
1086     r = scalbl (u.d, 0.0);
1087     if (!isnan (r))
1088       {
1089         puts ("scalbl(NaN, 0) does not return NaN");
1090         result = 1;
1091       }
1092     else if (memcmp (&r, &u.d, sizeof (double)) != 0)
1093       {
1094         puts ("scalbl(NaN, 0) does not return the same NaN");
1095         result = 1;
1096       }
1097   }
1098 #endif
1099
1100 #ifndef NO_LONG_DOUBLE
1101   {
1102     long double r;
1103
1104     feclearexcept (FE_ALL_EXCEPT);
1105     r = scalbl (LDBL_MIN, 2147483647);
1106     if (! isinf (r))
1107       {
1108         puts ("scalbl (LDBL_MIN, 2147483647) does not return Inf");
1109         result = 1;
1110       }
1111     else if (signbit (r) != 0)
1112       {
1113         puts ("scalbl (LDBL_MIN, 2147483647) returns -Inf");
1114         result = 1;
1115       }
1116     else if (fetestexcept (FE_UNDERFLOW))
1117       {
1118         puts ("scalbl(NaN, 0) raises underflow exception");
1119         result = 1;
1120       }
1121
1122     feclearexcept (FE_ALL_EXCEPT);
1123     r = scalbl (LDBL_MAX, -2147483647);
1124     if (r != 0.0)
1125       {
1126         puts ("scalbl (LDBL_MAX, -2147483647) does not return 0");
1127         result = 1;
1128       }
1129     else if (signbit (r) != 0)
1130       {
1131         puts ("scalbl (LDBL_MAX, -2147483647) returns -Inf");
1132         result = 1;
1133       }
1134     else if (fetestexcept (FE_OVERFLOW))
1135       {
1136         puts ("scalbl(NaN, 0) raises overflow exception");
1137         result = 1;
1138       }
1139   }
1140 #endif
1141
1142   return result;
1143 }