2 * ====================================================
3 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
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
9 * ====================================================
13 * from: @(#)fdlibm.h 5.1 93/09/24
17 #ifndef _MATH_PRIVATE_H_
18 #define _MATH_PRIVATE_H_
20 #include <machine/endian.h>
21 #include <sys/types.h>
23 /* The original fdlibm code used statements like:
24 n0 = ((*(int*)&one)>>29)^1; * index of high word *
25 ix0 = *(n0+(int*)&x); * high word of x *
26 ix1 = *((1-n0)+(int*)&x); * low word of x *
27 to dig two 32 bit words out of the 64 bit IEEE floating point
28 value. That is non-ANSI, and, moreover, the gcc instruction
29 scheduler gets it wrong. We instead use the following macros.
30 Unlike the original code, we determine the endianness at compile
31 time, not at run time; I don't see much benefit to selecting
32 endianness at run time. */
34 /* A union which permits us to convert between a double and two 32 bit
37 #if BYTE_ORDER == BIG_ENDIAN
47 } ieee_double_shape_type
;
51 #if BYTE_ORDER == LITTLE_ENDIAN
61 } ieee_double_shape_type
;
65 /* Get two 32 bit ints from a double. */
67 #define EXTRACT_WORDS(ix0,ix1,d) \
69 ieee_double_shape_type ew_u; \
71 (ix0) = ew_u.parts.msw; \
72 (ix1) = ew_u.parts.lsw; \
75 /* Get the more significant 32 bit int from a double. */
77 #define GET_HIGH_WORD(i,d) \
79 ieee_double_shape_type gh_u; \
81 (i) = gh_u.parts.msw; \
84 /* Get the less significant 32 bit int from a double. */
86 #define GET_LOW_WORD(i,d) \
88 ieee_double_shape_type gl_u; \
90 (i) = gl_u.parts.lsw; \
93 /* Set a double from two 32 bit ints. */
95 #define INSERT_WORDS(d,ix0,ix1) \
97 ieee_double_shape_type iw_u; \
98 iw_u.parts.msw = (ix0); \
99 iw_u.parts.lsw = (ix1); \
103 /* Set the more significant 32 bits of a double from an int. */
105 #define SET_HIGH_WORD(d,v) \
107 ieee_double_shape_type sh_u; \
109 sh_u.parts.msw = (v); \
113 /* Set the less significant 32 bits of a double from an int. */
115 #define SET_LOW_WORD(d,v) \
117 ieee_double_shape_type sl_u; \
119 sl_u.parts.lsw = (v); \
123 /* A union which permits us to convert between a float and a 32 bit
130 } ieee_float_shape_type
;
132 /* Get a 32 bit int from a float. */
134 #define GET_FLOAT_WORD(i,d) \
136 ieee_float_shape_type gf_u; \
141 /* Set a float from a 32 bit int. */
143 #define SET_FLOAT_WORD(d,i) \
145 ieee_float_shape_type sf_u; \
150 /* ieee style elementary functions */
151 extern double __ieee754_sqrt
__P((double));
152 extern double __ieee754_acos
__P((double));
153 extern double __ieee754_acosh
__P((double));
154 extern double __ieee754_log
__P((double));
155 extern double __ieee754_atanh
__P((double));
156 extern double __ieee754_asin
__P((double));
157 extern double __ieee754_atan2
__P((double,double));
158 extern double __ieee754_exp
__P((double));
159 extern double __ieee754_cosh
__P((double));
160 extern double __ieee754_fmod
__P((double,double));
161 extern double __ieee754_pow
__P((double,double));
162 extern double __ieee754_lgamma_r
__P((double,int *));
163 extern double __ieee754_gamma_r
__P((double,int *));
164 extern double __ieee754_lgamma
__P((double));
165 extern double __ieee754_gamma
__P((double));
166 extern double __ieee754_log10
__P((double));
167 extern double __ieee754_sinh
__P((double));
168 extern double __ieee754_hypot
__P((double,double));
169 extern double __ieee754_j0
__P((double));
170 extern double __ieee754_j1
__P((double));
171 extern double __ieee754_y0
__P((double));
172 extern double __ieee754_y1
__P((double));
173 extern double __ieee754_jn
__P((int,double));
174 extern double __ieee754_yn
__P((int,double));
175 extern double __ieee754_remainder
__P((double,double));
176 extern int __ieee754_rem_pio2
__P((double,double*));
177 extern double __ieee754_scalb
__P((double,double));
179 /* fdlibm kernel function */
180 extern double __kernel_standard
__P((double,double,int));
181 extern double __kernel_sin
__P((double,double,int));
182 extern double __kernel_cos
__P((double,double));
183 extern double __kernel_tan
__P((double,double,int));
184 extern int __kernel_rem_pio2
__P((double*,double*,int,int,int,const int*));
187 /* ieee style elementary float functions */
188 extern float __ieee754_sqrtf
__P((float));
189 extern float __ieee754_acosf
__P((float));
190 extern float __ieee754_acoshf
__P((float));
191 extern float __ieee754_logf
__P((float));
192 extern float __ieee754_atanhf
__P((float));
193 extern float __ieee754_asinf
__P((float));
194 extern float __ieee754_atan2f
__P((float,float));
195 extern float __ieee754_expf
__P((float));
196 extern float __ieee754_coshf
__P((float));
197 extern float __ieee754_fmodf
__P((float,float));
198 extern float __ieee754_powf
__P((float,float));
199 extern float __ieee754_lgammaf_r
__P((float,int *));
200 extern float __ieee754_gammaf_r
__P((float,int *));
201 extern float __ieee754_lgammaf
__P((float));
202 extern float __ieee754_gammaf
__P((float));
203 extern float __ieee754_log10f
__P((float));
204 extern float __ieee754_sinhf
__P((float));
205 extern float __ieee754_hypotf
__P((float,float));
206 extern float __ieee754_j0f
__P((float));
207 extern float __ieee754_j1f
__P((float));
208 extern float __ieee754_y0f
__P((float));
209 extern float __ieee754_y1f
__P((float));
210 extern float __ieee754_jnf
__P((int,float));
211 extern float __ieee754_ynf
__P((int,float));
212 extern float __ieee754_remainderf
__P((float,float));
213 extern int __ieee754_rem_pio2f
__P((float,float*));
214 extern float __ieee754_scalbf
__P((float,float));
216 /* float versions of fdlibm kernel functions */
217 extern float __kernel_sinf
__P((float,float,int));
218 extern float __kernel_cosf
__P((float,float));
219 extern float __kernel_tanf
__P((float,float,int));
220 extern int __kernel_rem_pio2f
__P((float*,float*,int,int,int,const int*));
222 #endif /* _MATH_PRIVATE_H_ */