math/bug-nexttoward.c

   1 #include <fenv.h>
   2 #include <math.h>
   3 #include <float.h>
   4 #include <stdlib.h>
   5 #include <stdio.h>
   6 #include <math-tests.h>
   7
   8 #if !defined(FE_OVERFLOW) && !defined(FE_UNDERFLOW)
   9 /* If there's no support for the exceptions this test is checking,
  10    then just return success and allow the test to be compiled.  */
  11 # define fetestexcept(e) 1
  12 #endif
  13
  14 float zero = 0.0;
  15 float inf = INFINITY;
  16
  17 int
  18 main (void)
  19 {
  20   int result = 0;
  21
  22   long double tl = (long double) FLT_MAX + 0x1.0p128L;
  23   float fi = INFINITY;
  24   float m = FLT_MAX;
  25   feclearexcept (FE_ALL_EXCEPT);
  26   if (nexttowardf (m, tl) != fi)
  27     {
  28       puts ("nexttowardf+ failed");
  29       ++result;
  30     }
  31   if (EXCEPTION_TESTS (float) && fetestexcept (FE_OVERFLOW) == 0)
  32     {
  33       puts ("nexttowardf+ did not overflow");
  34       ++result;
  35     }
  36   feclearexcept (FE_ALL_EXCEPT);
  37   if (nexttowardf (-m, -tl) != -fi)
  38     {
  39       puts ("nexttowardf- failed");
  40       ++result;
  41     }
  42   if (EXCEPTION_TESTS (float) && fetestexcept (FE_OVERFLOW) == 0)
  43     {
  44       puts ("nexttowardf- did not overflow");
  45       ++result;
  46     }
  47
  48   fi = 0;
  49   m = FLT_MIN;
  50   feclearexcept (FE_ALL_EXCEPT);
  51   fi = nexttowardf (m, fi);
  52   if (fi < 0 || fi >= FLT_MIN)
  53     {
  54       puts ("nexttowardf+ failed");
  55       ++result;
  56     }
  57   if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
  58     {
  59       puts ("nexttowardf+ did not underflow");
  60       ++result;
  61     }
  62   fi = 0;
  63   feclearexcept (FE_ALL_EXCEPT);
  64   fi = nexttowardf (-m, -fi);
  65   if (fi > 0 || fi <= -FLT_MIN)
  66     {
  67       puts ("nexttowardf- failed");
  68       ++result;
  69     }
  70   if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
  71     {
  72       puts ("nexttowardf- did not underflow");
  73       ++result;
  74     }
  75   fi = -INFINITY;
  76   feclearexcept (FE_ALL_EXCEPT);
  77   m = nexttowardf (zero, inf);
  78   if (m < 0.0 || m >= FLT_MIN)
  79     {
  80       puts ("nexttowardf+ failed");
  81       ++result;
  82     }
  83   if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
  84     {
  85       puts ("nexttowardf+ did not underflow");
  86       ++result;
  87     }
  88   feclearexcept (FE_ALL_EXCEPT);
  89   if (nexttowardf (m, fi) != 0.0)
  90     {
  91       puts ("nexttowardf+ failed");
  92       ++result;
  93     }
  94   if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
  95     {
  96       puts ("nexttowardf+ did not underflow");
  97       ++result;
  98     }
  99   feclearexcept (FE_ALL_EXCEPT);
 100   m = nexttowardf (copysignf (zero, -1.0), -inf);
 101   if (m > 0.0 || m <= -FLT_MIN)
 102     {
 103       puts ("nexttowardf- failed");
 104       ++result;
 105     }
 106   if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
 107     {
 108       puts ("nexttowardf- did not underflow");
 109       ++result;
 110     }
 111   feclearexcept (FE_ALL_EXCEPT);
 112   if (nexttowardf (m, -fi) != 0.0)
 113     {
 114       puts ("nexttowardf- failed");
 115       ++result;
 116     }
 117   if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
 118     {
 119       puts ("nexttowardf- did not underflow");
 120       ++result;
 121     }
 122
 123   tl = (long double) DBL_MAX + 1.0e305L;
 124   double di = INFINITY;
 125   double dm = DBL_MAX;
 126   feclearexcept (FE_ALL_EXCEPT);
 127   if (nexttoward (dm, tl) != di)
 128     {
 129       puts ("nexttoward+ failed");
 130       ++result;
 131     }
 132   if (EXCEPTION_TESTS (double) && fetestexcept (FE_OVERFLOW) == 0)
 133     {
 134       puts ("nexttoward+ did not overflow");
 135       ++result;
 136     }
 137   feclearexcept (FE_ALL_EXCEPT);
 138   if (nexttoward (-dm, -tl) != -di)
 139     {
 140       puts ("nexttoward- failed");
 141       ++result;
 142     }
 143   if (EXCEPTION_TESTS (double) && fetestexcept (FE_OVERFLOW) == 0)
 144     {
 145       puts ("nexttoward- did not overflow");
 146       ++result;
 147     }
 148
 149   di = 0;
 150   dm = DBL_MIN;
 151   feclearexcept (FE_ALL_EXCEPT);
 152   di = nexttoward (dm, di);
 153   if (di < 0 || di >= DBL_MIN)
 154     {
 155       puts ("nexttoward+ failed");
 156       ++result;
 157     }
 158   if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
 159     {
 160       puts ("nexttoward+ did not underflow");
 161       ++result;
 162     }
 163   di = 0;
 164   feclearexcept (FE_ALL_EXCEPT);
 165   di = nexttoward (-dm, -di);
 166   if (di > 0 || di <= -DBL_MIN)
 167     {
 168       puts ("nexttoward- failed");
 169       ++result;
 170     }
 171   if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
 172     {
 173       puts ("nexttoward- did not underflow");
 174       ++result;
 175     }
 176   di = -INFINITY;
 177   feclearexcept (FE_ALL_EXCEPT);
 178   dm = nexttoward (zero, inf);
 179   if (dm < 0.0 || dm >= DBL_MIN)
 180     {
 181       puts ("nexttoward+ failed");
 182       ++result;
 183     }
 184   if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
 185     {
 186       puts ("nexttoward+ did not underflow");
 187       ++result;
 188     }
 189   feclearexcept (FE_ALL_EXCEPT);
 190   if (nexttoward (dm, di) != 0.0)
 191     {
 192       puts ("nexttoward+ failed");
 193       ++result;
 194     }
 195   if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
 196     {
 197       puts ("nexttoward+ did not underflow");
 198       ++result;
 199     }
 200   feclearexcept (FE_ALL_EXCEPT);
 201   dm = nexttoward (copysign (zero, -1.0), -inf);
 202   if (dm > 0.0 || dm <= -DBL_MIN)
 203     {
 204       puts ("nexttoward- failed");
 205       ++result;
 206     }
 207   if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
 208     {
 209       puts ("nexttoward- did not underflow");
 210       ++result;
 211     }
 212   feclearexcept (FE_ALL_EXCEPT);
 213   if (nexttoward (dm, -di) != 0.0)
 214     {
 215       puts ("nexttoward- failed");
 216       ++result;
 217     }
 218   if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
 219     {
 220       puts ("nexttoward- did not underflow");
 221       ++result;
 222     }
 223
 224   long double li = INFINITY;
 225   long double lm = LDBL_MAX;
 226   feclearexcept (FE_ALL_EXCEPT);
 227   if (nexttowardl (lm, li) != li)
 228     {
 229       puts ("nexttowardl+ failed");
 230       ++result;
 231     }
 232   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_OVERFLOW) == 0)
 233     {
 234       puts ("nexttowardl+ did not overflow");
 235       ++result;
 236     }
 237   feclearexcept (FE_ALL_EXCEPT);
 238   if (nexttowardl (-lm, -li) != -li)
 239     {
 240       puts ("nexttowardl failed");
 241       ++result;
 242     }
 243   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_OVERFLOW) == 0)
 244     {
 245       puts ("nexttowardl- did not overflow");
 246       ++result;
 247     }
 248
 249   li = 0;
 250   lm = LDBL_MIN;
 251   feclearexcept (FE_ALL_EXCEPT);
 252   li = nexttowardl (lm, li);
 253   if (li < 0 || li >= LDBL_MIN)
 254     {
 255       puts ("nexttowardl+ failed");
 256       ++result;
 257     }
 258   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
 259     {
 260       puts ("nexttowardl+ did not underflow");
 261       ++result;
 262     }
 263   li = 0;
 264   feclearexcept (FE_ALL_EXCEPT);
 265   li = nexttowardl (-lm, -li);
 266   if (li > 0 || li <= -LDBL_MIN)
 267     {
 268       puts ("nexttowardl- failed");
 269       ++result;
 270     }
 271   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
 272     {
 273       puts ("nexttowardl- did not underflow");
 274       ++result;
 275     }
 276   li = -INFINITY;
 277   feclearexcept (FE_ALL_EXCEPT);
 278   lm = nexttowardl (zero, inf);
 279   if (lm < 0.0 || lm >= LDBL_MIN)
 280     {
 281       puts ("nexttowardl+ failed");
 282       ++result;
 283     }
 284   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
 285     {
 286       puts ("nexttowardl+ did not underflow");
 287       ++result;
 288     }
 289   feclearexcept (FE_ALL_EXCEPT);
 290   if (nexttowardl (lm, li) != 0.0)
 291     {
 292       puts ("nexttowardl+ failed");
 293       ++result;
 294     }
 295   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
 296     {
 297       puts ("nexttowardl+ did not underflow");
 298       ++result;
 299     }
 300   feclearexcept (FE_ALL_EXCEPT);
 301   lm = nexttowardl (copysign (zero, -1.0), -inf);
 302   if (lm > 0.0 || lm <= -LDBL_MIN)
 303     {
 304       puts ("nexttowardl- failed");
 305       ++result;
 306     }
 307   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
 308     {
 309       puts ("nexttowardl- did not underflow");
 310       ++result;
 311     }
 312   feclearexcept (FE_ALL_EXCEPT);
 313   if (nexttowardl (lm, -li) != 0.0)
 314     {
 315       puts ("nexttowardl- failed");
 316       ++result;
 317     }
 318   if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
 319     {
 320       puts ("nexttowardl- did not underflow");
 321       ++result;
 322     }
 323
 324   return result;
 325 }