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
73typedef unsigned int U_int32_t __attribute ((mode (SI)));
74typedef int Int32_t __attribute ((mode (SI)));
75typedef unsigned int U_int64_t __attribute ((mode (DI)));
76typedef 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
131extern "C" {
132#endif
133
134/* signbit is a macro in ISO C99.  */
135#ifndef signbit
136extern int __signbitf (float);
137extern int __signbit (double);
138extern 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
148typedef 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
159typedef 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)                                      \
171do {                                                            \
172  ieee_double_shape_type gh_u;                                  \
173  gh_u.value = (d);                                             \
174  (i) = gh_u.parts.msw;                                         \
175} while (0)
176
177
178typedef 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)                                     \
185do {                                                            \
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
193typedef 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
206typedef 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)                                  \
220do {                                                            \
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
227typedef 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
242typedef 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)				\
258do {								\
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
310