1*12b1df16SDavid du Colombier #ifndef __MATH 2*12b1df16SDavid du Colombier #define __MATH 3*12b1df16SDavid du Colombier #pragma lib "/$M/lib/ape/libap.a" 4*12b1df16SDavid du Colombier 5*12b1df16SDavid du Colombier /* a HUGE_VAL appropriate for IEEE double-precision */ 6*12b1df16SDavid du Colombier /* the correct value, 1.797693134862316e+308, causes a ken overflow */ 7*12b1df16SDavid du Colombier #define HUGE_VAL 1.79769313486231e+308 8*12b1df16SDavid du Colombier 9*12b1df16SDavid du Colombier #ifdef __cplusplus 10*12b1df16SDavid du Colombier extern "C" { 11*12b1df16SDavid du Colombier #endif 12*12b1df16SDavid du Colombier 13*12b1df16SDavid du Colombier extern double acos(double); 14*12b1df16SDavid du Colombier extern double asin(double); 15*12b1df16SDavid du Colombier extern double atan(double); 16*12b1df16SDavid du Colombier extern double atan2(double, double); 17*12b1df16SDavid du Colombier extern double cos(double); 18*12b1df16SDavid du Colombier extern double hypot(double, double); 19*12b1df16SDavid du Colombier extern double sin(double); 20*12b1df16SDavid du Colombier extern double tan(double); 21*12b1df16SDavid du Colombier extern double cosh(double); 22*12b1df16SDavid du Colombier extern double sinh(double); 23*12b1df16SDavid du Colombier extern double tanh(double); 24*12b1df16SDavid du Colombier extern double exp(double); 25*12b1df16SDavid du Colombier extern double frexp(double, int *); 26*12b1df16SDavid du Colombier extern double ldexp(double, int); 27*12b1df16SDavid du Colombier extern double log(double); 28*12b1df16SDavid du Colombier extern double log10(double); 29*12b1df16SDavid du Colombier extern double modf(double, double *); 30*12b1df16SDavid du Colombier extern double pow(double, double); 31*12b1df16SDavid du Colombier extern double sqrt(double); 32*12b1df16SDavid du Colombier extern double ceil(double); 33*12b1df16SDavid du Colombier extern double fabs(double); 34*12b1df16SDavid du Colombier extern double floor(double); 35*12b1df16SDavid du Colombier extern double fmod(double, double); 36*12b1df16SDavid du Colombier extern double NaN(void); 37*12b1df16SDavid du Colombier extern int isNaN(double); 38*12b1df16SDavid du Colombier extern double Inf(int); 39*12b1df16SDavid du Colombier extern int isInf(double, int); 40*12b1df16SDavid du Colombier 41*12b1df16SDavid du Colombier #ifdef _RESEARCH_SOURCE 42*12b1df16SDavid du Colombier /* does >> treat left operand as unsigned ? */ 43*12b1df16SDavid du Colombier #define Unsigned_Shifts 1 44*12b1df16SDavid du Colombier #define M_E 2.7182818284590452354 /* e */ 45*12b1df16SDavid du Colombier #define M_LOG2E 1.4426950408889634074 /* log 2e */ 46*12b1df16SDavid du Colombier #define M_LOG10E 0.43429448190325182765 /* log 10e */ 47*12b1df16SDavid du Colombier #define M_LN2 0.69314718055994530942 /* log e2 */ 48*12b1df16SDavid du Colombier #define M_LN10 2.30258509299404568402 /* log e10 */ 49*12b1df16SDavid du Colombier #define M_PI 3.14159265358979323846 /* pi */ 50*12b1df16SDavid du Colombier #define M_PI_2 1.57079632679489661923 /* pi/2 */ 51*12b1df16SDavid du Colombier #define M_PI_4 0.78539816339744830962 /* pi/4 */ 52*12b1df16SDavid du Colombier #define M_1_PI 0.31830988618379067154 /* 1/pi */ 53*12b1df16SDavid du Colombier #define M_2_PI 0.63661977236758134308 /* 2/pi */ 54*12b1df16SDavid du Colombier #define M_2_SQRTPI 1.12837916709551257390 /* 2/sqrt(pi) */ 55*12b1df16SDavid du Colombier #define M_SQRT2 1.41421356237309504880 /* sqrt(2) */ 56*12b1df16SDavid du Colombier #define M_SQRT1_2 0.70710678118654752440 /* 1/sqrt(2) */ 57*12b1df16SDavid du Colombier 58*12b1df16SDavid du Colombier extern double hypot(double, double); 59*12b1df16SDavid du Colombier extern double erf(double); 60*12b1df16SDavid du Colombier extern double erfc(double); 61*12b1df16SDavid du Colombier extern double j0(double); 62*12b1df16SDavid du Colombier extern double y0(double); 63*12b1df16SDavid du Colombier extern double j1(double); 64*12b1df16SDavid du Colombier extern double y1(double); 65*12b1df16SDavid du Colombier extern double jn(int, double); 66*12b1df16SDavid du Colombier extern double yn(int, double); 67*12b1df16SDavid du Colombier 68*12b1df16SDavid du Colombier #endif 69*12b1df16SDavid du Colombier 70*12b1df16SDavid du Colombier 71*12b1df16SDavid du Colombier #ifdef __cplusplus 72*12b1df16SDavid du Colombier } 73*12b1df16SDavid du Colombier #endif 74*12b1df16SDavid du Colombier 75*12b1df16SDavid du Colombier #define isnan(x) isNaN(x) 76*12b1df16SDavid du Colombier #define isinf(x) isInf(x, 0) 77*12b1df16SDavid du Colombier 78*12b1df16SDavid du Colombier #endif /* __MATH */ 79