math.h revision e76ee993ff30a639d24c7db6a080c14d1b1a10f2
1/* 2 * ==================================================== 3 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 4 * 5 * Developed at SunPro, a Sun Microsystems, Inc. business. 6 * Permission to use, copy, modify, and distribute this 7 * software is freely granted, provided that this notice 8 * is preserved. 9 * ==================================================== 10 */ 11 12/* 13 * from: @(#)fdlibm.h 5.1 93/09/24 14 * $FreeBSD$ 15 */ 16 17#ifndef _MATH_H_ 18#define _MATH_H_ 19 20#include <sys/cdefs.h> 21#include <limits.h> 22 23__BEGIN_DECLS 24 25#define HUGE_VAL __builtin_huge_val() 26 27#define FP_ILOGB0 (-INT_MAX) 28#define FP_ILOGBNAN INT_MAX 29 30#define HUGE_VALF __builtin_huge_valf() 31#define HUGE_VALL __builtin_huge_vall() 32#define INFINITY __builtin_inff() 33#define NAN __builtin_nanf("") 34 35#define MATH_ERRNO 1 36#define MATH_ERREXCEPT 2 37#define math_errhandling MATH_ERREXCEPT 38 39#if defined(__FP_FAST_FMA) 40#define FP_FAST_FMA 1 41#endif 42#if defined(__FP_FAST_FMAF) 43#define FP_FAST_FMAF 1 44#endif 45#if defined(__FP_FAST_FMAL) 46#define FP_FAST_FMAL 1 47#endif 48 49/* Symbolic constants to classify floating point numbers. */ 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#define fpclassify(x) \ 56 __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL, FP_SUBNORMAL, FP_ZERO, x) 57 58#define isfinite(x) __builtin_isfinite(x) 59#define isinf(x) __builtin_isinf(x) 60#define isnan(x) __builtin_isnan(x) 61#define isnormal(x) __builtin_isnormal(x) 62 63#define isgreater(x, y) __builtin_isgreater((x), (y)) 64#define isgreaterequal(x, y) __builtin_isgreaterequal((x), (y)) 65#define isless(x, y) __builtin_isless((x), (y)) 66#define islessequal(x, y) __builtin_islessequal((x), (y)) 67#define islessgreater(x, y) __builtin_islessgreater((x), (y)) 68#define isunordered(x, y) __builtin_isunordered((x), (y)) 69 70#define signbit(x) \ 71 ((sizeof(x) == sizeof(float)) ? __builtin_signbitf(x) \ 72 : (sizeof(x) == sizeof(double)) ? __builtin_signbit(x) \ 73 : __builtin_signbitl(x)) 74 75typedef double __double_t; 76typedef __double_t double_t; 77typedef float __float_t; 78typedef __float_t float_t; 79 80#if defined(__USE_BSD) 81#define HUGE MAXFLOAT 82#endif 83 84extern int signgam; 85 86/* 87 * Most of these functions depend on the rounding mode and have the side 88 * effect of raising floating-point exceptions, so they are not declared 89 * as __attribute_const__. In C99, FENV_ACCESS affects the purity of these functions. 90 */ 91 92int __fpclassifyd(double) __attribute_const__; 93int __fpclassifyf(float) __attribute_const__; 94int __fpclassifyl(long double) __attribute_const__; 95int __isfinitef(float) __attribute_const__; 96int __isfinite(double) __attribute_const__; 97int __isfinitel(long double) __attribute_const__; 98int __isinff(float) __attribute_const__; 99int __isinfl(long double) __attribute_const__; 100int __isnanf(float) __attribute_const__ __INTRODUCED_IN(21); 101int __isnanl(long double) __attribute_const__; 102int __isnormalf(float) __attribute_const__; 103int __isnormal(double) __attribute_const__; 104int __isnormall(long double) __attribute_const__; 105int __signbit(double) __attribute_const__; 106int __signbitf(float) __attribute_const__; 107int __signbitl(long double) __attribute_const__; 108 109double acos(double); 110double asin(double); 111double atan(double); 112double atan2(double, double); 113double cos(double); 114double sin(double); 115double tan(double); 116 117double cosh(double); 118double sinh(double); 119double tanh(double); 120 121double exp(double); 122double frexp(double, int *); /* fundamentally !__attribute_const__ */ 123double ldexp(double, int); 124double log(double); 125double log10(double); 126double modf(double, double *); /* fundamentally !__attribute_const__ */ 127 128double pow(double, double); 129double sqrt(double); 130 131double ceil(double); 132double fabs(double) __attribute_const__; 133double floor(double); 134double fmod(double, double); 135 136double acosh(double); 137double asinh(double); 138double atanh(double); 139double cbrt(double); 140double erf(double); 141double erfc(double); 142double exp2(double); 143double expm1(double); 144double fma(double, double, double); 145double hypot(double, double); 146int ilogb(double) __attribute_const__; 147double lgamma(double); 148long long llrint(double); 149long long llround(double); 150double log1p(double); 151double log2(double) __INTRODUCED_IN(18); 152double logb(double); 153long lrint(double); 154long lround(double); 155 156double nan(const char*) __attribute_const__ __INTRODUCED_IN_ARM(13) __INTRODUCED_IN_MIPS(13) 157 __INTRODUCED_IN_X86(9); 158 159double nextafter(double, double); 160double remainder(double, double); 161double remquo(double, double, int*); 162double rint(double); 163 164double copysign(double, double) __attribute_const__; 165double fdim(double, double); 166double fmax(double, double) __attribute_const__; 167double fmin(double, double) __attribute_const__; 168double nearbyint(double); 169double round(double); 170double scalbln(double, long) __INTRODUCED_IN_X86(18) __VERSIONER_NO_GUARD; 171double scalbn(double, int); 172double tgamma(double); 173double trunc(double); 174 175float acosf(float); 176float asinf(float); 177float atanf(float); 178float atan2f(float, float); 179float cosf(float); 180float sinf(float); 181float tanf(float); 182 183float coshf(float); 184float sinhf(float); 185float tanhf(float); 186 187float exp2f(float); 188float expf(float); 189float expm1f(float); 190float frexpf(float, int *); /* fundamentally !__attribute_const__ */ 191int ilogbf(float) __attribute_const__; 192float ldexpf(float, int); 193float log10f(float); 194float log1pf(float); 195float log2f(float) __INTRODUCED_IN(18); 196float logf(float); 197float modff(float, float *); /* fundamentally !__attribute_const__ */ 198 199float powf(float, float); 200float sqrtf(float); 201 202float ceilf(float); 203float fabsf(float) __attribute_const__; 204float floorf(float); 205float fmodf(float, float); 206float roundf(float); 207 208float erff(float); 209float erfcf(float); 210float hypotf(float, float); 211float lgammaf(float); 212float tgammaf(float) __INTRODUCED_IN_ARM(13) __INTRODUCED_IN_MIPS(13) __INTRODUCED_IN_X86(9); 213 214float acoshf(float); 215float asinhf(float); 216float atanhf(float); 217float cbrtf(float); 218float logbf(float); 219float copysignf(float, float) __attribute_const__; 220long long llrintf(float); 221long long llroundf(float); 222long lrintf(float); 223long lroundf(float); 224float nanf(const char*) __attribute_const__ __INTRODUCED_IN_ARM(13) __INTRODUCED_IN_MIPS(13) 225 __INTRODUCED_IN_X86(9); 226float nearbyintf(float); 227float nextafterf(float, float); 228float remainderf(float, float); 229float remquof(float, float, int *); 230float rintf(float); 231float scalblnf(float, long) __INTRODUCED_IN_X86(18) __VERSIONER_NO_GUARD; 232float scalbnf(float, int); 233float truncf(float); 234 235float fdimf(float, float); 236float fmaf(float, float, float); 237float fmaxf(float, float) __attribute_const__; 238float fminf(float, float) __attribute_const__; 239 240long double acoshl(long double) __INTRODUCED_IN(21); 241long double acosl(long double) __INTRODUCED_IN(21); 242long double asinhl(long double) __INTRODUCED_IN(21); 243long double asinl(long double) __INTRODUCED_IN(21); 244long double atan2l(long double, long double) __INTRODUCED_IN(21); 245long double atanhl(long double) __INTRODUCED_IN(21); 246long double atanl(long double) __INTRODUCED_IN(21); 247long double cbrtl(long double) __INTRODUCED_IN(21); 248long double ceill(long double); 249long double copysignl(long double, long double) __attribute_const__; 250long double coshl(long double) __INTRODUCED_IN(21); 251long double cosl(long double) __INTRODUCED_IN(21); 252long double erfcl(long double) __INTRODUCED_IN(21); 253long double erfl(long double) __INTRODUCED_IN(21); 254long double exp2l(long double) __INTRODUCED_IN(21); 255long double expl(long double) __INTRODUCED_IN(21); 256long double expm1l(long double) __INTRODUCED_IN(21); 257long double fabsl(long double) __attribute_const__; 258long double fdiml(long double, long double); 259long double floorl(long double); 260long double fmal(long double, long double, long double) __INTRODUCED_IN(21) __VERSIONER_NO_GUARD; 261long double fmaxl(long double, long double) __attribute_const__; 262long double fminl(long double, long double) __attribute_const__; 263long double fmodl(long double, long double) __INTRODUCED_IN(21); 264long double frexpl(long double value, int*) 265 __INTRODUCED_IN(21) __VERSIONER_NO_GUARD; /* fundamentally !__attribute_const__ */ 266long double hypotl(long double, long double) __INTRODUCED_IN(21); 267int ilogbl(long double) __attribute_const__; 268long double ldexpl(long double, int); 269long double lgammal(long double) __INTRODUCED_IN(21); 270long long llrintl(long double) __INTRODUCED_IN(21); 271long long llroundl(long double); 272long double log10l(long double) __INTRODUCED_IN(21); 273long double log1pl(long double) __INTRODUCED_IN(21); 274long double log2l(long double) __INTRODUCED_IN(18); 275long double logbl(long double) __INTRODUCED_IN(18); 276long double logl(long double) __INTRODUCED_IN(21); 277long lrintl(long double) __INTRODUCED_IN(21); 278long lroundl(long double); 279long double modfl(long double, long double*) __INTRODUCED_IN(21); /* fundamentally !__attribute_const__ */ 280long double nanl(const char*) __attribute_const__ __INTRODUCED_IN(13); 281long double nearbyintl(long double) __INTRODUCED_IN(21); 282long double nextafterl(long double, long double) __INTRODUCED_IN(21) __VERSIONER_NO_GUARD; 283double nexttoward(double, long double) __INTRODUCED_IN(18) __VERSIONER_NO_GUARD; 284float nexttowardf(float, long double); 285long double nexttowardl(long double, long double) __INTRODUCED_IN(18) __VERSIONER_NO_GUARD; 286long double powl(long double, long double) __INTRODUCED_IN(21); 287long double remainderl(long double, long double) __INTRODUCED_IN(21); 288long double remquol(long double, long double, int*) __INTRODUCED_IN(21); 289long double rintl(long double) __INTRODUCED_IN(21); 290long double roundl(long double); 291long double scalblnl(long double, long) __INTRODUCED_IN_X86(18) __VERSIONER_NO_GUARD; 292long double scalbnl(long double, int); 293long double sinhl(long double) __INTRODUCED_IN(21); 294long double sinl(long double) __INTRODUCED_IN(21); 295long double sqrtl(long double) __INTRODUCED_IN(21); 296long double tanhl(long double) __INTRODUCED_IN(21); 297long double tanl(long double) __INTRODUCED_IN(21); 298long double tgammal(long double) __INTRODUCED_IN(21); 299long double truncl(long double); 300 301double j0(double); 302double j1(double); 303double jn(int, double); 304double y0(double); 305double y1(double); 306double yn(int, double); 307 308#define M_E 2.7182818284590452354 /* e */ 309#define M_LOG2E 1.4426950408889634074 /* log 2e */ 310#define M_LOG10E 0.43429448190325182765 /* log 10e */ 311#define M_LN2 0.69314718055994530942 /* log e2 */ 312#define M_LN10 2.30258509299404568402 /* log e10 */ 313#define M_PI 3.14159265358979323846 /* pi */ 314#define M_PI_2 1.57079632679489661923 /* pi/2 */ 315#define M_PI_4 0.78539816339744830962 /* pi/4 */ 316#define M_1_PI 0.31830988618379067154 /* 1/pi */ 317#define M_2_PI 0.63661977236758134308 /* 2/pi */ 318#define M_2_SQRTPI 1.12837916709551257390 /* 2/sqrt(pi) */ 319#define M_SQRT2 1.41421356237309504880 /* sqrt(2) */ 320#define M_SQRT1_2 0.70710678118654752440 /* 1/sqrt(2) */ 321 322#define MAXFLOAT ((float)3.40282346638528860e+38) 323 324#if defined(__USE_BSD) || defined(__USE_GNU) 325double gamma(double); 326double scalb(double, double); 327double drem(double, double); 328int finite(double) __attribute_const__; 329int isnanf(float) __attribute_const__; 330double gamma_r(double, int*); 331double lgamma_r(double, int*); 332double significand(double); 333long double lgammal_r(long double, int*) __INTRODUCED_IN(23); 334long double significandl(long double) __INTRODUCED_IN(21); 335float dremf(float, float); 336int finitef(float) __attribute_const__; 337float gammaf(float); 338float j0f(float); 339float j1f(float); 340float jnf(int, float); 341float scalbf(float, float); 342float y0f(float); 343float y1f(float); 344float ynf(int, float); 345float gammaf_r(float, int *); 346float lgammaf_r(float, int *); 347float significandf(float); 348#endif 349 350#if defined(__USE_GNU) 351#define M_El 2.718281828459045235360287471352662498L /* e */ 352#define M_LOG2El 1.442695040888963407359924681001892137L /* log 2e */ 353#define M_LOG10El 0.434294481903251827651128918916605082L /* log 10e */ 354#define M_LN2l 0.693147180559945309417232121458176568L /* log e2 */ 355#define M_LN10l 2.302585092994045684017991454684364208L /* log e10 */ 356#define M_PIl 3.141592653589793238462643383279502884L /* pi */ 357#define M_PI_2l 1.570796326794896619231321691639751442L /* pi/2 */ 358#define M_PI_4l 0.785398163397448309615660845819875721L /* pi/4 */ 359#define M_1_PIl 0.318309886183790671537767526745028724L /* 1/pi */ 360#define M_2_PIl 0.636619772367581343075535053490057448L /* 2/pi */ 361#define M_2_SQRTPIl 1.128379167095512573896158903121545172L /* 2/sqrt(pi) */ 362#define M_SQRT2l 1.414213562373095048801688724209698079L /* sqrt(2) */ 363#define M_SQRT1_2l 0.707106781186547524400844362104849039L /* 1/sqrt(2) */ 364void sincos(double, double*, double*); 365void sincosf(float, float*, float*); 366void sincosl(long double, long double*, long double*); 367#endif 368 369__END_DECLS 370 371#endif /* !_MATH_H_ */ 372