tests/test-ldexpl.c

   1 /* Test of multiplying a 'long double' by a power of 2.
   2    Copyright (C) 2007-2010 Free Software Foundation, Inc.
   3
   4    This program is free software: you can redistribute it and/or modify
   5    it under the terms of the GNU General Public License as published by
   6    the Free Software Foundation; either version 3 of the License, or
   7    (at your option) any later version.
   8
   9    This program is distributed in the hope that it will be useful,
  10    but WITHOUT ANY WARRANTY; without even the implied warranty of
  11    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12    GNU General Public License for more details.
  13
  14    You should have received a copy of the GNU General Public License
  15    along with this program.  If not, see <http://www.gnu.org/licenses/>.  */
  16
  17 /* Written by Bruno Haible <bruno@clisp.org>, 2007.  */
  18
  19 #include <config.h>
  20
  21 #include <math.h>
  22
  23 #include "signature.h"
  24 SIGNATURE_CHECK (ldexpl, long double, (long double, int));
  25
  26 #include <float.h>
  27
  28 #include "fpucw.h"
  29 #include "isnanl-nolibm.h"
  30 #include "nan.h"
  31 #include "macros.h"
  32
  33 /* On HP-UX 10.20, negating 0.0L does not yield -0.0L.
  34    So we use minus_zero instead.
  35    IRIX cc can't put -0.0L into .data, but can compute at runtime.
  36    Note that the expression -LDBL_MIN * LDBL_MIN does not work on other
  37    platforms, such as when cross-compiling to PowerPC on MacOS X 10.5.  */
  38 #if defined __hpux || defined __sgi
  39 static long double
  40 compute_minus_zero (void)
  41 {
  42   return -LDBL_MIN * LDBL_MIN;
  43 }
  44 # define minus_zero compute_minus_zero ()
  45 #else
  46 long double minus_zero = -0.0L;
  47 #endif
  48
  49 int
  50 main ()
  51 {
  52   int i;
  53   long double x;
  54   long double y;
  55   DECL_LONG_DOUBLE_ROUNDING
  56
  57   BEGIN_LONG_DOUBLE_ROUNDING ();
  58
  59   { /* NaN.  */
  60     x = NaNl ();
  61     y = ldexpl (x, 0); ASSERT (isnanl (y));
  62     y = ldexpl (x, 5); ASSERT (isnanl (y));
  63     y = ldexpl (x, -5); ASSERT (isnanl (y));
  64   }
  65
  66   { /* Positive infinity.  */
  67     x = 1.0L / 0.0L;
  68     y = ldexpl (x, 0); ASSERT (y == x);
  69     y = ldexpl (x, 5); ASSERT (y == x);
  70     y = ldexpl (x, -5); ASSERT (y == x);
  71   }
  72
  73   { /* Negative infinity.  */
  74     x = -1.0L / 0.0L;
  75     y = ldexpl (x, 0); ASSERT (y == x);
  76     y = ldexpl (x, 5); ASSERT (y == x);
  77     y = ldexpl (x, -5); ASSERT (y == x);
  78   }
  79
  80   { /* Positive zero.  */
  81     x = 0.0L;
  82     y = ldexpl (x, 0); ASSERT (y == x); ASSERT (!signbit (x));
  83     y = ldexpl (x, 5); ASSERT (y == x); ASSERT (!signbit (x));
  84     y = ldexpl (x, -5); ASSERT (y == x); ASSERT (!signbit (x));
  85   }
  86
  87   { /* Negative zero.  */
  88     x = minus_zero;
  89     y = ldexpl (x, 0); ASSERT (y == x); ASSERT (signbit (x));
  90     y = ldexpl (x, 5); ASSERT (y == x); ASSERT (signbit (x));
  91     y = ldexpl (x, -5); ASSERT (y == x); ASSERT (signbit (x));
  92   }
  93
  94   { /* Positive finite number.  */
  95     x = 1.73205L;
  96     y = ldexpl (x, 0); ASSERT (y == x);
  97     y = ldexpl (x, 5); ASSERT (y == x * 32.0L);
  98     y = ldexpl (x, -5); ASSERT (y == x * 0.03125L);
  99   }
 100
 101   { /* Negative finite number.  */
 102     x = -20.085536923187667742L;
 103     y = ldexpl (x, 0); ASSERT (y == x);
 104     y = ldexpl (x, 5); ASSERT (y == x * 32.0L);
 105     y = ldexpl (x, -5); ASSERT (y == x * 0.03125L);
 106   }
 107
 108   for (i = 1, x = 1.73205L; i <= LDBL_MAX_EXP; i++, x *= 2.0L)
 109     {
 110       y = ldexpl (x, 0); ASSERT (y == x);
 111       y = ldexpl (x, 5); ASSERT (y == x * 32.0L);
 112       y = ldexpl (x, -5); ASSERT (y == x * 0.03125L);
 113     }
 114   for (i = 1, x = 1.73205L; i >= LDBL_MIN_EXP; i--, x *= 0.5L)
 115     {
 116       y = ldexpl (x, 0); ASSERT (y == x);
 117       y = ldexpl (x, 5); ASSERT (y == x * 32.0L);
 118       if (i - 5 >= LDBL_MIN_EXP)
 119         {
 120           y = ldexpl (x, -5); ASSERT (y == x * 0.03125L);
 121         }
 122     }
 123   for (; i >= LDBL_MIN_EXP - 100 && x > 0.0L; i--, x *= 0.5L)
 124     {
 125       y = ldexpl (x, 0); ASSERT (y == x);
 126       y = ldexpl (x, 5); ASSERT (y == x * 32.0L);
 127     }
 128
 129   return 0;
 130 }