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