libstdc++-v3/math/mathconf.h

   1 /* Configuration data for libmath subpart of libstdc++. */
   2
   3 /* Copyright (C) 1997-1999, 2000 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, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
  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 <bits/c++config.h>
  32
  33 #ifdef _GLIBCPP_HAVE_ENDIAN_H
  34 # include <endian.h>
  35 #else
  36 # ifdef _GLIBCPP_HAVE_MACHINE_ENDIAN_H
  37 #  include <machine/endian.h>
  38 # else
  39 #  ifdef _GLIBCPP_HAVE_SYS_MACHINE_H
  40 #   include <sys/machine.h>
  41 #  else
  42 #   if defined _GLIBCPP_HAVE_SYS_ISA_DEFS_H || defined _GLIBCPP_HAVE_MACHINE_PARAM_H
  43 /* This is on Solaris.  */
  44 #    ifdef _GLIBCPP_HAVE_SYS_ISA_DEFS_H
  45 #     include <sys/isa_defs.h>
  46 #    endif
  47 #    ifdef _GLIBCPP_HAVE_MACHINE_PARAM_H
  48 #     include <machine/param.h>
  49 #    endif
  50 #    ifdef _LITTLE_ENDIAN
  51 #     define LITTLE_ENDIAN 1
  52 #    endif
  53 #    ifdef _BIG_ENDIAN
  54 #     define BIG_ENDIAN 1
  55 #    endif
  56 #    define BYTE_ORDER 1
  57 #   else
  58 /* We have to rely on the AC_C_BIGENDIAN test.  */
  59 #    ifdef WORDS_BIGENDIAN
  60 #     define BIG_ENDIAN 1
  61 #    else
  62 #     define LITTLE_ENDIAN 1
  63 #    endif
  64 #    define BYTE_ORDER 1
  65 #   endif
  66 #  endif
  67 # endif
  68 #endif
  69
  70 typedef unsigned int U_int32_t __attribute ((mode (SI)));
  71 typedef int Int32_t __attribute ((mode (SI)));
  72 typedef unsigned int U_int64_t __attribute ((mode (DI)));
  73 typedef int Uint64_t __attribute ((mode (DI)));
  74
  75 #ifdef _GLIBCPP_HAVE_NAN_H
  76 # include <nan.h>
  77 #endif
  78
  79 #ifndef NAN
  80 # define NAN (nan())
  81 double nan (void);
  82 #endif
  83
  84 #ifdef _GLIBCPP_HAVE_IEEEFP_H
  85 # include <ieeefp.h>
  86 #endif
  87
  88 #ifdef _GLIBCPP_HAVE_FP_H
  89 # include <fp.h>
  90 #endif
  91
  92 #ifdef _GLIBCPP_HAVE_FLOAT_H
  93 # include <float.h>
  94 #endif
  95
  96 /* `float' variant of HUGE_VAL.  */
  97 #ifndef HUGE_VALF
  98 # ifdef HUGE_VALf
  99 #  define HUGE_VALF HUGE_VALf
 100 # else
 101 #  define HUGE_VALF HUGE_VAL
 102 # endif
 103 #endif
 104
 105 /* `long double' variant of HUGE_VAL.  */
 106 #ifndef HUGE_VALL
 107 # ifdef HUGE_VALl
 108 #  define HUGE_VALL HUGE_VALl
 109 # else
 110 #  define HUGE_VALL HUGE_VAL
 111 # endif
 112 #endif
 113
 114 /* Make sure that at least HUGE_VAL is defined.  */
 115 #ifndef HUGE_VAL
 116 # ifdef HUGE
 117 #  define HUGE_VAL HUGE
 118 # else
 119 #  ifdef MAXFLOAT
 120 #   define HUGE_VAL MAXFLOAT
 121 #  else
 122 #   error "We need HUGE_VAL!"
 123 #  endif
 124 # endif
 125 #endif
 126
 127 #ifndef M_PI
 128 # define M_PI 3.14159265358979323846
 129 #endif
 130
 131 /* Test whether number is finite.  */
 132 #ifdef isfinite
 133 /* This is an ISO C 9x function.  */
 134 # define FINITE_P(X) isfinite (X)
 135 # define FINITEF_P(X) isfinite (X)
 136 # define FINITEL_P(X) isfinite (X)
 137 #else
 138 # ifdef IsNANorINF
 139 /* This is for Solaris, which does not have special macros for other
 140    types.  */
 141 #  define FINITE_P(X) (!IsNANorINF (X))
 142 #  define FINITEF_P(X) (!IsNANorINF (X))
 143 #  define FINITEL_P(X) (!IsNANorINF (X))
 144 # else
 145 #  if defined _GLIBCPP_HAVE_ISINF && defined _GLIBCPP_HAVE_ISNAN
 146 #   define FINITE_P(X) ({ double __x = (X); !isinf (__x) && !isnan (__x); })
 147 #  else
 148 #   ifdef _GLIBCPP_HAVE_FINITE
 149 #    define FINITE_P(X) finite (X)
 150 #   else
 151 #    error "We need FINITE_P"
 152 #   endif
 153 #  endif
 154 #  if defined _GLIBCPP_HAVE_ISINFF && defined _GLIBCPP_HAVE_ISNANF
 155 #   define FINITEF_P(X) ({ float __x = (X); !isinff (__x) && !isnanf (__x); })
 156 #  else
 157 #   ifdef _GLIBCPP_HAVE_FINITE
 158 #    define FINITEF_P(X) finite (X)
 159 #   else
 160 #    define FINITEF_P(X) FINITE_P (X)
 161 #   endif
 162 #  endif
 163 #  if defined _GLIBCPP_HAVE_ISINFL && defined _GLIBCPP_HAVE_ISNANL
 164 #   define FINITEL_P(X) ({ long double __x = (X); \
 165                            !isinfl (__x) && !isnanl (__x); })
 166 #  else
 167 #   ifdef _GLIBCPP_HAVE_QFINITE
 168 #    define FINITEL_P(X) qfinite (X)
 169 #   else
 170 #    define FINITEL_P(X) FINITE_P (X)
 171 #   endif
 172 #  endif
 173 # endif
 174 #endif
 175
 176 /* Test whether number is infinite.  */
 177 #ifdef isinf
 178 /* This is an ISO C 9x macro.  */
 179 # define INFINITE_P(X) isinf (X)
 180 # define INFINITEF_P(X) isinf (X)
 181 # define INFINITEL_P(X) isinf (X)
 182 #else
 183 # ifdef IsINF
 184 /* This is for Solaris, which does not have special macros for other
 185    types.  */
 186 #  define INFINITE_P(X) IsINF (X)
 187 #  define INFINITEF_P(X) IsINF (X)
 188 #  define INFINITEL_P(X) IsINF (X)
 189 # else
 190 #  if defined _GLIBCPP_HAVE_ISINF
 191 #   define INFINITE_P(X) isinf (X)
 192 #  else
 193 #   ifdef _GLIBCPP_HAVE_FPCLASS
 194 #    ifdef _FPCLASS_PINF
 195 /* Mingw defines _FPCLASS_PINF.  */
 196 #     define FP_PINF _FPCLASS_PINF
 197 #    endif
 198 /* This is for Irix and Mingw.  */
 199 #    define INFINITE_P(X) (fpclass (fabs (X)) == FP_PINF)
 200 #   else
 201 #    ifdef IS_INF
 202 /* This is for AIX.  */
 203 #     define INFINITE_P(X) ({ double __d = (X); IS_INF (__d); })
 204 #    else
 205 #     error "We need INFINITE_P"
 206 #    endif
 207 #   endif
 208 #  endif
 209 #  if defined _GLIBCPP_HAVE_ISINFF
 210 #   define INFINITEF_P(X) isinff (X)
 211 #  else
 212 #   define INFINITEF_P(X) INFINITE_P (X)
 213 #  endif
 214 #  if defined _GLIBCPP_HAVE_ISINFL
 215 #   define INFINITEL_P(X) isinfl (X)
 216 #  else
 217 #   ifdef _GLIBCPP_HAVE_QFPCLASS
 218 #    define INFINITEL_P(X) (qfpclass (fabsl (X)) == FP_PINF)
 219 #   else
 220 #    define INFINITEL_P(X) INFINITE_P (X)
 221 #   endif
 222 #  endif
 223 # endif
 224 #endif
 225
 226 #ifdef __cplusplus
 227 extern "C" {
 228 #endif
 229
 230 #ifdef _GLIBCPP_HAVE_COMPLEX_H
 231 # include <complex.h>
 232 #else
 233 # include "complex-stub.h"
 234 #endif
 235 __complex__ double c_log (__complex__ double x);
 236 __complex__ float c_logf (__complex__ float x);
 237 __complex__ long double c_logl (__complex__ long double x);
 238
 239 /* signbit is a macro in ISO C 9x.  */
 240 #ifndef signbit
 241 extern int __signbitf (float);
 242 extern int __signbit (double);
 243 extern int __signbitl (long double);
 244
 245 # define signbit(x) \
 246      (sizeof (x) == sizeof (float) ?                                          \
 247         __signbitf (x)                                                        \
 248       : sizeof (x) == sizeof (double) ?                                       \
 249         __signbit (x) : __signbitl (x))
 250 #endif
 251
 252 #if BYTE_ORDER == BIG_ENDIAN
 253 typedef union
 254 {
 255   double value;
 256   struct
 257   {
 258     U_int32_t msw;
 259     U_int32_t lsw;
 260   } parts;
 261 } ieee_double_shape_type;
 262 #endif
 263 #if BYTE_ORDER == LITTLE_ENDIAN
 264 typedef union
 265 {
 266   double value;
 267   struct
 268   {
 269     U_int32_t lsw;
 270     U_int32_t msw;
 271   } parts;
 272 } ieee_double_shape_type;
 273 #endif
 274 /* Get the more significant 32 bit int from a double.  */
 275 #define GET_HIGH_WORD(i,d)                                      \
 276 do {                                                            \
 277   ieee_double_shape_type gh_u;                                  \
 278   gh_u.value = (d);                                             \
 279   (i) = gh_u.parts.msw;                                         \
 280 } while (0)
 281
 282
 283 typedef union
 284 {
 285   float value;
 286   U_int32_t word;
 287 } ieee_float_shape_type;
 288 /* Get a 32 bit int from a float.  */
 289 #define GET_FLOAT_WORD(i,d)                                     \
 290 do {                                                            \
 291   ieee_float_shape_type gf_u;                                   \
 292   gf_u.value = (d);                                             \
 293   (i) = gf_u.word;                                              \
 294 } while (0)
 295
 296
 297 #if BYTE_ORDER == BIG_ENDIAN
 298 typedef union
 299 {
 300   long double value;
 301   struct
 302   {
 303     unsigned int sign_exponent:16;
 304     unsigned int empty:16;
 305     U_int32_t msw;
 306     U_int32_t lsw;
 307   } parts;
 308 } ieee_long_double_shape_type;
 309 #endif
 310 #if BYTE_ORDER == LITTLE_ENDIAN
 311 typedef union
 312 {
 313   long double value;
 314   struct
 315   {
 316     U_int32_t lsw;
 317     U_int32_t msw;
 318     unsigned int sign_exponent:16;
 319     unsigned int empty:16;
 320   } parts;
 321 } ieee_long_double_shape_type;
 322 #endif
 323 /* Get int from the exponent of a long double.  */
 324 #define GET_LDOUBLE_EXP(exp,d)                                  \
 325 do {                                                            \
 326   ieee_long_double_shape_type ge_u;                             \
 327   ge_u.value = (d);                                             \
 328   (exp) = ge_u.parts.sign_exponent;                             \
 329 } while (0)
 330
 331 #if BYTE_ORDER == BIG_ENDIAN
 332 typedef union
 333 {
 334   long double value;
 335   struct
 336   {
 337     U_int64_t msw;
 338     U_int64_t lsw;
 339   } parts64;
 340   struct
 341   {
 342     U_int32_t w0, w1, w2, w3;
 343   } parts32;
 344 } ieee_quad_double_shape_type;
 345 #endif
 346 #if BYTE_ORDER == LITTLE_ENDIAN
 347 typedef union
 348 {
 349   long double value;
 350   struct
 351   {
 352     U_int64_t lsw;
 353     U_int64_t msw;
 354   } parts64;
 355   struct
 356   {
 357     U_int32_t w3, w2, w1, w0;
 358   } parts32;
 359 } ieee_quad_double_shape_type;
 360 #endif
 361 /* Get most significant 64 bit int from a quad long double.  */
 362 #define GET_LDOUBLE_MSW64(msw,d)                                \
 363 do {                                                            \
 364   ieee_quad_double_shape_type qw_u;                             \
 365   qw_u.value = (d);                                             \
 366   (ix0) = qw_u.parts64.msw;                                     \
 367 } while (0)
 368
 369
 370 /* Replacement for non-existing float functions.  */
 371 #if !defined(_GLIBCPP_HAVE_FABSF) && !defined(_GLIBCPP_HAVE___BUILTIN_FABSF)
 372 # define fabsf(x) fabs (x)
 373 #endif
 374 #if !defined(_GLIBCPP_HAVE_COSF) && !defined(_GLIBCPP_HAVE___BUILTIN_COSF)
 375 # define cosf(x) cos (x)
 376 #endif
 377 #ifndef _GLIBCPP_HAVE_COSHF
 378 # define coshf(x) cosh (x)
 379 #endif
 380 #ifndef _GLIBCPP_HAVE_EXPF
 381 # define expf(x) expf (x)
 382 #endif
 383 #ifndef _GLIBCPP_HAVE_LOGF
 384 # define logf(x) log(x)
 385 #endif
 386 #ifndef _GLIBCPP_HAVE_LOG10F
 387 # define log10f(x) log10 (x)
 388 #endif
 389 #ifndef _GLIBCPP_HAVE_POWF
 390 # define powf(x, y) pow (x, y)
 391 #endif
 392 #if !defined(_GLIBCPP_HAVE_SINF) && !defined(_GLIBCPP_HAVE___BUILTIN_SINF)
 393 # define sinf(x) sin (x)
 394 #endif
 395 #ifndef _GLIBCPP_HAVE_SINHF
 396 # define sinhf(x) sinh (x)
 397 #endif
 398 #if !defined(_GLIBCPP_HAVE_SQRTF) && !defined(_GLIBCPP_HAVE___BUILTIN_SQRTF)
 399 # define sqrtf(x) sqrt (x)
 400 #endif
 401 #ifndef _GLIBCPP_HAVE_TANF
 402 # define tanf(x) tan (x)
 403 #endif
 404 #ifndef _GLIBCPP_HAVE_TANHF
 405 # define tanhf(x) tanh (x)
 406 #endif
 407 #ifndef _GLIBCPP_HAVE_STRTOF
 408 # define strtof(s, e) strtod (s, e)
 409 #endif
 410
 411 #ifdef __cplusplus
 412 }
 413 #endif