1 2/* @(#)fdlibm.h 1.5 04/04/22 */ 3/* 4 * ==================================================== 5 * Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved. 6 * 7 * Permission to use, copy, modify, and distribute this 8 * software is freely granted, provided that this notice 9 * is preserved. 10 * ==================================================== 11 */ 12 13/* Sometimes it's necessary to define __LITTLE_ENDIAN explicitly 14 but these catch some common cases. */ 15 16#if defined(i386) || defined(i486) || \ 17 defined(intel) || defined(x86) || defined(i86pc) || \ 18 defined(__alpha) || defined(__osf__) 19#define __LITTLE_ENDIAN 20#endif 21 22#ifdef __LITTLE_ENDIAN 23#define __HI(x) *(1+(int*)&x) 24#define __LO(x) *(int*)&x 25#define __HIp(x) *(1+(int*)x) 26#define __LOp(x) *(int*)x 27#else 28#define __HI(x) *(int*)&x 29#define __LO(x) *(1+(int*)&x) 30#define __HIp(x) *(int*)x 31#define __LOp(x) *(1+(int*)x) 32#endif 33 34#ifdef __STDC__ 35#define __P(p) p 36#else 37#define __P(p) () 38#endif 39 40/* 41 * ANSI/POSIX 42 */ 43 44extern int signgam; 45 46#define MAXFLOAT ((float)3.40282346638528860e+38) 47 48enum fdversion {fdlibm_ieee = -1, fdlibm_svid, fdlibm_xopen, fdlibm_posix}; 49 50#define _LIB_VERSION_TYPE enum fdversion 51#define _LIB_VERSION _fdlib_version 52 53/* if global variable _LIB_VERSION is not desirable, one may 54 * change the following to be a constant by: 55 * #define _LIB_VERSION_TYPE const enum version 56 * In that case, after one initializes the value _LIB_VERSION (see 57 * s_lib_version.c) during compile time, it cannot be modified 58 * in the middle of a program 59 */ 60extern _LIB_VERSION_TYPE _LIB_VERSION; 61 62#define _IEEE_ fdlibm_ieee 63#define _SVID_ fdlibm_svid 64#define _XOPEN_ fdlibm_xopen 65#define _POSIX_ fdlibm_posix 66 67struct exception { 68 int type; 69 char *name; 70 double arg1; 71 double arg2; 72 double retval; 73}; 74 75#define HUGE MAXFLOAT 76 77/* 78 * set X_TLOSS = pi*2**52, which is possibly defined in <values.h> 79 * (one may replace the following line by "#include <values.h>") 80 */ 81 82#define X_TLOSS 1.41484755040568800000e+16 83 84#define DOMAIN 1 85#define SING 2 86#define OVERFLOW 3 87#define UNDERFLOW 4 88#define TLOSS 5 89#define PLOSS 6 90 91/* 92 * ANSI/POSIX 93 */ 94extern double ieee_acos __P((double)); 95extern double ieee_asin __P((double)); 96extern double ieee_atan __P((double)); 97extern double ieee_atan2 __P((double, double)); 98extern double ieee_cos __P((double)); 99extern double ieee_sin __P((double)); 100extern double ieee_tan __P((double)); 101 102extern double ieee_cosh __P((double)); 103extern double ieee_sinh __P((double)); 104extern double ieee_tanh __P((double)); 105 106extern double ieee_exp __P((double)); 107extern double ieee_frexp __P((double, int *)); 108extern double ieee_ldexp __P((double, int)); 109extern double ieee_log __P((double)); 110extern double ieee_log10 __P((double)); 111extern double ieee_modf __P((double, double *)); 112 113extern double ieee_pow __P((double, double)); 114extern double ieee_sqrt __P((double)); 115 116extern double ieee_ceil __P((double)); 117extern double ieee_fabs __P((double)); 118extern double ieee_floor __P((double)); 119extern double ieee_fmod __P((double, double)); 120 121extern double ieee_erf __P((double)); 122extern double ieee_erfc __P((double)); 123extern double ieee_gamma __P((double)); 124extern double ieee_hypot __P((double, double)); 125extern int ieee_isnan __P((double)); 126extern int ieee_finite __P((double)); 127extern double ieee_j0 __P((double)); 128extern double ieee_j1 __P((double)); 129extern double ieee_jn __P((int, double)); 130extern double ieee_lgamma __P((double)); 131extern double ieee_y0 __P((double)); 132extern double ieee_y1 __P((double)); 133extern double ieee_yn __P((int, double)); 134 135extern double acosh __P((double)); 136extern double asinh __P((double)); 137extern double atanh __P((double)); 138extern double ieee_cbrt __P((double)); 139extern double ieee_logb __P((double)); 140extern double ieee_nextafter __P((double, double)); 141extern double ieee_remainder __P((double, double)); 142#ifdef _SCALB_INT 143extern double ieee_scalb __P((double, int)); 144#else 145extern double ieee_scalb __P((double, double)); 146#endif 147 148extern int ieee_matherr __P((struct exception *)); 149 150/* 151 * IEEE Test Vector 152 */ 153extern double ieee_significand __P((double)); 154 155/* 156 * Functions callable from C, intended to support IEEE arithmetic. 157 */ 158extern double ieee_copysign __P((double, double)); 159extern int ieee_ilogb __P((double)); 160extern double ieee_rint __P((double)); 161extern double ieee_scalbn __P((double, int)); 162 163/* 164 * BSD math library entry points 165 */ 166extern double ieee_expm1 __P((double)); 167extern double ieee_log1p __P((double)); 168 169/* 170 * Reentrant version of gamma & lgamma; passes signgam back by reference 171 * as the second argument; user must allocate space for signgam. 172 */ 173#ifdef _REENTRANT 174extern double ieee_gamma_r __P((double, int *)); 175extern double ieee_lgamma_r __P((double, int *)); 176#endif /* _REENTRANT */ 177 178/* ieee style elementary functions */ 179extern double __ieee754_sqrt __P((double)); 180extern double __ieee754_acos __P((double)); 181extern double __ieee754_acosh __P((double)); 182extern double __ieee754_log __P((double)); 183extern double __ieee754_atanh __P((double)); 184extern double __ieee754_asin __P((double)); 185extern double __ieee754_atan2 __P((double,double)); 186extern double __ieee754_exp __P((double)); 187extern double __ieee754_cosh __P((double)); 188extern double __ieee754_fmod __P((double,double)); 189extern double __ieee754_pow __P((double,double)); 190extern double __ieee754_lgamma_r __P((double,int *)); 191extern double __ieee754_gamma_r __P((double,int *)); 192extern double __ieee754_lgamma __P((double)); 193extern double __ieee754_gamma __P((double)); 194extern double __ieee754_log10 __P((double)); 195extern double __ieee754_sinh __P((double)); 196extern double __ieee754_hypot __P((double,double)); 197extern double __ieee754_j0 __P((double)); 198extern double __ieee754_j1 __P((double)); 199extern double __ieee754_y0 __P((double)); 200extern double __ieee754_y1 __P((double)); 201extern double __ieee754_jn __P((int,double)); 202extern double __ieee754_yn __P((int,double)); 203extern double __ieee754_remainder __P((double,double)); 204extern int __ieee754_rem_pio2 __P((double,double*)); 205#ifdef _SCALB_INT 206extern double __ieee754_scalb __P((double,int)); 207#else 208extern double __ieee754_scalb __P((double,double)); 209#endif 210 211/* fdlibm kernel function */ 212extern double __kernel_standard __P((double,double,int)); 213extern double __kernel_sin __P((double,double,int)); 214extern double __kernel_cos __P((double,double)); 215extern double __kernel_tan __P((double,double,int)); 216extern int __kernel_rem_pio2 __P((double*,double*,int,int,int,const int*)); 217