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