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 #if defined(__EXTENSIONS__) || defined(_XOPEN_SOURCE) || \ 70 !defined(_STRICT_STDC) && !defined(_POSIX_C_SOURCE) 71 /* 72 * SVID & X/Open 73 */ 74 #define M_E 2.7182818284590452354 75 #define M_LOG2E 1.4426950408889634074 76 #define M_LOG10E 0.43429448190325182765 77 #define M_LN2 0.69314718055994530942 78 #define M_LN10 2.30258509299404568402 79 #define M_PI 3.14159265358979323846 80 #define M_PI_2 1.57079632679489661923 81 #define M_PI_4 0.78539816339744830962 82 #define M_1_PI 0.31830988618379067154 83 #define M_2_PI 0.63661977236758134308 84 #define M_2_SQRTPI 1.12837916709551257390 85 #define M_SQRT2 1.41421356237309504880 86 #define M_SQRT1_2 0.70710678118654752440 87 88 extern int signgam; 89 90 #define MAXFLOAT ((float)3.40282346638528860e+38) 91 92 #if defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) 93 /* 94 * SVID 95 */ 96 enum version {libm_ieee = -1, c_issue_4, ansi_1, strict_ansi}; 97 98 #ifdef __STDC__ 99 extern const enum version _lib_version; 100 #else 101 extern enum version _lib_version; 102 #endif 103 104 struct exception { 105 int type; 106 char *name; 107 double arg1; 108 double arg2; 109 double retval; 110 }; 111 112 #define HUGE MAXFLOAT 113 114 #define _ABS(x) ((x) < 0 ? -(x) : (x)) 115 116 #define _REDUCE(TYPE, X, XN, C1, C2) { \ 117 double x1 = (double)(TYPE)X, x2 = X - x1; \ 118 X = x1 - (XN) * (C1); X += x2; X -= (XN) * (C2); } 119 120 #define DOMAIN 1 121 #define SING 2 122 #define OVERFLOW 3 123 #define UNDERFLOW 4 124 #define TLOSS 5 125 #define PLOSS 6 126 127 #define _POLY1(x, c) ((c)[0] * (x) + (c)[1]) 128 #define _POLY2(x, c) (_POLY1((x), (c)) * (x) + (c)[2]) 129 #define _POLY3(x, c) (_POLY2((x), (c)) * (x) + (c)[3]) 130 #define _POLY4(x, c) (_POLY3((x), (c)) * (x) + (c)[4]) 131 #define _POLY5(x, c) (_POLY4((x), (c)) * (x) + (c)[5]) 132 #define _POLY6(x, c) (_POLY5((x), (c)) * (x) + (c)[6]) 133 #define _POLY7(x, c) (_POLY6((x), (c)) * (x) + (c)[7]) 134 #define _POLY8(x, c) (_POLY7((x), (c)) * (x) + (c)[8]) 135 #define _POLY9(x, c) (_POLY8((x), (c)) * (x) + (c)[9]) 136 #endif /* defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) */ 137 138 /* 139 * SVID & X/Open 140 */ 141 /* BEGIN adopted by C99 */ 142 extern double erf(double); 143 extern double erfc(double); 144 extern double hypot(double, double); 145 extern double lgamma(double); 146 147 #if defined(__MATHERR_ERRNO_DONTCARE) 148 #pragma does_not_read_global_data(erf, erfc, hypot) 149 #pragma does_not_write_global_data(erf, erfc, hypot) 150 #pragma no_side_effect(erf, erfc, hypot) 151 #endif 152 153 #if !defined(_STDC_C99) && _XOPEN_SOURCE - 0 < 600 && !defined(__C99FEATURES__) 154 extern int isnan(double); 155 156 #pragma does_not_read_global_data(isnan) 157 #pragma does_not_write_global_data(isnan) 158 #pragma no_side_effect(isnan) 159 #endif 160 /* END adopted by C99 */ 161 162 #if defined(__EXTENSIONS__) || _XOPEN_SOURCE - 0 < 600 163 extern double gamma(double); /* deprecated; use lgamma */ 164 #endif 165 extern double j0(double); 166 extern double j1(double); 167 extern double jn(int, double); 168 extern double y0(double); 169 extern double y1(double); 170 extern double yn(int, double); 171 172 #if defined(__MATHERR_ERRNO_DONTCARE) 173 #pragma does_not_read_global_data(j0, j1, jn, y0, y1, yn) 174 #pragma does_not_write_global_data(j0, j1, jn, y0, y1, yn) 175 #pragma no_side_effect(j0, j1, jn, y0, y1, yn) 176 #endif 177 #if defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) || \ 178 _XOPEN_SOURCE - 0 >= 500 || \ 179 defined(_XOPEN_SOURCE) && _XOPEN_SOURCE_EXTENDED - 0 == 1 180 /* 181 * SVID & XPG 4.2/5 182 */ 183 extern double scalb(double, double); 184 185 #if defined(__MATHERR_ERRNO_DONTCARE) 186 #pragma does_not_read_global_data(scalb) 187 #pragma does_not_write_global_data(scalb) 188 #pragma no_side_effect(scalb) 189 #endif 190 191 /* BEGIN adopted by C99 */ 192 extern double acosh(double); 193 extern double asinh(double); 194 extern double atanh(double); 195 extern double cbrt(double); 196 extern double logb(double); 197 extern double nextafter(double, double); 198 extern double remainder(double, double); 199 200 /* 201 * XPG 4.2/5 202 */ 203 extern double expm1(double); 204 extern int ilogb(double); 205 extern double log1p(double); 206 extern double rint(double); 207 208 #if defined(__MATHERR_ERRNO_DONTCARE) 209 #pragma does_not_read_global_data(acosh, asinh, atanh, cbrt) 210 #pragma does_not_read_global_data(logb, nextafter, remainder) 211 #pragma does_not_read_global_data(expm1, ilogb, log1p, rint) 212 #pragma does_not_write_global_data(acosh, asinh, atanh, cbrt) 213 #pragma does_not_write_global_data(logb, nextafter, remainder) 214 #pragma does_not_write_global_data(expm1, ilogb, log1p, rint) 215 #pragma no_side_effect(acosh, asinh, atanh, cbrt) 216 #pragma no_side_effect(logb, nextafter, remainder) 217 #pragma no_side_effect(expm1, ilogb, log1p, rint) 218 #endif 219 /* END adopted by C99 */ 220 #endif /* defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) || ... */ 221 222 #if defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) 223 /* 224 * SVID 225 */ 226 extern int matherr(struct exception *); 227 228 /* 229 * IEEE Test Vector 230 */ 231 extern double significand(double); 232 233 #if defined(__MATHERR_ERRNO_DONTCARE) 234 #pragma does_not_read_global_data(significand) 235 #pragma does_not_write_global_data(significand) 236 #pragma no_side_effect(significand) 237 #endif 238 239 extern int signgamf; /* deprecated; use signgam */ 240 extern int signgaml; /* deprecated; use signgam */ 241 242 extern int isnanf(float); 243 extern int isnanl(long double); 244 extern float gammaf(float); /* deprecated; use lgammaf */ 245 extern float gammaf_r(float, int *); /* deprecated; use lgammaf_r */ 246 extern float j0f(float); 247 extern float j1f(float); 248 extern float jnf(int, float); 249 extern float lgammaf_r(float, int *); 250 extern float scalbf(float, float); 251 extern float significandf(float); 252 extern float y0f(float); 253 extern float y1f(float); 254 extern float ynf(int, float); 255 extern long double gammal(long double); /* deprecated; use lgammal */ 256 extern long double gammal_r(long double, int *); /* deprecated */ 257 extern long double j0l(long double); 258 extern long double j1l(long double); 259 extern long double jnl(int, long double); 260 extern long double lgammal_r(long double, int *); 261 extern long double scalbl(long double, long double); 262 extern long double significandl(long double); 263 extern long double y0l(long double); 264 extern long double y1l(long double); 265 extern long double ynl(int, long double); 266 267 #if defined(__MATHERR_ERRNO_DONTCARE) 268 #pragma does_not_read_global_data(isnanf, isnanl) 269 #pragma does_not_write_global_data(isnanf, isnanl) 270 #pragma no_side_effect(isnanf, isnanl) 271 #pragma does_not_read_global_data(gammaf_r, j0f, j1f, jnf, lgammaf_r, scalbf) 272 #pragma does_not_read_global_data(significandf, y0f, y1f, ynf) 273 #pragma does_not_write_global_data(j0f, j1f, jnf, scalbf) 274 #pragma does_not_write_global_data(significandf, y0f, y1f, ynf) 275 #pragma no_side_effect(j0f, j1f, jnf, scalbf) 276 #pragma no_side_effect(significandf, y0f, y1f, ynf) 277 #pragma does_not_read_global_data(gammal_r, j0l, j1l, jnl, lgammal_r, scalbl) 278 #pragma does_not_read_global_data(significandl, y0l, y1l, ynl) 279 #pragma does_not_write_global_data(j0l, j1l, jnl, scalbl) 280 #pragma does_not_write_global_data(significandl, y0l, y1l, ynl) 281 #pragma no_side_effect(j0l, j1l, jnl, scalbl) 282 #pragma no_side_effect(significandl, y0l, y1l, ynl) 283 #endif 284 285 /* 286 * for sin+cos->sincos transformation 287 */ 288 extern void sincos(double, double *, double *); 289 extern void sincosf(float, float *, float *); 290 extern void sincosl(long double, long double *, long double *); 291 292 #if defined(__MATHERR_ERRNO_DONTCARE) 293 #pragma does_not_read_global_data(sincos, sincosf, sincosl) 294 #endif 295 296 /* BEGIN adopted by C99 */ 297 /* 298 * Functions callable from C, intended to support IEEE arithmetic. 299 */ 300 extern double copysign(double, double); 301 extern double scalbn(double, int); 302 303 #if defined(__MATHERR_ERRNO_DONTCARE) 304 #pragma does_not_read_global_data(copysign, scalbn) 305 #pragma does_not_write_global_data(copysign, scalbn) 306 #pragma no_side_effect(copysign, scalbn) 307 #endif 308 /* END adopted by C99 */ 309 310 /* 311 * Reentrant version of gamma & lgamma; passes signgam back by reference 312 * as the second argument; user must allocate space for signgam. 313 */ 314 extern double gamma_r(double, int *); /* deprecated; use lgamma_r */ 315 extern double lgamma_r(double, int *); 316 317 #if defined(__MATHERR_ERRNO_DONTCARE) 318 #pragma does_not_read_global_data(gamma_r, lgamma_r) 319 #endif 320 321 /* BEGIN adopted by C99 */ 322 extern float modff(float, float *); 323 324 #if defined(__MATHERR_ERRNO_DONTCARE) 325 #pragma does_not_read_global_data(modff) 326 #endif 327 /* END adopted by C99 */ 328 329 #if defined(__EXTENSIONS__) || !defined(__cplusplus) 330 #include <floatingpoint.h> 331 #endif 332 #endif /* defined(__EXTENSIONS__) || !defined(_XOPEN_SOURCE) */ 333 #endif /* defined(__EXTENSIONS__) || defined(_XOPEN_SOURCE) || ... */ 334 335 #if defined(__cplusplus) && defined(__GNUC__) 336 #undef exception 337 #endif 338 339 #ifdef __cplusplus 340 } 341 #endif 342 343 #endif /* _MATH_H */ 344