Annotation of src/include/math.h, Revision 1.20
1.20 ! martynas 1: /* $OpenBSD: math.h,v 1.19 2008/07/24 09:35:30 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: extern double exp2(double);
215: extern double expm1(double);
216: extern int ilogb(double);
217: extern double log1p(double);
218: #if 0
219: extern double log2(double);
220: #endif
221: extern double logb(double);
222: extern double scalbn(double, int);
223: #if 0
224: extern double scalbln(double, long int);
225: #endif
226:
227: extern double cbrt(double);
228: extern double hypot(double, double);
229:
230: extern double erf(double);
231: extern double erfc(double);
232: extern double lgamma(double);
233: extern double tgamma(double);
234:
235: #if 0
236: extern double nearbyint(double);
237: #endif
238: extern double rint(double);
239: extern long int lrint(double);
240: extern long long int llrint(double);
1.14 brad 241: extern double round(double);
1.18 martynas 242: extern long int lround(double);
243: extern long long int llround(double);
1.13 brad 244: extern double trunc(double);
245:
1.18 martynas 246: extern double remainder(double, double);
247: extern double remquo(double, double, int *);
248:
249: extern double copysign(double, double);
250: extern double nan(const char *);
251: extern double nextafter(double, double);
252: #if 0
253: extern double nexttoward(double, long double);
254: #endif
255:
256: #if 0
257: extern double fdim(double, double);
258: extern double fmax(double, double);
259: extern double fmin(double, double);
260: #endif
261:
262: #if 0
263: extern double fma(double, double, double);
264: #endif
265: #endif /* __BSD_VISIBLE || __ISO_C_VISIBLE >= 1999 || __XPG_VISIBLE */
266:
1.11 millert 267: #if __BSD_VISIBLE || __XPG_VISIBLE
1.6 millert 268: extern double j0(double);
269: extern double j1(double);
270: extern double jn(int, double);
1.18 martynas 271: extern double scalb(double, double);
1.6 millert 272: extern double y0(double);
273: extern double y1(double);
274: extern double yn(int, double);
1.18 martynas 275: #endif /* __BSD_VISIBLE || __XPG_VISIBLE */
1.1 deraadt 276:
1.18 martynas 277: #if __BSD_VISIBLE || __XPG_VISIBLE <= 500
278: extern double gamma(double);
279: #endif /* __BSD_VISIBLE || __XPG_VISIBLE <= 500 */
1.1 deraadt 280:
281: /*
282: * BSD math library entry points
283: */
1.18 martynas 284: #if __BSD_VISIBLE
1.6 millert 285: extern double drem(double, double);
1.18 martynas 286: extern int finite(double);
1.1 deraadt 287:
288: /*
289: * Reentrant version of gamma & lgamma; passes signgam back by reference
290: * as the second argument; user must allocate space for signgam.
291: */
1.6 millert 292: extern double gamma_r(double, int *);
293: extern double lgamma_r(double, int *);
1.1 deraadt 294:
1.18 martynas 295: #ifdef __LIBM_PRIVATE
296: extern int matherr(struct exception *);
297: #endif /* __LIBM_PRIVATE */
298:
299: /*
300: * IEEE Test Vector
301: */
302: extern double significand(double);
303: #endif /* __BSD_VISIBLE */
1.1 deraadt 304:
1.18 martynas 305: /*
306: * Float versions of C99 functions
307: */
308: #if __ISO_C_VISIBLE >= 1999
1.6 millert 309: extern float acosf(float);
310: extern float asinf(float);
311: extern float atanf(float);
312: extern float atan2f(float, float);
313: extern float cosf(float);
314: extern float sinf(float);
315: extern float tanf(float);
316:
1.18 martynas 317: extern float acoshf(float);
318: extern float asinhf(float);
319: extern float atanhf(float);
1.6 millert 320: extern float coshf(float);
321: extern float sinhf(float);
322: extern float tanhf(float);
323:
324: extern float expf(float);
1.18 martynas 325: extern float exp2f(float);
326: extern float expm1f(float);
1.6 millert 327: extern float frexpf(float, int *);
1.18 martynas 328: extern int ilogbf(float);
1.6 millert 329: extern float ldexpf(float, int);
330: extern float logf(float);
331: extern float log10f(float);
1.18 martynas 332: extern float log1pf(float);
333: #if 0
334: extern float log2f(float);
335: #endif
336: extern float logbf(float);
1.6 millert 337: extern float modff(float, float *);
1.18 martynas 338: extern float scalbnf(float, int);
339: #if 0
340: extern float scalblnf(float, long int);
341: #endif
1.6 millert 342:
1.18 martynas 343: extern float cbrtf(float);
344: extern float fabsf(float);
345: extern float hypotf(float, float);
1.6 millert 346: extern float powf(float, float);
347: extern float sqrtf(float);
348:
1.18 martynas 349: extern float erff(float);
350: extern float erfcf(float);
351: extern float lgammaf(float);
352: extern float tgammaf(float);
353:
1.6 millert 354: extern float ceilf(float);
355: extern float floorf(float);
1.18 martynas 356: #if 0
357: extern float nearbyintf(float);
358: #endif
359: extern float rintf(float);
360: extern long int lrintf(float);
361: extern long long int llrintf(float);
362: extern float roundf(float);
363: extern long int lroundf(float);
364: extern long long int llroundf(float);
365: extern float truncf(float);
366:
1.6 millert 367: extern float fmodf(float, float);
1.18 martynas 368: extern float remainderf(float, float);
369: extern float remquof(float, float, int *);
370:
371: extern float copysignf(float, float);
372: extern float nanf(const char *);
373: extern float nextafterf(float, float);
374: #if 0
375: extern float nexttowardf(float, long double);
376: #endif
1.13 brad 377:
1.18 martynas 378: #if 0
379: extern float fdimf(float, float);
380: extern float fmaxf(float, float);
381: extern float fminf(float, float);
382: #endif
383:
384: #if 0
385: extern float fmaf(float, float, float);
386: #endif
387: #endif /* __ISO_C_VISIBLE >= 1999 */
1.6 millert 388:
1.18 martynas 389: #if __BSD_VISIBLE || __XPG_VISIBLE
1.6 millert 390: extern float j0f(float);
391: extern float j1f(float);
392: extern float jnf(int, float);
1.18 martynas 393: extern float scalbf(float, float);
1.6 millert 394: extern float y0f(float);
395: extern float y1f(float);
396: extern float ynf(int, float);
1.18 martynas 397: #endif /* __BSD_VISIBLE || __XPG_VISIBLE */
1.6 millert 398:
1.18 martynas 399: #if __BSD_VISIBLE || __XPG_VISIBLE <= 500
400: extern float gammaf(float);
401: #endif /* __BSD_VISIBLE || __XPG_VISIBLE <= 500 */
1.1 deraadt 402:
403: /*
1.18 martynas 404: * Float versions of BSD math library entry points
1.1 deraadt 405: */
1.18 martynas 406: #if __BSD_VISIBLE
407: extern float dremf(float, float);
408: extern int finitef(float);
409: extern int isinff(float);
410: extern int isnanf(float);
1.1 deraadt 411:
412: /*
1.18 martynas 413: * Float versions of reentrant version of gamma & lgamma; passes
414: * signgam back by reference as the second argument; user must
415: * allocate space for signgam.
1.1 deraadt 416: */
1.18 martynas 417: extern float gammaf_r(float, int *);
418: extern float lgammaf_r(float, int *);
1.1 deraadt 419:
420: /*
1.18 martynas 421: * Float version of IEEE Test Vector
1.1 deraadt 422: */
1.18 martynas 423: extern float significandf(float);
424: #endif /* __BSD_VISIBLE */
1.1 deraadt 425:
426: /*
1.18 martynas 427: * Long double versions of C99 functions
1.1 deraadt 428: */
1.18 martynas 429: #if __ISO_C_VISIBLE >= 1999
430: #if 0
431: extern long double acosl(long double);
432: extern long double asinl(long double);
433: extern long double atanl(long double);
434: extern long double atan2l(long double, long double);
435: extern long double cosl(long double);
436: extern long double sinl(long double);
437: extern long double tanl(long double);
438:
439: extern long double acoshl(long double);
440: extern long double asinhl(long double);
441: extern long double atanhl(long double);
442: extern long double coshl(long double);
443: extern long double sinhl(long double);
444: extern long double tanhl(long double);
445:
446: extern long double expl(long double);
447: extern long double exp2l(long double);
448: extern long double expm1l(long double);
449: extern long double frexpl(long double, int *);
450: extern int ilogbl(long double);
451: extern long double ldexpl(long double, int);
452: extern long double logl(long double);
453: extern long double log10l(long double);
454: extern long double log1pl(long double);
455: extern long double log2l(long double);
456: extern long double logbl(long double);
457: extern long double modfl(long double, long double *);
458: extern long double scalbnl(long double, int);
459: extern long double scalblnl(long double, long int);
460:
461: extern long double cbrtl(long double);
462: extern long double fabsl(long double);
463: extern long double hypotl(long double, long double);
464: extern long double powl(long double, long double);
465: extern long double sqrtl(long double);
466:
467: extern long double erfl(long double);
468: extern long double erfcl(long double);
469: extern long double lgammal(long double);
470: extern long double tgammal(long double);
471:
472: extern long double ceill(long double);
473: extern long double floorl(long double);
474: extern long double nearbyintl(long double);
475: extern long double rintl(long double);
476: extern long int lrintl(long double);
477: extern long long int llrintl(long double);
478: extern long double roundl(long double);
479: extern long int lroundl(long double);
480: extern long long int llroundl(long double);
481: extern long double truncl(long double);
482:
483: extern long double fmodl(long double, long double);
484: extern long double remainderl(long double, long double);
485: extern long double remquol(long double, long double, int *);
486:
487: extern long double copysignl(long double, long double);
488: extern long double nanl(const char *);
489: extern long double nextafterl(long double, long double);
490: extern long double nexttowardl(long double, long double);
491:
492: extern long double fdiml(long double, long double);
493: extern long double fmaxl(long double, long double);
494: extern long double fminl(long double, long double);
495:
496: extern long double fmal(long double, long double, long double);
497: #endif
498: #endif /* __ISO_C_VISIBLE >= 1999 */
1.1 deraadt 499:
1.18 martynas 500: /*
501: * Library implementation
502: */
503: extern int __fpclassify(double);
504: extern int __fpclassifyf(float);
505: extern int __fpclassifyl(long double);
506: extern int __isfinite(double);
507: extern int __isfinitef(float);
508: extern int __isfinitel(long double);
509: extern int __isinf(double);
510: extern int __isinfl(long double);
511: extern int __isnan(double);
512: extern int __isnanl(long double);
513: extern int __isnormal(double);
514: extern int __isnormalf(float);
515: extern int __isnormall(long double);
516: extern int __signbit(double);
517: extern int __signbitf(float);
518: extern int __signbitl(long double);
1.16 martynas 519:
520: #if __BSD_VISIBLE && defined(__vax__)
521: extern double infnan(int);
522: #endif /* __BSD_VISIBLE && defined(__vax__) */
1.1 deraadt 523: __END_DECLS
524:
1.18 martynas 525: #endif /* !_MATH_H_ */