math/math_private.h

   1 /*
   2  * ====================================================
   3  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
   4  *
   5  * Developed at SunPro, a Sun Microsystems, Inc. business.
   6  * Permission to use, copy, modify, and distribute this
   7  * software is freely granted, provided that this notice
   8  * is preserved.
   9  * ====================================================
  10  */
  11
  12 /*
  13  * from: @(#)fdlibm.h 5.1 93/09/24
  14  */
  15
  16 #ifndef _MATH_PRIVATE_H_
  17 #define _MATH_PRIVATE_H_
  18
  19 #include <endian.h>
  20 #include <sys/types.h>
  21
  22 /* The original fdlibm code used statements like:
  23         n0 = ((*(int*)&one)>>29)^1;             * index of high word *
  24         ix0 = *(n0+(int*)&x);                   * high word of x *
  25         ix1 = *((1-n0)+(int*)&x);               * low word of x *
  26    to dig two 32 bit words out of the 64 bit IEEE floating point
  27    value.  That is non-ANSI, and, moreover, the gcc instruction
  28    scheduler gets it wrong.  We instead use the following macros.
  29    Unlike the original code, we determine the endianness at compile
  30    time, not at run time; I don't see much benefit to selecting
  31    endianness at run time.  */
  32
  33 /* A union which permits us to convert between a double and two 32 bit
  34    ints.  */
  35
  36 #if __FLOAT_WORD_ORDER == BIG_ENDIAN
  37
  38 typedef union
  39 {
  40   double value;
  41   struct
  42   {
  43     u_int32_t msw;
  44     u_int32_t lsw;
  45   } parts;
  46 } ieee_double_shape_type;
  47
  48 #endif
  49
  50 #if __FLOAT_WORD_ORDER == LITTLE_ENDIAN
  51
  52 typedef union
  53 {
  54   double value;
  55   struct
  56   {
  57     u_int32_t lsw;
  58     u_int32_t msw;
  59   } parts;
  60 } ieee_double_shape_type;
  61
  62 #endif
  63
  64 /* Get two 32 bit ints from a double.  */
  65
  66 #define EXTRACT_WORDS(ix0,ix1,d)                                \
  67 do {                                                            \
  68   ieee_double_shape_type ew_u;                                  \
  69   ew_u.value = (d);                                             \
  70   (ix0) = ew_u.parts.msw;                                       \
  71   (ix1) = ew_u.parts.lsw;                                       \
  72 } while (0)
  73
  74 /* Get the more significant 32 bit int from a double.  */
  75
  76 #define GET_HIGH_WORD(i,d)                                      \
  77 do {                                                            \
  78   ieee_double_shape_type gh_u;                                  \
  79   gh_u.value = (d);                                             \
  80   (i) = gh_u.parts.msw;                                         \
  81 } while (0)
  82
  83 /* Get the less significant 32 bit int from a double.  */
  84
  85 #define GET_LOW_WORD(i,d)                                       \
  86 do {                                                            \
  87   ieee_double_shape_type gl_u;                                  \
  88   gl_u.value = (d);                                             \
  89   (i) = gl_u.parts.lsw;                                         \
  90 } while (0)
  91
  92 /* Set a double from two 32 bit ints.  */
  93
  94 #define INSERT_WORDS(d,ix0,ix1)                                 \
  95 do {                                                            \
  96   ieee_double_shape_type iw_u;                                  \
  97   iw_u.parts.msw = (ix0);                                       \
  98   iw_u.parts.lsw = (ix1);                                       \
  99   (d) = iw_u.value;                                             \
 100 } while (0)
 101
 102 /* Set the more significant 32 bits of a double from an int.  */
 103
 104 #define SET_HIGH_WORD(d,v)                                      \
 105 do {                                                            \
 106   ieee_double_shape_type sh_u;                                  \
 107   sh_u.value = (d);                                             \
 108   sh_u.parts.msw = (v);                                         \
 109   (d) = sh_u.value;                                             \
 110 } while (0)
 111
 112 /* Set the less significant 32 bits of a double from an int.  */
 113
 114 #define SET_LOW_WORD(d,v)                                       \
 115 do {                                                            \
 116   ieee_double_shape_type sl_u;                                  \
 117   sl_u.value = (d);                                             \
 118   sl_u.parts.lsw = (v);                                         \
 119   (d) = sl_u.value;                                             \
 120 } while (0)
 121
 122 /* A union which permits us to convert between a float and a 32 bit
 123    int.  */
 124
 125 typedef union
 126 {
 127   float value;
 128   u_int32_t word;
 129 } ieee_float_shape_type;
 130
 131 /* Get a 32 bit int from a float.  */
 132
 133 #define GET_FLOAT_WORD(i,d)                                     \
 134 do {                                                            \
 135   ieee_float_shape_type gf_u;                                   \
 136   gf_u.value = (d);                                             \
 137   (i) = gf_u.word;                                              \
 138 } while (0)
 139
 140 /* Set a float from a 32 bit int.  */
 141
 142 #define SET_FLOAT_WORD(d,i)                                     \
 143 do {                                                            \
 144   ieee_float_shape_type sf_u;                                   \
 145   sf_u.word = (i);                                              \
 146   (d) = sf_u.value;                                             \
 147 } while (0)
 148
 149 /* Get long double macros from a separate header.  */
 150 #include <math_ldbl.h>
 151
 152 /* ieee style elementary functions */
 153 extern double __ieee754_sqrt (double);
 154 extern double __ieee754_acos (double);
 155 extern double __ieee754_acosh (double);
 156 extern double __ieee754_log (double);
 157 extern double __ieee754_atanh (double);
 158 extern double __ieee754_asin (double);
 159 extern double __ieee754_atan2 (double,double);
 160 extern double __ieee754_exp (double);
 161 extern double __ieee754_exp2 (double);
 162 extern double __ieee754_exp10 (double);
 163 extern double __ieee754_cosh (double);
 164 extern double __ieee754_fmod (double,double);
 165 extern double __ieee754_pow (double,double);
 166 extern double __ieee754_lgamma_r (double,int *);
 167 extern double __ieee754_gamma_r (double,int *);
 168 extern double __ieee754_lgamma (double);
 169 extern double __ieee754_gamma (double);
 170 extern double __ieee754_log10 (double);
 171 extern double __ieee754_log2 (double);
 172 extern double __ieee754_sinh (double);
 173 extern double __ieee754_hypot (double,double);
 174 extern double __ieee754_j0 (double);
 175 extern double __ieee754_j1 (double);
 176 extern double __ieee754_y0 (double);
 177 extern double __ieee754_y1 (double);
 178 extern double __ieee754_jn (int,double);
 179 extern double __ieee754_yn (int,double);
 180 extern double __ieee754_remainder (double,double);
 181 extern int32_t __ieee754_rem_pio2 (double,double*);
 182 extern double __ieee754_scalb (double,double);
 183
 184 /* fdlibm kernel function */
 185 extern double __kernel_standard (double,double,int);
 186 extern double __kernel_sin (double,double,int);
 187 extern double __kernel_cos (double,double);
 188 extern double __kernel_tan (double,double,int);
 189 extern int    __kernel_rem_pio2 (double*,double*,int,int,int, const int32_t*);
 190
 191 /* internal functions.  */
 192 extern double __copysign (double x, double __y);
 193
 194 #if __GNUC_PREREQ (4, 0)
 195 extern inline double __copysign (double x, double y)
 196 { return __builtin_copysign (x, y); }
 197 #endif
 198
 199 /* ieee style elementary float functions */
 200 extern float __ieee754_sqrtf (float);
 201 extern float __ieee754_acosf (float);
 202 extern float __ieee754_acoshf (float);
 203 extern float __ieee754_logf (float);
 204 extern float __ieee754_atanhf (float);
 205 extern float __ieee754_asinf (float);
 206 extern float __ieee754_atan2f (float,float);
 207 extern float __ieee754_expf (float);
 208 extern float __ieee754_exp2f (float);
 209 extern float __ieee754_exp10f (float);
 210 extern float __ieee754_coshf (float);
 211 extern float __ieee754_fmodf (float,float);
 212 extern float __ieee754_powf (float,float);
 213 extern float __ieee754_lgammaf_r (float,int *);
 214 extern float __ieee754_gammaf_r (float,int *);
 215 extern float __ieee754_lgammaf (float);
 216 extern float __ieee754_gammaf (float);
 217 extern float __ieee754_log10f (float);
 218 extern float __ieee754_log2f (float);
 219 extern float __ieee754_sinhf (float);
 220 extern float __ieee754_hypotf (float,float);
 221 extern float __ieee754_j0f (float);
 222 extern float __ieee754_j1f (float);
 223 extern float __ieee754_y0f (float);
 224 extern float __ieee754_y1f (float);
 225 extern float __ieee754_jnf (int,float);
 226 extern float __ieee754_ynf (int,float);
 227 extern float __ieee754_remainderf (float,float);
 228 extern int32_t __ieee754_rem_pio2f (float,float*);
 229 extern float __ieee754_scalbf (float,float);
 230
 231
 232 /* float versions of fdlibm kernel functions */
 233 extern float __kernel_sinf (float,float,int);
 234 extern float __kernel_cosf (float,float);
 235 extern float __kernel_tanf (float,float,int);
 236 extern int   __kernel_rem_pio2f (float*,float*,int,int,int, const int32_t*);
 237
 238 /* internal functions.  */
 239 extern float __copysignf (float x, float __y);
 240
 241 #if __GNUC_PREREQ (4, 0)
 242 extern inline float __copysignf (float x, float y)
 243 { return __builtin_copysignf (x, y); }
 244 #endif
 245
 246 /* ieee style elementary long double functions */
 247 extern long double __ieee754_sqrtl (long double);
 248 extern long double __ieee754_acosl (long double);
 249 extern long double __ieee754_acoshl (long double);
 250 extern long double __ieee754_logl (long double);
 251 extern long double __ieee754_atanhl (long double);
 252 extern long double __ieee754_asinl (long double);
 253 extern long double __ieee754_atan2l (long double,long double);
 254 extern long double __ieee754_expl (long double);
 255 extern long double __ieee754_exp2l (long double);
 256 extern long double __ieee754_exp10l (long double);
 257 extern long double __ieee754_coshl (long double);
 258 extern long double __ieee754_fmodl (long double,long double);
 259 extern long double __ieee754_powl (long double,long double);
 260 extern long double __ieee754_lgammal_r (long double,int *);
 261 extern long double __ieee754_gammal_r (long double,int *);
 262 extern long double __ieee754_lgammal (long double);
 263 extern long double __ieee754_gammal (long double);
 264 extern long double __ieee754_log10l (long double);
 265 extern long double __ieee754_log2l (long double);
 266 extern long double __ieee754_sinhl (long double);
 267 extern long double __ieee754_hypotl (long double,long double);
 268 extern long double __ieee754_j0l (long double);
 269 extern long double __ieee754_j1l (long double);
 270 extern long double __ieee754_y0l (long double);
 271 extern long double __ieee754_y1l (long double);
 272 extern long double __ieee754_jnl (int,long double);
 273 extern long double __ieee754_ynl (int,long double);
 274 extern long double __ieee754_remainderl (long double,long double);
 275 extern int   __ieee754_rem_pio2l (long double,long double*);
 276 extern long double __ieee754_scalbl (long double,long double);
 277
 278 /* long double versions of fdlibm kernel functions */
 279 extern long double __kernel_sinl (long double,long double,int);
 280 extern long double __kernel_cosl (long double,long double);
 281 extern long double __kernel_tanl (long double,long double,int);
 282 extern void __kernel_sincosl (long double,long double,
 283                               long double *,long double *, int);
 284 extern int   __kernel_rem_pio2l (long double*,long double*,int,int,
 285                                  int,const int*);
 286
 287 #ifndef NO_LONG_DOUBLE
 288 /* prototypes required to compile the ldbl-96 support without warnings */
 289 extern int __finitel (long double);
 290 extern int __ilogbl (long double);
 291 extern int __isinfl (long double);
 292 extern int __isnanl (long double);
 293 extern long double __atanl (long double);
 294 extern long double __copysignl (long double, long double);
 295 extern long double __expm1l (long double);
 296 extern long double __floorl (long double);
 297 extern long double __frexpl (long double, int *);
 298 extern long double __ldexpl (long double, int);
 299 extern long double __log1pl (long double);
 300 extern long double __nanl (const char *);
 301 extern long double __rintl (long double);
 302 extern long double __scalbnl (long double, int);
 303 extern long double __sqrtl (long double x);
 304 extern long double fabsl (long double x);
 305 extern void __sincosl (long double, long double *, long double *);
 306 extern long double __logbl (long double x);
 307 extern long double __significandl (long double x);
 308
 309 #if __GNUC_PREREQ (4, 0)
 310 extern inline long double __copysignl (long double x, long double y)
 311 { return __builtin_copysignl (x, y); }
 312 #endif
 313
 314 #endif
 315
 316 /* Prototypes for functions of the IBM Accurate Mathematical Library.  */
 317 extern double __exp1 (double __x, double __xx, double __error);
 318 extern double __sin (double __x);
 319 extern double __cos (double __x);
 320 extern int __branred (double __x, double *__a, double *__aa);
 321 extern void __doasin (double __x, double __dx, double __v[]);
 322 extern void __dubsin (double __x, double __dx, double __v[]);
 323 extern void __dubcos (double __x, double __dx, double __v[]);
 324 extern double __halfulp (double __x, double __y);
 325 extern double __sin32 (double __x, double __res, double __res1);
 326 extern double __cos32 (double __x, double __res, double __res1);
 327 extern double __mpsin (double __x, double __dx);
 328 extern double __mpcos (double __x, double __dx);
 329 extern double __mpsin1 (double __x);
 330 extern double __mpcos1 (double __x);
 331 extern double __slowexp (double __x);
 332 extern double __slowpow (double __x, double __y, double __z);
 333 extern void __docos (double __x, double __dx, double __v[]);
 334
 335 #ifndef math_opt_barrier
 336 #define math_opt_barrier(x) \
 337 ({ __typeof (x) __x = x; __asm ("" : "+m" (__x)); __x; })
 338 #define math_force_eval(x) __asm __volatile ("" : : "m" (x))
 339 #endif
 340
 341 #endif /* _MATH_PRIVATE_H_ */