Annotation of src/include/math.h, Revision 1.19
1.19 ! martynas 1: /* $OpenBSD: math.h,v 1.18 2008/07/22 16:01:46 martynas Exp $ */
1.1 deraadt 2: /*
3: * ====================================================
4: * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
5: *
6: * Developed at SunPro, a Sun Microsystems, Inc. business.
7: * Permission to use, copy, modify, and distribute this
8: * software is freely granted, provided that this notice
9: * is preserved.
10: * ====================================================
11: */
12:
13: /*
1.3 millert 14: * from: @(#)fdlibm.h 5.1 93/09/24
1.1 deraadt 15: */
16:
17: #ifndef _MATH_H_
18: #define _MATH_H_
19:
1.17 martynas 20: #include <sys/_types.h>
1.11 millert 21: #include <sys/cdefs.h>
1.18 martynas 22: #include <sys/limits.h>
1.11 millert 23:
1.1 deraadt 24: /*
25: * ANSI/POSIX
26: */
27: extern char __infinity[];
1.12 otto 28: #define HUGE_VAL (*(double *)(void *)__infinity)
1.9 espie 29:
1.18 martynas 30: /*
1.9 espie 31: * C99
32: */
1.18 martynas 33: #if __ISO_C_VISIBLE >= 1999
34: typedef __double_t double_t;
35: typedef __float_t float_t;
36:
37: #define HUGE_VALF ((float)HUGE_VAL)
38: #define HUGE_VALL ((long double)HUGE_VAL)
39: #define INFINITY HUGE_VALF
1.19 ! martynas 40: #ifndef __vax__
! 41: extern char __nan[];
! 42: #define NAN (*(float *)(void *)__nan)
! 43: #endif /* !__vax__ */
1.18 martynas 44:
45: #define FP_INFINITE 0x01
46: #define FP_NAN 0x02
47: #define FP_NORMAL 0x04
48: #define FP_SUBNORMAL 0x08
49: #define FP_ZERO 0x10
50:
51: #define FP_ILOGB0 (-INT_MAX)
52: #define FP_ILOGBNAN INT_MAX
53:
54: #define fpclassify(x) \
55: ((sizeof (x) == sizeof (float)) ? \
56: __fpclassifyf(x) \
57: : (sizeof (x) == sizeof (double)) ? \
58: __fpclassify(x) \
59: : __fpclassifyl(x))
60: #define isfinite(x) \
61: ((sizeof (x) == sizeof (float)) ? \
62: __isfinitef(x) \
63: : (sizeof (x) == sizeof (double)) ? \
64: __isfinite(x) \
65: : __isfinitel(x))
66: #define isinf(x) \
67: ((sizeof (x) == sizeof (float)) ? \
68: isinff(x) \
69: : (sizeof (x) == sizeof (double)) ? \
70: __isinf(x) \
71: : __isinfl(x))
72: #define isnan(x) \
73: ((sizeof (x) == sizeof (float)) ? \
74: isnanf(x) \
75: : (sizeof (x) == sizeof (double)) ? \
76: __isnan(x) \
77: : __isnanl(x))
78: #define isnormal(x) \
79: ((sizeof (x) == sizeof (float)) ? \
80: __isnormalf(x) \
81: : (sizeof (x) == sizeof (double)) ? \
82: __isnormal(x) \
83: : __isnormall(x))
84: #define signbit(x) \
85: ((sizeof (x) == sizeof (float)) ? \
86: __signbitf(x) \
87: : (sizeof (x) == sizeof (double)) ? \
88: __signbit(x) \
89: : __signbitl(x))
90:
91: #define isgreater(x, y) (!isunordered((x), (y)) && (x) > (y))
92: #define isgreaterequal(x, y) (!isunordered((x), (y)) && (x) >= (y))
93: #define isless(x, y) (!isunordered((x), (y)) && (x) < (y))
94: #define islessequal(x, y) (!isunordered((x), (y)) && (x) <= (y))
95: #define islessgreater(x, y) (!isunordered((x), (y)) && \
96: ((x) > (y) || (y) > (x)))
97: #define isunordered(x, y) (isnan(x) || isnan(y))
98: #endif /* __ISO_C_VISIBLE >= 1999 */
1.1 deraadt 99:
100: /*
101: * XOPEN/SVID
102: */
1.11 millert 103: #if __BSD_VISIBLE || __XPG_VISIBLE
1.1 deraadt 104: #define M_E 2.7182818284590452354 /* e */
105: #define M_LOG2E 1.4426950408889634074 /* log 2e */
106: #define M_LOG10E 0.43429448190325182765 /* log 10e */
107: #define M_LN2 0.69314718055994530942 /* log e2 */
108: #define M_LN10 2.30258509299404568402 /* log e10 */
109: #define M_PI 3.14159265358979323846 /* pi */
110: #define M_PI_2 1.57079632679489661923 /* pi/2 */
111: #define M_PI_4 0.78539816339744830962 /* pi/4 */
112: #define M_1_PI 0.31830988618379067154 /* 1/pi */
113: #define M_2_PI 0.63661977236758134308 /* 2/pi */
114: #define M_2_SQRTPI 1.12837916709551257390 /* 2/sqrt(pi) */
115: #define M_SQRT2 1.41421356237309504880 /* sqrt(2) */
116: #define M_SQRT1_2 0.70710678118654752440 /* 1/sqrt(2) */
117:
1.7 jason 118: #ifdef __vax__
1.18 martynas 119: #define MAXFLOAT ((float)1.70141173319264430e+38)
1.7 jason 120: #else
1.1 deraadt 121: #define MAXFLOAT ((float)3.40282346638528860e+38)
1.18 martynas 122: #endif /* __vax__ */
1.7 jason 123:
1.1 deraadt 124: extern int signgam;
1.11 millert 125: #endif /* __BSD_VISIBLE || __XPG_VISIBLE */
1.1 deraadt 126:
1.11 millert 127: #if __BSD_VISIBLE
1.1 deraadt 128: enum fdversion {fdlibm_ieee = -1, fdlibm_svid, fdlibm_xopen, fdlibm_posix};
129:
130: #define _LIB_VERSION_TYPE enum fdversion
1.18 martynas 131: #define _LIB_VERSION _fdlib_version
1.1 deraadt 132:
1.18 martynas 133: /* if global variable _LIB_VERSION is not desirable, one may
134: * change the following to be a constant by:
1.1 deraadt 135: * #define _LIB_VERSION_TYPE const enum version
136: * In that case, after one initializes the value _LIB_VERSION (see
137: * s_lib_version.c) during compile time, it cannot be modified
138: * in the middle of a program
1.18 martynas 139: */
140: extern _LIB_VERSION_TYPE _LIB_VERSION;
1.1 deraadt 141:
1.18 martynas 142: #define _IEEE_ fdlibm_ieee
143: #define _SVID_ fdlibm_svid
144: #define _XOPEN_ fdlibm_xopen
145: #define _POSIX_ fdlibm_posix
1.1 deraadt 146:
1.5 millert 147: #ifndef __cplusplus
1.1 deraadt 148: struct exception {
149: int type;
150: char *name;
151: double arg1;
152: double arg2;
153: double retval;
154: };
1.18 martynas 155: #endif /* !__cplusplus */
1.1 deraadt 156:
157: #define HUGE MAXFLOAT
158:
1.18 martynas 159: /*
1.1 deraadt 160: * set X_TLOSS = pi*2**52, which is possibly defined in <values.h>
161: * (one may replace the following line by "#include <values.h>")
162: */
163:
1.18 martynas 164: #define X_TLOSS 1.41484755040568800000e+16
1.1 deraadt 165:
166: #define DOMAIN 1
167: #define SING 2
168: #define OVERFLOW 3
169: #define UNDERFLOW 4
170: #define TLOSS 5
171: #define PLOSS 6
172:
1.11 millert 173: #endif /* __BSD_VISIBLE */
1.1 deraadt 174:
175: __BEGIN_DECLS
176: /*
177: * ANSI/POSIX
178: */
1.6 millert 179: extern double acos(double);
180: extern double asin(double);
181: extern double atan(double);
182: extern double atan2(double, double);
183: extern double cos(double);
184: extern double sin(double);
185: extern double tan(double);
186:
187: extern double cosh(double);
188: extern double sinh(double);
189: extern double tanh(double);
190:
191: extern double exp(double);
192: extern double frexp(double, int *);
193: extern double ldexp(double, int);
194: extern double log(double);
195: extern double log10(double);
196: extern double modf(double, double *);
197:
198: extern double pow(double, double);
199: extern double sqrt(double);
200:
201: extern double ceil(double);
202: extern double fabs(double);
203: extern double floor(double);
204: extern double fmod(double, double);
1.1 deraadt 205:
1.18 martynas 206: /*
207: * C99
208: */
209: #if __BSD_VISIBLE || __ISO_C_VISIBLE >= 1999 || __XPG_VISIBLE
210: extern double acosh(double);
211: extern double asinh(double);
212: extern double atanh(double);
213:
214: #if 0 /* XXX */
215: extern double exp2(double);
216: #endif /* XXX */
217: extern double expm1(double);
218: extern int ilogb(double);
219: extern double log1p(double);
220: #if 0
221: extern double log2(double);
222: #endif
223: extern double logb(double);
224: extern double scalbn(double, int);
225: #if 0
226: extern double scalbln(double, long int);
227: #endif
228:
229: extern double cbrt(double);
230: extern double hypot(double, double);
231:
232: extern double erf(double);
233: extern double erfc(double);
234: extern double lgamma(double);
235: extern double tgamma(double);
236:
237: #if 0
238: extern double nearbyint(double);
239: #endif
240: extern double rint(double);
241: extern long int lrint(double);
242: extern long long int llrint(double);
1.14 brad 243: extern double round(double);
1.18 martynas 244: extern long int lround(double);
245: extern long long int llround(double);
1.13 brad 246: extern double trunc(double);
247:
1.18 martynas 248: extern double remainder(double, double);
249: #if 0 /* XXX */
250: extern double remquo(double, double, int *);
251: #endif /* XXX */
252:
253: extern double copysign(double, double);
254: #if 0 /* XXX */
255: extern double nan(const char *);
256: #endif /* XXX */
257: extern double nextafter(double, double);
258: #if 0
259: extern double nexttoward(double, long double);
260: #endif
261:
262: #if 0
263: extern double fdim(double, double);
264: extern double fmax(double, double);
265: extern double fmin(double, double);
266: #endif
267:
268: #if 0
269: extern double fma(double, double, double);
270: #endif
271: #endif /* __BSD_VISIBLE || __ISO_C_VISIBLE >= 1999 || __XPG_VISIBLE */
272:
1.11 millert 273: #if __BSD_VISIBLE || __XPG_VISIBLE
1.6 millert 274: extern double j0(double);
275: extern double j1(double);
276: extern double jn(int, double);
1.18 martynas 277: extern double scalb(double, double);
1.6 millert 278: extern double y0(double);
279: extern double y1(double);
280: extern double yn(int, double);
1.18 martynas 281: #endif /* __BSD_VISIBLE || __XPG_VISIBLE */
1.1 deraadt 282:
1.18 martynas 283: #if __BSD_VISIBLE || __XPG_VISIBLE <= 500
284: extern double gamma(double);
285: #endif /* __BSD_VISIBLE || __XPG_VISIBLE <= 500 */
1.1 deraadt 286:
287: /*
288: * BSD math library entry points
289: */
1.18 martynas 290: #if __BSD_VISIBLE
1.6 millert 291: extern double drem(double, double);
1.18 martynas 292: extern int finite(double);
1.1 deraadt 293:
294: /*
295: * Reentrant version of gamma & lgamma; passes signgam back by reference
296: * as the second argument; user must allocate space for signgam.
297: */
1.6 millert 298: extern double gamma_r(double, int *);
299: extern double lgamma_r(double, int *);
1.1 deraadt 300:
1.18 martynas 301: #ifdef __LIBM_PRIVATE
302: extern int matherr(struct exception *);
303: #endif /* __LIBM_PRIVATE */
304:
305: /*
306: * IEEE Test Vector
307: */
308: extern double significand(double);
309: #endif /* __BSD_VISIBLE */
1.1 deraadt 310:
1.18 martynas 311: /*
312: * Float versions of C99 functions
313: */
314: #if __ISO_C_VISIBLE >= 1999
1.6 millert 315: extern float acosf(float);
316: extern float asinf(float);
317: extern float atanf(float);
318: extern float atan2f(float, float);
319: extern float cosf(float);
320: extern float sinf(float);
321: extern float tanf(float);
322:
1.18 martynas 323: extern float acoshf(float);
324: extern float asinhf(float);
325: extern float atanhf(float);
1.6 millert 326: extern float coshf(float);
327: extern float sinhf(float);
328: extern float tanhf(float);
329:
330: extern float expf(float);
1.18 martynas 331: #if 0 /* XXX */
332: extern float exp2f(float);
333: #endif /* XXX */
334: extern float expm1f(float);
1.6 millert 335: extern float frexpf(float, int *);
1.18 martynas 336: extern int ilogbf(float);
1.6 millert 337: extern float ldexpf(float, int);
338: extern float logf(float);
339: extern float log10f(float);
1.18 martynas 340: extern float log1pf(float);
341: #if 0
342: extern float log2f(float);
343: #endif
344: extern float logbf(float);
1.6 millert 345: extern float modff(float, float *);
1.18 martynas 346: extern float scalbnf(float, int);
347: #if 0
348: extern float scalblnf(float, long int);
349: #endif
1.6 millert 350:
1.18 martynas 351: extern float cbrtf(float);
352: extern float fabsf(float);
353: extern float hypotf(float, float);
1.6 millert 354: extern float powf(float, float);
355: extern float sqrtf(float);
356:
1.18 martynas 357: extern float erff(float);
358: extern float erfcf(float);
359: extern float lgammaf(float);
360: extern float tgammaf(float);
361:
1.6 millert 362: extern float ceilf(float);
363: extern float floorf(float);
1.18 martynas 364: #if 0
365: extern float nearbyintf(float);
366: #endif
367: extern float rintf(float);
368: extern long int lrintf(float);
369: extern long long int llrintf(float);
370: extern float roundf(float);
371: extern long int lroundf(float);
372: extern long long int llroundf(float);
373: extern float truncf(float);
374:
1.6 millert 375: extern float fmodf(float, float);
1.18 martynas 376: extern float remainderf(float, float);
377: #if 0 /* XXX */
378: extern float remquof(float, float, int *);
379: #endif /* XXX */
380:
381: extern float copysignf(float, float);
382: #if 0 /* XXX */
383: extern float nanf(const char *);
384: #endif /* XXX */
385: extern float nextafterf(float, float);
386: #if 0
387: extern float nexttowardf(float, long double);
388: #endif
1.13 brad 389:
1.18 martynas 390: #if 0
391: extern float fdimf(float, float);
392: extern float fmaxf(float, float);
393: extern float fminf(float, float);
394: #endif
395:
396: #if 0
397: extern float fmaf(float, float, float);
398: #endif
399: #endif /* __ISO_C_VISIBLE >= 1999 */
1.6 millert 400:
1.18 martynas 401: #if __BSD_VISIBLE || __XPG_VISIBLE
1.6 millert 402: extern float j0f(float);
403: extern float j1f(float);
404: extern float jnf(int, float);
1.18 martynas 405: extern float scalbf(float, float);
1.6 millert 406: extern float y0f(float);
407: extern float y1f(float);
408: extern float ynf(int, float);
1.18 martynas 409: #endif /* __BSD_VISIBLE || __XPG_VISIBLE */
1.6 millert 410:
1.18 martynas 411: #if __BSD_VISIBLE || __XPG_VISIBLE <= 500
412: extern float gammaf(float);
413: #endif /* __BSD_VISIBLE || __XPG_VISIBLE <= 500 */
1.1 deraadt 414:
415: /*
1.18 martynas 416: * Float versions of BSD math library entry points
1.1 deraadt 417: */
1.18 martynas 418: #if __BSD_VISIBLE
419: extern float dremf(float, float);
420: extern int finitef(float);
421: extern int isinff(float);
422: extern int isnanf(float);
1.1 deraadt 423:
424: /*
1.18 martynas 425: * Float versions of reentrant version of gamma & lgamma; passes
426: * signgam back by reference as the second argument; user must
427: * allocate space for signgam.
1.1 deraadt 428: */
1.18 martynas 429: extern float gammaf_r(float, int *);
430: extern float lgammaf_r(float, int *);
1.1 deraadt 431:
432: /*
1.18 martynas 433: * Float version of IEEE Test Vector
1.1 deraadt 434: */
1.18 martynas 435: extern float significandf(float);
436: #endif /* __BSD_VISIBLE */
1.1 deraadt 437:
438: /*
1.18 martynas 439: * Long double versions of C99 functions
1.1 deraadt 440: */
1.18 martynas 441: #if __ISO_C_VISIBLE >= 1999
442: #if 0
443: extern long double acosl(long double);
444: extern long double asinl(long double);
445: extern long double atanl(long double);
446: extern long double atan2l(long double, long double);
447: extern long double cosl(long double);
448: extern long double sinl(long double);
449: extern long double tanl(long double);
450:
451: extern long double acoshl(long double);
452: extern long double asinhl(long double);
453: extern long double atanhl(long double);
454: extern long double coshl(long double);
455: extern long double sinhl(long double);
456: extern long double tanhl(long double);
457:
458: extern long double expl(long double);
459: extern long double exp2l(long double);
460: extern long double expm1l(long double);
461: extern long double frexpl(long double, int *);
462: extern int ilogbl(long double);
463: extern long double ldexpl(long double, int);
464: extern long double logl(long double);
465: extern long double log10l(long double);
466: extern long double log1pl(long double);
467: extern long double log2l(long double);
468: extern long double logbl(long double);
469: extern long double modfl(long double, long double *);
470: extern long double scalbnl(long double, int);
471: extern long double scalblnl(long double, long int);
472:
473: extern long double cbrtl(long double);
474: extern long double fabsl(long double);
475: extern long double hypotl(long double, long double);
476: extern long double powl(long double, long double);
477: extern long double sqrtl(long double);
478:
479: extern long double erfl(long double);
480: extern long double erfcl(long double);
481: extern long double lgammal(long double);
482: extern long double tgammal(long double);
483:
484: extern long double ceill(long double);
485: extern long double floorl(long double);
486: extern long double nearbyintl(long double);
487: extern long double rintl(long double);
488: extern long int lrintl(long double);
489: extern long long int llrintl(long double);
490: extern long double roundl(long double);
491: extern long int lroundl(long double);
492: extern long long int llroundl(long double);
493: extern long double truncl(long double);
494:
495: extern long double fmodl(long double, long double);
496: extern long double remainderl(long double, long double);
497: extern long double remquol(long double, long double, int *);
498:
499: extern long double copysignl(long double, long double);
500: extern long double nanl(const char *);
501: extern long double nextafterl(long double, long double);
502: extern long double nexttowardl(long double, long double);
503:
504: extern long double fdiml(long double, long double);
505: extern long double fmaxl(long double, long double);
506: extern long double fminl(long double, long double);
507:
508: extern long double fmal(long double, long double, long double);
509: #endif
510: #endif /* __ISO_C_VISIBLE >= 1999 */
1.1 deraadt 511:
1.18 martynas 512: /*
513: * Library implementation
514: */
515: extern int __fpclassify(double);
516: extern int __fpclassifyf(float);
517: extern int __fpclassifyl(long double);
518: extern int __isfinite(double);
519: extern int __isfinitef(float);
520: extern int __isfinitel(long double);
521: extern int __isinf(double);
522: extern int __isinfl(long double);
523: extern int __isnan(double);
524: extern int __isnanl(long double);
525: extern int __isnormal(double);
526: extern int __isnormalf(float);
527: extern int __isnormall(long double);
528: extern int __signbit(double);
529: extern int __signbitf(float);
530: extern int __signbitl(long double);
1.16 martynas 531:
532: #if __BSD_VISIBLE && defined(__vax__)
533: extern double infnan(int);
534: #endif /* __BSD_VISIBLE && defined(__vax__) */
1.1 deraadt 535: __END_DECLS
536:
1.18 martynas 537: #endif /* !_MATH_H_ */