blob: f375ad258e78405d35727f5b623408a234067bba [file] [log] [blame]
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -08001
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
Elliott Hughes26ea4292010-08-24 16:50:43 -070013#ifdef __cplusplus
14extern "C" {
15#endif
16
17#ifndef __LITTLE_ENDIAN
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -080018/* Sometimes it's necessary to define __LITTLE_ENDIAN explicitly
19 but these catch some common cases. */
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -080020#if defined(i386) || defined(i486) || \
21 defined(intel) || defined(x86) || defined(i86pc) || \
Chris Dearman2b584d82011-04-14 17:11:20 -070022 defined(__alpha) || defined(__osf__) || defined(__MIPSEL__)
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -080023#define __LITTLE_ENDIAN
24#endif
Elliott Hughes26ea4292010-08-24 16:50:43 -070025#endif
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -080026
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
Elliott Hughes26ea4292010-08-24 16:50:43 -070039#ifndef __P
40#define __FDLIBM_P_DEFINED
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -080041#ifdef __STDC__
42#define __P(p) p
43#else
44#define __P(p) ()
45#endif
Elliott Hughes26ea4292010-08-24 16:50:43 -070046#endif
The Android Open Source Projectb07e1d92009-03-03 19:29:30 -080047
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*));
Elliott Hughes26ea4292010-08-24 16:50:43 -0700225
226#ifdef __FDLIBM_P_DEFINED
227#undef __P
228#endif
229
230#ifdef __cplusplus
231}; /* extern "C" */
232#endif