J. Duke | 319a3b9 | 2007-12-01 00:00:00 +0000 | [diff] [blame^] | 1 | |
| 2 | /* |
| 3 | * Copyright 1998-2002 Sun Microsystems, Inc. All Rights Reserved. |
| 4 | * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER. |
| 5 | * |
| 6 | * This code is free software; you can redistribute it and/or modify it |
| 7 | * under the terms of the GNU General Public License version 2 only, as |
| 8 | * published by the Free Software Foundation. Sun designates this |
| 9 | * particular file as subject to the "Classpath" exception as provided |
| 10 | * by Sun in the LICENSE file that accompanied this code. |
| 11 | * |
| 12 | * This code is distributed in the hope that it will be useful, but WITHOUT |
| 13 | * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
| 14 | * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| 15 | * version 2 for more details (a copy is included in the LICENSE file that |
| 16 | * accompanied this code). |
| 17 | * |
| 18 | * You should have received a copy of the GNU General Public License version |
| 19 | * 2 along with this work; if not, write to the Free Software Foundation, |
| 20 | * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. |
| 21 | * |
| 22 | * Please contact Sun Microsystems, Inc., 4150 Network Circle, Santa Clara, |
| 23 | * CA 95054 USA or visit www.sun.com if you need additional information or |
| 24 | * have any questions. |
| 25 | */ |
| 26 | |
| 27 | #include "jfdlibm.h" |
| 28 | |
| 29 | #ifdef __NEWVALID /* special setup for Sun test regime */ |
| 30 | #if defined(i386) || defined(i486) || \ |
| 31 | defined(intel) || defined(x86) || \ |
| 32 | defined(i86pc) || defined(_M_IA64) |
| 33 | #define _LITTLE_ENDIAN |
| 34 | #endif |
| 35 | #endif |
| 36 | |
| 37 | #ifdef _LITTLE_ENDIAN |
| 38 | #define __HI(x) *(1+(int*)&x) |
| 39 | #define __LO(x) *(int*)&x |
| 40 | #define __HIp(x) *(1+(int*)x) |
| 41 | #define __LOp(x) *(int*)x |
| 42 | #else |
| 43 | #define __HI(x) *(int*)&x |
| 44 | #define __LO(x) *(1+(int*)&x) |
| 45 | #define __HIp(x) *(int*)x |
| 46 | #define __LOp(x) *(1+(int*)x) |
| 47 | #endif |
| 48 | |
| 49 | #ifdef __STDC__ |
| 50 | #define __P(p) p |
| 51 | #else |
| 52 | #define __P(p) () |
| 53 | #endif |
| 54 | |
| 55 | /* |
| 56 | * ANSI/POSIX |
| 57 | */ |
| 58 | |
| 59 | extern int signgam; |
| 60 | |
| 61 | #define MAXFLOAT ((float)3.40282346638528860e+38) |
| 62 | |
| 63 | enum fdversion {fdlibm_ieee = -1, fdlibm_svid, fdlibm_xopen, fdlibm_posix}; |
| 64 | |
| 65 | #define _LIB_VERSION_TYPE enum fdversion |
| 66 | #define _LIB_VERSION _fdlib_version |
| 67 | |
| 68 | /* if global variable _LIB_VERSION is not desirable, one may |
| 69 | * change the following to be a constant by: |
| 70 | * #define _LIB_VERSION_TYPE const enum version |
| 71 | * In that case, after one initializes the value _LIB_VERSION (see |
| 72 | * s_lib_version.c) during compile time, it cannot be modified |
| 73 | * in the middle of a program |
| 74 | */ |
| 75 | extern _LIB_VERSION_TYPE _LIB_VERSION; |
| 76 | |
| 77 | #define _IEEE_ fdlibm_ieee |
| 78 | #define _SVID_ fdlibm_svid |
| 79 | #define _XOPEN_ fdlibm_xopen |
| 80 | #define _POSIX_ fdlibm_posix |
| 81 | |
| 82 | struct exception { |
| 83 | int type; |
| 84 | char *name; |
| 85 | double arg1; |
| 86 | double arg2; |
| 87 | double retval; |
| 88 | }; |
| 89 | |
| 90 | #define HUGE MAXFLOAT |
| 91 | |
| 92 | /* |
| 93 | * set X_TLOSS = pi*2**52, which is possibly defined in <values.h> |
| 94 | * (one may replace the following line by "#include <values.h>") |
| 95 | */ |
| 96 | |
| 97 | #define X_TLOSS 1.41484755040568800000e+16 |
| 98 | |
| 99 | #define DOMAIN 1 |
| 100 | #define SING 2 |
| 101 | #define OVERFLOW 3 |
| 102 | #define UNDERFLOW 4 |
| 103 | #define TLOSS 5 |
| 104 | #define PLOSS 6 |
| 105 | |
| 106 | /* |
| 107 | * ANSI/POSIX |
| 108 | */ |
| 109 | extern double acos __P((double)); |
| 110 | extern double asin __P((double)); |
| 111 | extern double atan __P((double)); |
| 112 | extern double atan2 __P((double, double)); |
| 113 | extern double cos __P((double)); |
| 114 | extern double sin __P((double)); |
| 115 | extern double tan __P((double)); |
| 116 | |
| 117 | extern double cosh __P((double)); |
| 118 | extern double sinh __P((double)); |
| 119 | extern double tanh __P((double)); |
| 120 | |
| 121 | extern double exp __P((double)); |
| 122 | extern double frexp __P((double, int *)); |
| 123 | extern double ldexp __P((double, int)); |
| 124 | extern double log __P((double)); |
| 125 | extern double log10 __P((double)); |
| 126 | extern double modf __P((double, double *)); |
| 127 | |
| 128 | extern double pow __P((double, double)); |
| 129 | extern double sqrt __P((double)); |
| 130 | |
| 131 | extern double ceil __P((double)); |
| 132 | extern double fabs __P((double)); |
| 133 | extern double floor __P((double)); |
| 134 | extern double fmod __P((double, double)); |
| 135 | |
| 136 | extern double erf __P((double)); |
| 137 | extern double erfc __P((double)); |
| 138 | extern double gamma __P((double)); |
| 139 | extern double hypot __P((double, double)); |
| 140 | extern int isnan __P((double)); |
| 141 | extern int finite __P((double)); |
| 142 | extern double j0 __P((double)); |
| 143 | extern double j1 __P((double)); |
| 144 | extern double jn __P((int, double)); |
| 145 | extern double lgamma __P((double)); |
| 146 | extern double y0 __P((double)); |
| 147 | extern double y1 __P((double)); |
| 148 | extern double yn __P((int, double)); |
| 149 | |
| 150 | extern double acosh __P((double)); |
| 151 | extern double asinh __P((double)); |
| 152 | extern double atanh __P((double)); |
| 153 | extern double cbrt __P((double)); |
| 154 | extern double logb __P((double)); |
| 155 | extern double nextafter __P((double, double)); |
| 156 | extern double remainder __P((double, double)); |
| 157 | #ifdef _SCALB_INT |
| 158 | extern double scalb __P((double, int)); |
| 159 | #else |
| 160 | extern double scalb __P((double, double)); |
| 161 | #endif |
| 162 | |
| 163 | extern int matherr __P((struct exception *)); |
| 164 | |
| 165 | /* |
| 166 | * IEEE Test Vector |
| 167 | */ |
| 168 | extern double significand __P((double)); |
| 169 | |
| 170 | /* |
| 171 | * Functions callable from C, intended to support IEEE arithmetic. |
| 172 | */ |
| 173 | extern double copysign __P((double, double)); |
| 174 | extern int ilogb __P((double)); |
| 175 | extern double rint __P((double)); |
| 176 | extern double scalbn __P((double, int)); |
| 177 | |
| 178 | /* |
| 179 | * BSD math library entry points |
| 180 | */ |
| 181 | extern double expm1 __P((double)); |
| 182 | extern double log1p __P((double)); |
| 183 | |
| 184 | /* |
| 185 | * Reentrant version of gamma & lgamma; passes signgam back by reference |
| 186 | * as the second argument; user must allocate space for signgam. |
| 187 | */ |
| 188 | #ifdef _REENTRANT |
| 189 | extern double gamma_r __P((double, int *)); |
| 190 | extern double lgamma_r __P((double, int *)); |
| 191 | #endif /* _REENTRANT */ |
| 192 | |
| 193 | /* ieee style elementary functions */ |
| 194 | extern double __ieee754_sqrt __P((double)); |
| 195 | extern double __ieee754_acos __P((double)); |
| 196 | extern double __ieee754_acosh __P((double)); |
| 197 | extern double __ieee754_log __P((double)); |
| 198 | extern double __ieee754_atanh __P((double)); |
| 199 | extern double __ieee754_asin __P((double)); |
| 200 | extern double __ieee754_atan2 __P((double,double)); |
| 201 | extern double __ieee754_exp __P((double)); |
| 202 | extern double __ieee754_cosh __P((double)); |
| 203 | extern double __ieee754_fmod __P((double,double)); |
| 204 | extern double __ieee754_pow __P((double,double)); |
| 205 | extern double __ieee754_lgamma_r __P((double,int *)); |
| 206 | extern double __ieee754_gamma_r __P((double,int *)); |
| 207 | extern double __ieee754_lgamma __P((double)); |
| 208 | extern double __ieee754_gamma __P((double)); |
| 209 | extern double __ieee754_log10 __P((double)); |
| 210 | extern double __ieee754_sinh __P((double)); |
| 211 | extern double __ieee754_hypot __P((double,double)); |
| 212 | extern double __ieee754_j0 __P((double)); |
| 213 | extern double __ieee754_j1 __P((double)); |
| 214 | extern double __ieee754_y0 __P((double)); |
| 215 | extern double __ieee754_y1 __P((double)); |
| 216 | extern double __ieee754_jn __P((int,double)); |
| 217 | extern double __ieee754_yn __P((int,double)); |
| 218 | extern double __ieee754_remainder __P((double,double)); |
| 219 | extern int __ieee754_rem_pio2 __P((double,double*)); |
| 220 | #ifdef _SCALB_INT |
| 221 | extern double __ieee754_scalb __P((double,int)); |
| 222 | #else |
| 223 | extern double __ieee754_scalb __P((double,double)); |
| 224 | #endif |
| 225 | |
| 226 | /* fdlibm kernel function */ |
| 227 | extern double __kernel_standard __P((double,double,int)); |
| 228 | extern double __kernel_sin __P((double,double,int)); |
| 229 | extern double __kernel_cos __P((double,double)); |
| 230 | extern double __kernel_tan __P((double,double,int)); |
| 231 | extern int __kernel_rem_pio2 __P((double*,double*,int,int,int,const int*)); |