nkeynes@185: /* math.h -- Definitions for the math floating point package.  */
nkeynes@185: 
nkeynes@185: #ifndef  _MATH_H_
nkeynes@185: #ifdef __cplusplus
nkeynes@185: extern "C" {
nkeynes@185: #endif
nkeynes@185: #define  _MATH_H_
nkeynes@185: 
nkeynes@185: #include <sys/reent.h>
nkeynes@185: #include <machine/ieeefp.h>
nkeynes@185: #include "_ansi.h"
nkeynes@185: 
nkeynes@185: #ifndef HUGE_VAL
nkeynes@185: 
nkeynes@185: /* Define HUGE_VAL as infinity, unless HUGE_VAL is already defined
nkeynes@185:    (which might have been done by something like math-68881.h).  */
nkeynes@185: 
nkeynes@185: union __dmath
nkeynes@185: {
nkeynes@185:   __uint32_t i[2];
nkeynes@185:   double d;
nkeynes@185: };
nkeynes@185: 
nkeynes@185: extern const union __dmath __infinity;
nkeynes@185: 
nkeynes@185: #define HUGE_VAL (__infinity.d)
nkeynes@185: 
nkeynes@185: #endif /* ! defined (HUGE_VAL) */
nkeynes@185: 
nkeynes@185: /* Reentrant ANSI C functions.  */
nkeynes@185: 
nkeynes@185: #ifndef __math_68881
nkeynes@185: extern double atan _PARAMS((double));
nkeynes@185: extern double cos _PARAMS((double));
nkeynes@185: extern double sin _PARAMS((double));
nkeynes@185: extern double tan _PARAMS((double));
nkeynes@185: extern double tanh _PARAMS((double));
nkeynes@185: extern double frexp _PARAMS((double, int *));
nkeynes@185: extern double modf _PARAMS((double, double *));
nkeynes@185: extern double ceil _PARAMS((double));
nkeynes@185: extern double fabs _PARAMS((double));
nkeynes@185: extern double floor _PARAMS((double));
nkeynes@185: #endif /* ! defined (__math_68881) */
nkeynes@185: 
nkeynes@185: /* Non reentrant ANSI C functions.  */
nkeynes@185: 
nkeynes@185: #ifndef _REENT_ONLY
nkeynes@185: #ifndef __math_6881
nkeynes@185: extern double acos _PARAMS((double));
nkeynes@185: extern double asin _PARAMS((double));
nkeynes@185: extern double atan2 _PARAMS((double, double));
nkeynes@185: extern double cosh _PARAMS((double));
nkeynes@185: extern double sinh _PARAMS((double));
nkeynes@185: extern double exp _PARAMS((double));
nkeynes@185: extern double ldexp _PARAMS((double, int));
nkeynes@185: extern double log _PARAMS((double));
nkeynes@185: extern double log10 _PARAMS((double));
nkeynes@185: extern double pow _PARAMS((double, double));
nkeynes@185: extern double sqrt _PARAMS((double));
nkeynes@185: extern double fmod _PARAMS((double, double));
nkeynes@185: #endif /* ! defined (__math_68881) */
nkeynes@185: #endif /* ! defined (_REENT_ONLY) */
nkeynes@185: 
nkeynes@185: #ifndef __STRICT_ANSI__
nkeynes@185: 
nkeynes@185: /* Non ANSI double precision functions.  */
nkeynes@185: 
nkeynes@185: extern double infinity _PARAMS((void));
nkeynes@185: extern double nan _PARAMS((void));
nkeynes@185: extern int isnan _PARAMS((double));
nkeynes@185: extern int isinf _PARAMS((double));
nkeynes@185: extern int finite _PARAMS((double));
nkeynes@185: extern double copysign _PARAMS((double, double));
nkeynes@185: extern int ilogb _PARAMS((double));
nkeynes@185: 
nkeynes@185: extern double asinh _PARAMS((double));
nkeynes@185: extern double cbrt _PARAMS((double));
nkeynes@185: extern double nextafter _PARAMS((double, double));
nkeynes@185: extern double rint _PARAMS((double));
nkeynes@185: extern double scalbn _PARAMS((double, int));
nkeynes@185: 
nkeynes@185: #ifndef __math_68881
nkeynes@185: extern double log1p _PARAMS((double));
nkeynes@185: extern double expm1 _PARAMS((double));
nkeynes@185: #endif /* ! defined (__math_68881) */
nkeynes@185: 
nkeynes@185: #ifndef _REENT_ONLY
nkeynes@185: extern double acosh _PARAMS((double));
nkeynes@185: extern double atanh _PARAMS((double));
nkeynes@185: extern double remainder _PARAMS((double, double));
nkeynes@185: extern double gamma _PARAMS((double));
nkeynes@185: extern double gamma_r _PARAMS((double, int *));
nkeynes@185: extern double lgamma _PARAMS((double));
nkeynes@185: extern double lgamma_r _PARAMS((double, int *));
nkeynes@185: extern double erf _PARAMS((double));
nkeynes@185: extern double erfc _PARAMS((double));
nkeynes@185: extern double y0 _PARAMS((double));
nkeynes@185: extern double y1 _PARAMS((double));
nkeynes@185: extern double yn _PARAMS((int, double));
nkeynes@185: extern double j0 _PARAMS((double));
nkeynes@185: extern double j1 _PARAMS((double));
nkeynes@185: extern double jn _PARAMS((int, double));
nkeynes@185: #define log2(x) (log (x) / M_LOG2_E)
nkeynes@185: 
nkeynes@185: #ifndef __math_68881
nkeynes@185: extern double hypot _PARAMS((double, double));
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: extern double cabs();
nkeynes@185: extern double drem _PARAMS((double, double));
nkeynes@185: 
nkeynes@185: #endif /* ! defined (_REENT_ONLY) */
nkeynes@185: 
nkeynes@185: /* Single precision versions of ANSI functions.  */
nkeynes@185: 
nkeynes@185: extern float atanf _PARAMS((float));
nkeynes@185: extern float cosf _PARAMS((float));
nkeynes@185: extern float sinf _PARAMS((float));
nkeynes@185: extern float tanf _PARAMS((float));
nkeynes@185: extern float tanhf _PARAMS((float));
nkeynes@185: extern float frexpf _PARAMS((float, int *));
nkeynes@185: extern float modff _PARAMS((float, float *));
nkeynes@185: extern float ceilf _PARAMS((float));
nkeynes@185: extern float fabsf _PARAMS((float));
nkeynes@185: extern float floorf _PARAMS((float));
nkeynes@185: 
nkeynes@185: #ifndef _REENT_ONLY
nkeynes@185: extern float acosf _PARAMS((float));
nkeynes@185: extern float asinf _PARAMS((float));
nkeynes@185: extern float atan2f _PARAMS((float, float));
nkeynes@185: extern float coshf _PARAMS((float));
nkeynes@185: extern float sinhf _PARAMS((float));
nkeynes@185: extern float expf _PARAMS((float));
nkeynes@185: extern float ldexpf _PARAMS((float, int));
nkeynes@185: extern float logf _PARAMS((float));
nkeynes@185: extern float log10f _PARAMS((float));
nkeynes@185: extern float powf _PARAMS((float, float));
nkeynes@185: extern float sqrtf _PARAMS((float));
nkeynes@185: extern float fmodf _PARAMS((float, float));
nkeynes@185: #endif /* ! defined (_REENT_ONLY) */
nkeynes@185: 
nkeynes@185: /* Other single precision functions.  */
nkeynes@185: 
nkeynes@185: extern float infinityf _PARAMS((void));
nkeynes@185: extern float nanf _PARAMS((void));
nkeynes@185: extern int isnanf _PARAMS((float));
nkeynes@185: extern int isinff _PARAMS((float));
nkeynes@185: extern int finitef _PARAMS((float));
nkeynes@185: extern float copysignf _PARAMS((float, float));
nkeynes@185: extern int ilogbf _PARAMS((float));
nkeynes@185: 
nkeynes@185: extern float asinhf _PARAMS((float));
nkeynes@185: extern float cbrtf _PARAMS((float));
nkeynes@185: extern float nextafterf _PARAMS((float, float));
nkeynes@185: extern float rintf _PARAMS((float));
nkeynes@185: extern float scalbnf _PARAMS((float, int));
nkeynes@185: extern float log1pf _PARAMS((float));
nkeynes@185: extern float expm1f _PARAMS((float));
nkeynes@185: 
nkeynes@185: #ifndef _REENT_ONLY
nkeynes@185: extern float acoshf _PARAMS((float));
nkeynes@185: extern float atanhf _PARAMS((float));
nkeynes@185: extern float remainderf _PARAMS((float, float));
nkeynes@185: extern float gammaf _PARAMS((float));
nkeynes@185: extern float gammaf_r _PARAMS((float, int *));
nkeynes@185: extern float lgammaf _PARAMS((float));
nkeynes@185: extern float lgammaf_r _PARAMS((float, int *));
nkeynes@185: extern float erff _PARAMS((float));
nkeynes@185: extern float erfcf _PARAMS((float));
nkeynes@185: extern float y0f _PARAMS((float));
nkeynes@185: extern float y1f _PARAMS((float));
nkeynes@185: extern float ynf _PARAMS((int, float));
nkeynes@185: extern float j0f _PARAMS((float));
nkeynes@185: extern float j1f _PARAMS((float));
nkeynes@185: extern float jnf _PARAMS((int, float));
nkeynes@185: #define log2f(x) (logf (x) / (float) M_LOG2_E)
nkeynes@185: extern float hypotf _PARAMS((float, float));
nkeynes@185: 
nkeynes@185: extern float cabsf();
nkeynes@185: extern float dremf _PARAMS((float, float));
nkeynes@185: 
nkeynes@185: #endif /* ! defined (_REENT_ONLY) */
nkeynes@185: 
nkeynes@185: /* The gamma functions use a global variable, signgam.  */
nkeynes@185: 
nkeynes@185: extern int signgam;
nkeynes@185: 
nkeynes@185: /* The exception structure passed to the matherr routine.  */
nkeynes@185: 
nkeynes@185: #ifdef __cplusplus
nkeynes@185: struct __exception 
nkeynes@185: #else
nkeynes@185: struct exception 
nkeynes@185: #endif
nkeynes@185: {
nkeynes@185:   int type;
nkeynes@185:   char *name;
nkeynes@185:   double arg1;
nkeynes@185:   double arg2;
nkeynes@185:   double retval;
nkeynes@185:   int err;
nkeynes@185: };
nkeynes@185: 
nkeynes@185: #ifdef __cplusplus
nkeynes@185: extern int matherr _PARAMS((struct __exception *e));
nkeynes@185: #else
nkeynes@185: extern int matherr _PARAMS((struct exception *e));
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: /* Values for the type field of struct exception.  */
nkeynes@185: 
nkeynes@185: #define DOMAIN 1
nkeynes@185: #define SING 2
nkeynes@185: #define OVERFLOW 3
nkeynes@185: #define UNDERFLOW 4
nkeynes@185: #define TLOSS 5
nkeynes@185: #define PLOSS 6
nkeynes@185: 
nkeynes@185: /* Useful constants.  */
nkeynes@185: 
nkeynes@185: #define M_E		2.7182818284590452354
nkeynes@185: #define M_LOG2E		1.4426950408889634074
nkeynes@185: #define M_LOG10E	0.43429448190325182765
nkeynes@185: #define M_LN2		0.69314718055994530942
nkeynes@185: #define M_LN10		2.30258509299404568402
nkeynes@185: #define M_PI		3.14159265358979323846
nkeynes@185: #define M_TWOPI         (M_PI * 2.0)
nkeynes@185: #define M_PI_2		1.57079632679489661923
nkeynes@185: #define M_PI_4		0.78539816339744830962
nkeynes@185: #define M_3PI_4		2.3561944901923448370E0
nkeynes@185: #define M_SQRTPI        1.77245385090551602792981
nkeynes@185: #define M_1_PI		0.31830988618379067154
nkeynes@185: #define M_2_PI		0.63661977236758134308
nkeynes@185: #define M_2_SQRTPI	1.12837916709551257390
nkeynes@185: #define M_SQRT2		1.41421356237309504880
nkeynes@185: #define M_SQRT1_2	0.70710678118654752440
nkeynes@185: #define M_LN2LO         1.9082149292705877000E-10
nkeynes@185: #define M_LN2HI         6.9314718036912381649E-1
nkeynes@185: #define M_SQRT3   	1.73205080756887719000
nkeynes@185: #define M_IVLN10        0.43429448190325182765 /* 1 / log(10) */
nkeynes@185: #define M_LOG2_E        0.693147180559945309417
nkeynes@185: #define M_INVLN2        1.4426950408889633870E0  /* 1 / log(2) */
nkeynes@185: 
nkeynes@185: /* Global control over fdlibm error handling.  */
nkeynes@185: 
nkeynes@185: enum __fdlibm_version
nkeynes@185: {
nkeynes@185:   __fdlibm_ieee = -1,
nkeynes@185:   __fdlibm_svid,
nkeynes@185:   __fdlibm_xopen,
nkeynes@185:   __fdlibm_posix
nkeynes@185: };
nkeynes@185: 
nkeynes@185: #define _LIB_VERSION_TYPE enum __fdlibm_version
nkeynes@185: #define _LIB_VERSION __fdlib_version
nkeynes@185: 
nkeynes@185: extern _CONST _LIB_VERSION_TYPE _LIB_VERSION;
nkeynes@185: 
nkeynes@185: #define _IEEE_  __fdlibm_ieee
nkeynes@185: #define _SVID_  __fdlibm_svid
nkeynes@185: #define _XOPEN_ __fdlibm_xopen
nkeynes@185: #define _POSIX_ __fdlibm_posix
nkeynes@185: 
nkeynes@185: #endif /* ! defined (__STRICT_ANSI__) */
nkeynes@185: 
nkeynes@185: #ifdef __FAST_MATH__
nkeynes@185: #include <machine/fastmath.h>
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: #ifdef __cplusplus
nkeynes@185: }
nkeynes@185: #endif
nkeynes@185: #endif /* _MATH_H_ */