1 /*
2 * CDDL HEADER START
3 *
4 * The contents of this file are subject to the terms of the
5 * Common Development and Distribution License (the "License").
6 * You may not use this file except in compliance with the License.
7 *
8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9 * or http://www.opensolaris.org/os/licensing.
10 * See the License for the specific language governing permissions
11 * and limitations under the License.
12 *
13 * When distributing Covered Code, include this CDDL HEADER in each
14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15 * If applicable, add the following below this CDDL HEADER, with the
16 * fields enclosed by brackets "[]" replaced with your own identifying
17 * information: Portions Copyright [yyyy] [name of copyright owner]
18 *
19 * CDDL HEADER END
20 */
21 /*
22 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
23 */
24 /*
25 * Copyright 2004 Sun Microsystems, Inc. All rights reserved.
26 * Use is subject to license terms.
27 */
28
29 #ifndef _MATH_H
30 #define _MATH_H
31
32 #include <iso/math_iso.h>
33 #include <iso/math_c99.h>
34
35 #if __cplusplus >= 199711L
36 using std::abs;
37 using std::acos;
38 using std::asin;
39 using std::atan2;
40 using std::atan;
41 using std::ceil;
42 using std::cos;
43 using std::cosh;
44 using std::exp;
45 using std::fabs;
46 using std::floor;
47 using std::fmod;
48 using std::frexp;
49 using std::ldexp;
50 using std::log10;
51 using std::log;
52 using std::modf;
53 using std::pow;
54 using std::sin;
55 using std::sinh;
56 using std::sqrt;
57 using std::tan;
58 using std::tanh;
59 #endif
60
61 #ifdef __cplusplus
62 extern "C" {
63 #endif
64
65 #if defined(__cplusplus)
66 #define exception __math_exception
67 #endif
68
69 #ifndef __P
70 #ifdef __STDC__
71 #define __P(p) p
72 #else
73 #define __P(p) ()
74 #endif
75 #endif /* !defined(__P) */
76
77 #if defined(__EXTENSIONS__) || defined(_XOPEN_SOURCE) || \
78 !defined(_STRICT_STDC) && !defined(_POSIX_C_SOURCE)
79 /*
80 * SVID & X/Open
81 */
82 #define M_E 2.7182818284590452354
83 #define M_LOG2E 1.4426950408889634074
84 #define M_LOG10E 0.43429448190325182765
85 #define M_LN2 0.69314718055994530942
86 #define M_LN10 2.30258509299404568402
87 #define M_PI 3.14159265358979323846
88 #define M_PI_2 1.57079632679489661923
89 #define M_PI_4 0.78539816339744830962
90 #define M_1_PI 0.31830988618379067154
91 #define M_2_PI 0.63661977236758134308
92 #define M_2_SQRTPI 1.12837916709551257390
93 #define M_SQRT2 1.41421356237309504880
94 #define M_SQRT1_2 0.70710678118654752440
95
96 extern int signgam;
97
98 #define MAXFLOAT ((float)3.40282346638528860e+38)
99
100 #if defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE)
101 /*
102 * SVID
103 */
104 enum version {libm_ieee = -1, c_issue_4, ansi_1, strict_ansi};
105
106 #ifdef __STDC__
107 extern const enum version _lib_version;
108 #else
109 extern enum version _lib_version;
110 #endif
111
112 struct exception {
113 int type;
114 char *name;
115 double arg1;
116 double arg2;
117 double retval;
118 };
119
120 #define HUGE MAXFLOAT
121
122 #define _ABS(x) ((x) < 0 ? -(x) : (x))
123
124 #define _REDUCE(TYPE, X, XN, C1, C2) { \
125 double x1 = (double)(TYPE)X, x2 = X - x1; \
126 X = x1 - (XN) * (C1); X += x2; X -= (XN) * (C2); }
127
128 #define DOMAIN 1
129 #define SING 2
130 #define OVERFLOW 3
131 #define UNDERFLOW 4
132 #define TLOSS 5
133 #define PLOSS 6
134
135 #define _POLY1(x, c) ((c)[0] * (x) + (c)[1])
136 #define _POLY2(x, c) (_POLY1((x), (c)) * (x) + (c)[2])
137 #define _POLY3(x, c) (_POLY2((x), (c)) * (x) + (c)[3])
138 #define _POLY4(x, c) (_POLY3((x), (c)) * (x) + (c)[4])
139 #define _POLY5(x, c) (_POLY4((x), (c)) * (x) + (c)[5])
140 #define _POLY6(x, c) (_POLY5((x), (c)) * (x) + (c)[6])
141 #define _POLY7(x, c) (_POLY6((x), (c)) * (x) + (c)[7])
142 #define _POLY8(x, c) (_POLY7((x), (c)) * (x) + (c)[8])
143 #define _POLY9(x, c) (_POLY8((x), (c)) * (x) + (c)[9])
144 #endif /* defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) */
145
146 /*
147 * SVID & X/Open
148 */
149 /* BEGIN adopted by C99 */
150 extern double erf __P((double));
151 extern double erfc __P((double));
152 extern double hypot __P((double, double));
153 extern double lgamma __P((double));
154
155 #if defined(__MATHERR_ERRNO_DONTCARE)
156 #pragma does_not_read_global_data(erf, erfc, hypot)
157 #pragma does_not_write_global_data(erf, erfc, hypot)
158 #pragma no_side_effect(erf, erfc, hypot)
159 #endif
160
161 #if !defined(_STDC_C99) && _XOPEN_SOURCE - 0 < 600 && !defined(__C99FEATURES__)
162 extern int isnan __P((double));
163
164 #pragma does_not_read_global_data(isnan)
165 #pragma does_not_write_global_data(isnan)
166 #pragma no_side_effect(isnan)
167 #endif
168 /* END adopted by C99 */
169
170 #if defined(__EXTENSIONS__) || _XOPEN_SOURCE - 0 < 600
171 extern double gamma __P((double)); /* deprecated; use lgamma */
172 #endif
173 extern double j0 __P((double));
174 extern double j1 __P((double));
175 extern double jn __P((int, double));
176 extern double y0 __P((double));
177 extern double y1 __P((double));
178 extern double yn __P((int, double));
179
180 #if defined(__MATHERR_ERRNO_DONTCARE)
181 #pragma does_not_read_global_data(j0, j1, jn, y0, y1, yn)
182 #pragma does_not_write_global_data(j0, j1, jn, y0, y1, yn)
183 #pragma no_side_effect(j0, j1, jn, y0, y1, yn)
184 #endif
185 #if defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) || \
186 _XOPEN_SOURCE - 0 >= 500 || \
187 defined(_XOPEN_SOURCE) && _XOPEN_SOURCE_EXTENDED - 0 == 1
188 /*
189 * SVID & XPG 4.2/5
190 */
191 extern double scalb __P((double, double));
192
193 #if defined(__MATHERR_ERRNO_DONTCARE)
194 #pragma does_not_read_global_data(scalb)
195 #pragma does_not_write_global_data(scalb)
196 #pragma no_side_effect(scalb)
197 #endif
198
199 /* BEGIN adopted by C99 */
200 extern double acosh __P((double));
201 extern double asinh __P((double));
202 extern double atanh __P((double));
203 extern double cbrt __P((double));
204 extern double logb __P((double));
205 extern double nextafter __P((double, double));
206 extern double remainder __P((double, double));
207
208 /*
209 * XPG 4.2/5
210 */
211 extern double expm1 __P((double));
212 extern int ilogb __P((double));
213 extern double log1p __P((double));
214 extern double rint __P((double));
215
216 #if defined(__MATHERR_ERRNO_DONTCARE)
217 #pragma does_not_read_global_data(acosh, asinh, atanh, cbrt)
218 #pragma does_not_read_global_data(logb, nextafter, remainder)
219 #pragma does_not_read_global_data(expm1, ilogb, log1p, rint)
220 #pragma does_not_write_global_data(acosh, asinh, atanh, cbrt)
221 #pragma does_not_write_global_data(logb, nextafter, remainder)
222 #pragma does_not_write_global_data(expm1, ilogb, log1p, rint)
223 #pragma no_side_effect(acosh, asinh, atanh, cbrt)
224 #pragma no_side_effect(logb, nextafter, remainder)
225 #pragma no_side_effect(expm1, ilogb, log1p, rint)
226 #endif
227 /* END adopted by C99 */
228 #endif /* defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) || ... */
229
230 #if defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE)
231 /*
232 * SVID
233 */
234 extern int matherr __P((struct exception *));
235
236 /*
237 * IEEE Test Vector
238 */
239 extern double significand __P((double));
240
241 #if defined(__MATHERR_ERRNO_DONTCARE)
242 #pragma does_not_read_global_data(significand)
243 #pragma does_not_write_global_data(significand)
244 #pragma no_side_effect(significand)
245 #endif
246
247 extern int signgamf; /* deprecated; use signgam */
248 extern int signgaml; /* deprecated; use signgam */
249
250 extern int isnanf __P((float));
251 extern int isnanl __P((long double));
252 extern float gammaf __P((float)); /* deprecated; use lgammaf */
253 extern float gammaf_r __P((float, int *)); /* deprecated; use lgammaf_r */
254 extern float j0f __P((float));
255 extern float j1f __P((float));
256 extern float jnf __P((int, float));
257 extern float lgammaf_r __P((float, int *));
258 extern float scalbf __P((float, float));
259 extern float significandf __P((float));
260 extern float y0f __P((float));
261 extern float y1f __P((float));
262 extern float ynf __P((int, float));
263 extern long double gammal __P((long double)); /* deprecated; use lgammal */
264 extern long double gammal_r __P((long double, int *)); /* deprecated */
265 extern long double j0l __P((long double));
266 extern long double j1l __P((long double));
267 extern long double jnl __P((int, long double));
268 extern long double lgammal_r __P((long double, int *));
269 extern long double scalbl __P((long double, long double));
270 extern long double significandl __P((long double));
271 extern long double y0l __P((long double));
272 extern long double y1l __P((long double));
273 extern long double ynl __P((int, long double));
274
275 #if defined(__MATHERR_ERRNO_DONTCARE)
276 #pragma does_not_read_global_data(isnanf, isnanl)
277 #pragma does_not_write_global_data(isnanf, isnanl)
278 #pragma no_side_effect(isnanf, isnanl)
279 #pragma does_not_read_global_data(gammaf_r, j0f, j1f, jnf, lgammaf_r, scalbf)
280 #pragma does_not_read_global_data(significandf, y0f, y1f, ynf)
281 #pragma does_not_write_global_data(j0f, j1f, jnf, scalbf)
282 #pragma does_not_write_global_data(significandf, y0f, y1f, ynf)
283 #pragma no_side_effect(j0f, j1f, jnf, scalbf)
284 #pragma no_side_effect(significandf, y0f, y1f, ynf)
285 #pragma does_not_read_global_data(gammal_r, j0l, j1l, jnl, lgammal_r, scalbl)
286 #pragma does_not_read_global_data(significandl, y0l, y1l, ynl)
287 #pragma does_not_write_global_data(j0l, j1l, jnl, scalbl)
288 #pragma does_not_write_global_data(significandl, y0l, y1l, ynl)
289 #pragma no_side_effect(j0l, j1l, jnl, scalbl)
290 #pragma no_side_effect(significandl, y0l, y1l, ynl)
291 #endif
292
293 /*
294 * for sin+cos->sincos transformation
295 */
296 extern void sincos __P((double, double *, double *));
297 extern void sincosf __P((float, float *, float *));
298 extern void sincosl __P((long double, long double *, long double *));
299
300 #if defined(__MATHERR_ERRNO_DONTCARE)
301 #pragma does_not_read_global_data(sincos, sincosf, sincosl)
302 #endif
303
304 /* BEGIN adopted by C99 */
305 /*
306 * Functions callable from C, intended to support IEEE arithmetic.
307 */
308 extern double copysign __P((double, double));
309 extern double scalbn __P((double, int));
310
311 #if defined(__MATHERR_ERRNO_DONTCARE)
312 #pragma does_not_read_global_data(copysign, scalbn)
313 #pragma does_not_write_global_data(copysign, scalbn)
314 #pragma no_side_effect(copysign, scalbn)
315 #endif
316 /* END adopted by C99 */
317
318 /*
319 * Reentrant version of gamma & lgamma; passes signgam back by reference
320 * as the second argument; user must allocate space for signgam.
321 */
322 extern double gamma_r __P((double, int *)); /* deprecated; use lgamma_r */
323 extern double lgamma_r __P((double, int *));
324
325 #if defined(__MATHERR_ERRNO_DONTCARE)
326 #pragma does_not_read_global_data(gamma_r, lgamma_r)
327 #endif
328
329 /* BEGIN adopted by C99 */
330 extern float modff __P((float, float *));
331
332 #if defined(__MATHERR_ERRNO_DONTCARE)
333 #pragma does_not_read_global_data(modff)
334 #endif
335 /* END adopted by C99 */
336
337 #if defined(__EXTENSIONS__) || !defined(__cplusplus)
338 #include <floatingpoint.h>
339 #endif
340 #endif /* defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) */
341 #endif /* defined(__EXTENSIONS__) || defined(_XOPEN_SOURCE) || ... */
342
343 #if defined(__cplusplus) && defined(__GNUC__)
344 #undef exception
345 #endif
346
347 #ifdef __cplusplus
348 }
349 #endif
350
351 #endif /* _MATH_H */