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