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