libstdc++-v3/libmath/mathconf.h

   1 /* Configuration data for libmath subpart of libstdc++. */
   2
   3 /* Copyright (C) 1997-1999, 2000, 2001 Free Software Foundation, Inc.
   4
   5    This file is part of the GNU ISO C++ Library.  This library is free
   6    software; you can redistribute it and/or modify it under the
   7    terms of the GNU General Public License as published by the
   8    Free Software Foundation; either version 2, or (at your option)
   9    any later version.
  10
  11    This library is distributed in the hope that it will be useful,
  12    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14    GNU General Public License for more details.
  15
  16    You should have received a copy of the GNU General Public License along
  17    with this library; see the file COPYING.  If not, write to the Free
  18    Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301,
  19    USA.
  20
  21    As a special exception, you may use this file as part of a free software
  22    library without restriction.  Specifically, if other files instantiate
  23    templates or use macros or inline functions from this file, or you compile
  24    this file and link it with other files to produce an executable, this
  25    file does not by itself cause the resulting executable to be covered by
  26    the GNU General Public License.  This exception does not however
  27    invalidate any other reasons why the executable file might be covered by
  28    the GNU General Public License.  */
  29
  30
  31 #include <config.h>
  32
  33 #ifdef HAVE_ENDIAN_H
  34 # include <endian.h>
  35 #else
  36 # ifdef HAVE_MACHINE_ENDIAN_H
  37 #  ifdef HAVE_SYS_TYPES_H
  38 #   include <sys/types.h>
  39 #  endif
  40 #  include <machine/endian.h>
  41 # else
  42 #  ifdef HAVE_SYS_MACHINE_H
  43 #   include <sys/machine.h>
  44 #  else
  45 #   if defined HAVE_SYS_ISA_DEFS_H || defined HAVE_MACHINE_PARAM_H
  46 /* This is on Solaris.  */
  47 #    ifdef HAVE_SYS_ISA_DEFS_H
  48 #     include <sys/isa_defs.h>
  49 #    endif
  50 #    ifdef HAVE_MACHINE_PARAM_H
  51 #     include <machine/param.h>
  52 #    endif
  53 #    ifdef _LITTLE_ENDIAN
  54 #     define LITTLE_ENDIAN 1
  55 #    endif
  56 #    ifdef _BIG_ENDIAN
  57 #     define BIG_ENDIAN 1
  58 #    endif
  59 #    define BYTE_ORDER 1
  60 #   else
  61 /* We have to rely on the AC_C_BIGENDIAN test.  */
  62 #    ifdef WORDS_BIGENDIAN
  63 #     define BIG_ENDIAN 1
  64 #    else
  65 #     define LITTLE_ENDIAN 1
  66 #    endif
  67 #    define BYTE_ORDER 1
  68 #   endif
  69 #  endif
  70 # endif
  71 #endif
  72
  73 typedef unsigned int U_int32_t __attribute ((mode (SI)));
  74 typedef int Int32_t __attribute ((mode (SI)));
  75 typedef unsigned int U_int64_t __attribute ((mode (DI)));
  76 typedef int Int64_t __attribute ((mode (DI)));
  77
  78 #ifdef HAVE_NAN_H
  79 # include <nan.h>
  80 #endif
  81
  82 #ifdef HAVE_IEEEFP_H
  83 # include <ieeefp.h>
  84 #endif
  85
  86 #ifdef HAVE_FP_H
  87 # include <fp.h>
  88 #endif
  89
  90 #ifdef HAVE_FLOAT_H
  91 # include <float.h>
  92 #endif
  93
  94 /* `float' variant of HUGE_VAL.  */
  95 #ifndef HUGE_VALF
  96 # ifdef HUGE_VALf
  97 #  define HUGE_VALF HUGE_VALf
  98 # else
  99 #  define HUGE_VALF HUGE_VAL
 100 # endif
 101 #endif
 102
 103 /* `long double' variant of HUGE_VAL.  */
 104 #ifndef HUGE_VALL
 105 # ifdef HUGE_VALl
 106 #  define HUGE_VALL HUGE_VALl
 107 # else
 108 #  define HUGE_VALL HUGE_VAL
 109 # endif
 110 #endif
 111
 112 /* Make sure that at least HUGE_VAL is defined.  */
 113 #ifndef HUGE_VAL
 114 # ifdef HUGE
 115 #  define HUGE_VAL HUGE
 116 # else
 117 #  ifdef MAXFLOAT
 118 #   define HUGE_VAL MAXFLOAT
 119 #  else
 120 #   error "We need HUGE_VAL!"
 121 #  endif
 122 # endif
 123 #endif
 124
 125 #ifndef M_PI
 126 # define M_PI 3.14159265358979323846
 127 #endif
 128
 129
 130 #ifdef __cplusplus
 131 extern "C" {
 132 #endif
 133
 134 /* signbit is a macro in ISO C99.  */
 135 #ifndef signbit
 136 extern int __signbitf (float);
 137 extern int __signbit (double);
 138 extern int __signbitl (long double);
 139
 140 # define signbit(x) \
 141      (sizeof (x) == sizeof (float) ?                                          \
 142         __signbitf (x)                                                        \
 143       : sizeof (x) == sizeof (double) ?                                       \
 144         __signbit (x) : __signbitl (x))
 145 #endif
 146
 147 #if BYTE_ORDER == BIG_ENDIAN
 148 typedef union
 149 {
 150   double value;
 151   struct
 152   {
 153     U_int32_t msw;
 154     U_int32_t lsw;
 155   } parts;
 156 } ieee_double_shape_type;
 157 #endif
 158 #if BYTE_ORDER == LITTLE_ENDIAN
 159 typedef union
 160 {
 161   double value;
 162   struct
 163   {
 164     U_int32_t lsw;
 165     U_int32_t msw;
 166   } parts;
 167 } ieee_double_shape_type;
 168 #endif
 169 /* Get the more significant 32 bit int from a double.  */
 170 #define GET_HIGH_WORD(i,d)                                      \
 171 do {                                                            \
 172   ieee_double_shape_type gh_u;                                  \
 173   gh_u.value = (d);                                             \
 174   (i) = gh_u.parts.msw;                                         \
 175 } while (0)
 176
 177
 178 typedef union
 179 {
 180   float value;
 181   U_int32_t word;
 182 } ieee_float_shape_type;
 183 /* Get a 32 bit int from a float.  */
 184 #define GET_FLOAT_WORD(i,d)                                     \
 185 do {                                                            \
 186   ieee_float_shape_type gf_u;                                   \
 187   gf_u.value = (d);                                             \
 188   (i) = gf_u.word;                                              \
 189 } while (0)
 190
 191
 192 #if BYTE_ORDER == BIG_ENDIAN
 193 typedef union
 194 {
 195   long double value;
 196   struct
 197   {
 198     unsigned int sign_exponent:16;
 199     unsigned int empty:16;
 200     U_int32_t msw;
 201     U_int32_t lsw;
 202   } parts;
 203 } ieee_long_double_shape_type;
 204 #endif
 205 #if BYTE_ORDER == LITTLE_ENDIAN
 206 typedef union
 207 {
 208   long double value;
 209   struct
 210   {
 211     U_int32_t lsw;
 212     U_int32_t msw;
 213     unsigned int sign_exponent:16;
 214     unsigned int empty:16;
 215   } parts;
 216 } ieee_long_double_shape_type;
 217 #endif
 218 /* Get int from the exponent of a long double.  */
 219 #define GET_LDOUBLE_EXP(exp,d)                                  \
 220 do {                                                            \
 221   ieee_long_double_shape_type ge_u;                             \
 222   ge_u.value = (d);                                             \
 223   (exp) = ge_u.parts.sign_exponent;                             \
 224 } while (0)
 225
 226 #if BYTE_ORDER == BIG_ENDIAN
 227 typedef union
 228 {
 229   long double value;
 230   struct
 231   {
 232     U_int64_t msw;
 233     U_int64_t lsw;
 234   } parts64;
 235   struct
 236   {
 237     U_int32_t w0, w1, w2, w3;
 238   } parts32;
 239 } ieee_quad_double_shape_type;
 240 #endif
 241 #if BYTE_ORDER == LITTLE_ENDIAN
 242 typedef union
 243 {
 244   long double value;
 245   struct
 246   {
 247     U_int64_t lsw;
 248     U_int64_t msw;
 249   } parts64;
 250   struct
 251   {
 252     U_int32_t w3, w2, w1, w0;
 253   } parts32;
 254 } ieee_quad_double_shape_type;
 255 #endif
 256 /* Get most significant 64 bit int from a quad long double.  */
 257 #define GET_LDOUBLE_MSW64(msw,d)                                \
 258 do {                                                            \
 259   ieee_quad_double_shape_type qw_u;                             \
 260   qw_u.value = (d);                                             \
 261   (msw) = qw_u.parts64.msw;                                     \
 262 } while (0)
 263
 264
 265 /* Replacement for non-existing float functions.  */
 266 #if !defined(HAVE_FABSF) && !defined(HAVE___BUILTIN_FABSF)
 267 # define fabsf(x) fabs (x)
 268 #endif
 269 #if !defined(HAVE_COSF) && !defined(HAVE___BUILTIN_COSF)
 270 # define cosf(x) cos (x)
 271 #endif
 272 #ifndef HAVE_COSHF
 273 # define coshf(x) cosh (x)
 274 #endif
 275 #ifndef HAVE_EXPF
 276 # define expf(x) expf (x)
 277 #endif
 278 #ifndef HAVE_LOGF
 279 # define logf(x) log(x)
 280 #endif
 281 #ifndef HAVE_LOG10F
 282 # define log10f(x) log10 (x)
 283 #endif
 284 #ifndef HAVE_POWF
 285 # define powf(x, y) pow (x, y)
 286 #endif
 287 #if !defined(HAVE_SINF) && !defined(HAVE___BUILTIN_SINF)
 288 # define sinf(x) sin (x)
 289 #endif
 290 #ifndef HAVE_SINHF
 291 # define sinhf(x) sinh (x)
 292 #endif
 293 #if !defined(HAVE_SQRTF) && !defined(HAVE___BUILTIN_SQRTF)
 294 # define sqrtf(x) sqrt (x)
 295 #endif
 296 #ifndef HAVE_TANF
 297 # define tanf(x) tan (x)
 298 #endif
 299 #ifndef HAVE_TANHF
 300 # define tanhf(x) tanh (x)
 301 #endif
 302 #ifndef HAVE_STRTOF
 303 # define strtof(s, e) strtod (s, e)
 304 #endif
 305
 306 #ifdef __cplusplus
 307 }
 308 #endif
 309