include/math.h

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