math/test-misc.c

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