xref: /netbsd-src/include/math.h (revision 4472dbe5e3bd91ef2540bada7a7ca7384627ff9b)
1 /*	$NetBSD: math.h,v 1.21 2000/01/04 14:20:05 kleink Exp $	*/
2 
3 /*
4  * ====================================================
5  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6  *
7  * Developed at SunPro, a Sun Microsystems, Inc. business.
8  * Permission to use, copy, modify, and distribute this
9  * software is freely granted, provided that this notice
10  * is preserved.
11  * ====================================================
12  */
13 
14 /*
15  * @(#)fdlibm.h 5.1 93/09/24
16  */
17 
18 #ifndef _MATH_H_
19 #define _MATH_H_
20 
21 #include <sys/cdefs.h>
22 #include <sys/featuretest.h>
23 #include <machine/math.h>
24 
25 /*
26  * ANSI/POSIX
27  */
28 extern __const char __infinity[];
29 #define HUGE_VAL	(*(__const double *)(__const void *)__infinity)
30 
31 /*
32  * XOPEN/SVID
33  */
34 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) || \
35     defined(_XOPEN_SOURCE)
36 #define	M_E		2.7182818284590452354	/* e */
37 #define	M_LOG2E		1.4426950408889634074	/* log 2e */
38 #define	M_LOG10E	0.43429448190325182765	/* log 10e */
39 #define	M_LN2		0.69314718055994530942	/* log e2 */
40 #define	M_LN10		2.30258509299404568402	/* log e10 */
41 #define	M_PI		3.14159265358979323846	/* pi */
42 #define	M_PI_2		1.57079632679489661923	/* pi/2 */
43 #define	M_PI_4		0.78539816339744830962	/* pi/4 */
44 #define	M_1_PI		0.31830988618379067154	/* 1/pi */
45 #define	M_2_PI		0.63661977236758134308	/* 2/pi */
46 #define	M_2_SQRTPI	1.12837916709551257390	/* 2/sqrt(pi) */
47 #define	M_SQRT2		1.41421356237309504880	/* sqrt(2) */
48 #define	M_SQRT1_2	0.70710678118654752440	/* 1/sqrt(2) */
49 
50 #define	MAXFLOAT	((float)3.40282346638528860e+38)
51 extern int signgam;
52 #endif /* !_ANSI_SOURCE && !_POSIX_C_SOURCE || _XOPEN_SOURCE */
53 
54 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) && \
55     !defined(_XOPEN_SOURCE)
56 enum 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  */
68 extern  _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 
75 #ifndef __cplusplus
76 struct exception {
77 	int type;
78 	char *name;
79 	double arg1;
80 	double arg2;
81 	double retval;
82 };
83 #endif
84 
85 #define	HUGE		MAXFLOAT
86 
87 /*
88  * set X_TLOSS = pi*2**52, which is possibly defined in <values.h>
89  * (one may replace the following line by "#include <values.h>")
90  */
91 
92 #define X_TLOSS		1.41484755040568800000e+16
93 
94 #define	DOMAIN		1
95 #define	SING		2
96 #define	OVERFLOW	3
97 #define	UNDERFLOW	4
98 #define	TLOSS		5
99 #define	PLOSS		6
100 
101 #endif /* !_ANSI_SOURCE && !_POSIX_C_SOURCE && !_XOPEN_SOURCE */
102 
103 __BEGIN_DECLS
104 /*
105  * ANSI/POSIX
106  */
107 extern double acos __P((double));
108 extern double asin __P((double));
109 extern double atan __P((double));
110 extern double atan2 __P((double, double));
111 extern double cos __P((double));
112 extern double sin __P((double));
113 extern double tan __P((double));
114 
115 extern double cosh __P((double));
116 extern double sinh __P((double));
117 extern double tanh __P((double));
118 
119 extern double exp __P((double));
120 extern double frexp __P((double, int *));
121 extern double ldexp __P((double, int));
122 extern double log __P((double));
123 extern double log10 __P((double));
124 extern double modf __P((double, double *));
125 
126 extern double pow __P((double, double));
127 extern double sqrt __P((double));
128 
129 extern double ceil __P((double));
130 extern double fabs __P((double));
131 extern double floor __P((double));
132 extern double fmod __P((double, double));
133 
134 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) || \
135     defined(_XOPEN_SOURCE)
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 #if !defined(_XOPEN_SOURCE) || (_XOPEN_SOURCE - 0) >= 500
151 extern double acosh __P((double));
152 extern double asinh __P((double));
153 extern double atanh __P((double));
154 extern double cbrt __P((double));
155 extern double expm1 __P((double));
156 extern int ilogb __P((double));
157 extern double log1p __P((double));
158 extern double logb __P((double));
159 extern double nextafter __P((double, double));
160 extern double remainder __P((double, double));
161 extern double rint __P((double));
162 extern double scalb __P((double, double));
163 #endif /* !defined(_XOPEN_SOURCE) || (_XOPEN_SOURCE - 0) >= 500 */
164 #endif /* !_ANSI_SOURCE) && !_POSIX_C_SOURCE || _XOPEN_SOURCE */
165 
166 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) && \
167     !defined(_XOPEN_SOURCE)
168 #ifndef __cplusplus
169 extern int matherr __P((struct exception *));
170 #endif
171 
172 /*
173  * IEEE Test Vector
174  */
175 extern double significand __P((double));
176 
177 /*
178  * Functions callable from C, intended to support IEEE arithmetic.
179  */
180 extern double copysign __P((double, double));
181 extern double scalbn __P((double, int));
182 
183 /*
184  * BSD math library entry points
185  */
186 extern double cabs __P((/* struct complex { double r; double i; } */));
187 extern double drem __P((double, double));
188 
189 #endif /* !_ANSI_SOURCE && !_POSIX_C_SOURCE && !_XOPEN_SOURCE */
190 
191 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) && \
192     !defined(_XOPEN_SOURCE) || defined(_REENTRANT)
193 /*
194  * Reentrant version of gamma & lgamma; passes signgam back by reference
195  * as the second argument; user must allocate space for signgam.
196  */
197 extern double gamma_r __P((double, int *));
198 extern double lgamma_r __P((double, int *));
199 #endif /* !... || _REENTRANT */
200 
201 
202 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) && \
203     !defined(_XOPEN_SOURCE)
204 extern int isinf __P((double));
205 
206 /* float versions of ANSI/POSIX functions */
207 extern float acosf __P((float));
208 extern float asinf __P((float));
209 extern float atanf __P((float));
210 extern float atan2f __P((float, float));
211 extern float cosf __P((float));
212 extern float sinf __P((float));
213 extern float tanf __P((float));
214 
215 extern float coshf __P((float));
216 extern float sinhf __P((float));
217 extern float tanhf __P((float));
218 
219 extern float expf __P((float));
220 extern float frexpf __P((float, int *));
221 extern float ldexpf __P((float, int));
222 extern float logf __P((float));
223 extern float log10f __P((float));
224 extern float modff __P((float, float *));
225 
226 extern float powf __P((float, float));
227 extern float sqrtf __P((float));
228 
229 extern float ceilf __P((float));
230 extern float fabsf __P((float));
231 extern float floorf __P((float));
232 extern float fmodf __P((float, float));
233 
234 extern float erff __P((float));
235 extern float erfcf __P((float));
236 extern float gammaf __P((float));
237 extern float hypotf __P((float, float));
238 extern int isinff __P((float));
239 extern int isnanf __P((float));
240 extern int finitef __P((float));
241 extern float j0f __P((float));
242 extern float j1f __P((float));
243 extern float jnf __P((int, float));
244 extern float lgammaf __P((float));
245 extern float y0f __P((float));
246 extern float y1f __P((float));
247 extern float ynf __P((int, float));
248 
249 extern float acoshf __P((float));
250 extern float asinhf __P((float));
251 extern float atanhf __P((float));
252 extern float cbrtf __P((float));
253 extern float logbf __P((float));
254 extern float nextafterf __P((float, float));
255 extern float remainderf __P((float, float));
256 extern float scalbf __P((float, float));
257 
258 /*
259  * float version of IEEE Test Vector
260  */
261 extern float significandf __P((float));
262 
263 /*
264  * Float versions of functions callable from C, intended to support
265  * IEEE arithmetic.
266  */
267 extern float copysignf __P((float, float));
268 extern int ilogbf __P((float));
269 extern float rintf __P((float));
270 extern float scalbnf __P((float, int));
271 
272 /*
273  * float versions of BSD math library entry points
274  */
275 extern float cabsf __P((/* struct complex { float r; float i; } */));
276 extern float dremf __P((float, float));
277 extern float expm1f __P((float));
278 extern float log1pf __P((float));
279 #endif /* !_ANSI_SOURCE && !_POSIX_C_SOURCE && !_XOPEN_SOURCE */
280 
281 #if !defined(_ANSI_SOURCE) && !defined(_POSIX_C_SOURCE) && \
282     !defined(_XOPEN_SOURCE) || defined(_REENTRANT)
283 /*
284  * Float versions of reentrant version of gamma & lgamma; passes
285  * signgam back by reference as the second argument; user must
286  * allocate space for signgam.
287  */
288 extern float gammaf_r __P((float, int *));
289 extern float lgammaf_r __P((float, int *));
290 #endif /* !... || _REENTRANT */
291 
292 __END_DECLS
293 
294 #endif /* _MATH_H_ */
295