Home | History | Annotate | Line # | Download | only in include
math.h revision 1.25
      1 /*	$NetBSD: math.h,v 1.25 2003/04/28 23:16:13 bjh21 Exp $	*/
      2 
      3 /*
      4  * ====================================================
      5  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
      6  *
      7  * Developed at SunPro, a Sun Microsystems, Inc. business.
      8  * Permission to use, copy, modify, and distribute this
      9  * software is freely granted, provided that this notice
     10  * is preserved.
     11  * ====================================================
     12  */
     13 
     14 /*
     15  * @(#)fdlibm.h 5.1 93/09/24
     16  */
     17 
     18 #ifndef _MATH_H_
     19 #define _MATH_H_
     20 
     21 #include <sys/cdefs.h>
     22 #include <sys/featuretest.h>
     23 
     24 union __float_u {
     25 	unsigned char __dummy[sizeof(float)];
     26 	float __val;
     27 };
     28 
     29 union __double_u {
     30 	unsigned char __dummy[sizeof(double)];
     31 	double __val;
     32 };
     33 
     34 #include <machine/math.h>		/* may use __float_u or __double_u */
     35 
     36 /*
     37  * ANSI/POSIX
     38  */
     39 extern __const union __double_u __infinity;
     40 #define HUGE_VAL	__infinity.__val
     41 
     42 /*
     43  * ISO C99
     44  */
     45 #if defined(__HAVE_NANF) && \
     46     (!defined(_ANSI_SOURCE) && \
     47      (defined(_ISOC99_SOURCE) || (__STDC_VERSION__ - 0) >= 199901L || \
     48       defined(_NETBSD_SOURCE)))
     49 extern __const union __float_u __nanf;
     50 #define	NAN		__nanf.__val
     51 #endif /* __HAVE_NANF && (!_ANSI_SOURCE && ....) */
     52 
     53 /*
     54  * XOPEN/SVID
     55  */
     56 #if defined(_XOPEN_SOURCE) || defined(_NETBSD_SOURCE)
     57 #define	M_E		2.7182818284590452354	/* e */
     58 #define	M_LOG2E		1.4426950408889634074	/* log 2e */
     59 #define	M_LOG10E	0.43429448190325182765	/* log 10e */
     60 #define	M_LN2		0.69314718055994530942	/* log e2 */
     61 #define	M_LN10		2.30258509299404568402	/* log e10 */
     62 #define	M_PI		3.14159265358979323846	/* pi */
     63 #define	M_PI_2		1.57079632679489661923	/* pi/2 */
     64 #define	M_PI_4		0.78539816339744830962	/* pi/4 */
     65 #define	M_1_PI		0.31830988618379067154	/* 1/pi */
     66 #define	M_2_PI		0.63661977236758134308	/* 2/pi */
     67 #define	M_2_SQRTPI	1.12837916709551257390	/* 2/sqrt(pi) */
     68 #define	M_SQRT2		1.41421356237309504880	/* sqrt(2) */
     69 #define	M_SQRT1_2	0.70710678118654752440	/* 1/sqrt(2) */
     70 
     71 #define	MAXFLOAT	((float)3.40282346638528860e+38)
     72 extern int signgam;
     73 #endif /* _XOPEN_SOURCE || _NETBSD_SOURCE */
     74 
     75 #if defined(_NETBSD_SOURCE)
     76 enum fdversion {fdlibm_ieee = -1, fdlibm_svid, fdlibm_xopen, fdlibm_posix};
     77 
     78 #define _LIB_VERSION_TYPE enum fdversion
     79 #define _LIB_VERSION _fdlib_version
     80 
     81 /* if global variable _LIB_VERSION is not desirable, one may
     82  * change the following to be a constant by:
     83  *	#define _LIB_VERSION_TYPE const enum version
     84  * In that case, after one initializes the value _LIB_VERSION (see
     85  * s_lib_version.c) during compile time, it cannot be modified
     86  * in the middle of a program
     87  */
     88 extern  _LIB_VERSION_TYPE  _LIB_VERSION;
     89 
     90 #define _IEEE_  fdlibm_ieee
     91 #define _SVID_  fdlibm_svid
     92 #define _XOPEN_ fdlibm_xopen
     93 #define _POSIX_ fdlibm_posix
     94 
     95 #ifndef __cplusplus
     96 struct exception {
     97 	int type;
     98 	char *name;
     99 	double arg1;
    100 	double arg2;
    101 	double retval;
    102 };
    103 #endif
    104 
    105 #define	HUGE		MAXFLOAT
    106 
    107 /*
    108  * set X_TLOSS = pi*2**52, which is possibly defined in <values.h>
    109  * (one may replace the following line by "#include <values.h>")
    110  */
    111 
    112 #define X_TLOSS		1.41484755040568800000e+16
    113 
    114 #define	DOMAIN		1
    115 #define	SING		2
    116 #define	OVERFLOW	3
    117 #define	UNDERFLOW	4
    118 #define	TLOSS		5
    119 #define	PLOSS		6
    120 
    121 #endif /* _NETBSD_SOURCE */
    122 
    123 __BEGIN_DECLS
    124 /*
    125  * ANSI/POSIX
    126  */
    127 double	acos __P((double));
    128 double	asin __P((double));
    129 double	atan __P((double));
    130 double	atan2 __P((double, double));
    131 double	cos __P((double));
    132 double	sin __P((double));
    133 double	tan __P((double));
    134 
    135 double	cosh __P((double));
    136 double	sinh __P((double));
    137 double	tanh __P((double));
    138 
    139 double	exp __P((double));
    140 double	frexp __P((double, int *));
    141 double	ldexp __P((double, int));
    142 double	log __P((double));
    143 double	log10 __P((double));
    144 double	modf __P((double, double *));
    145 
    146 double	pow __P((double, double));
    147 double	sqrt __P((double));
    148 
    149 double	ceil __P((double));
    150 double	fabs __P((double));
    151 double	floor __P((double));
    152 double	fmod __P((double, double));
    153 
    154 #if defined(_XOPEN_SOURCE) || defined(_NETBSD_SOURCE)
    155 double	erf __P((double));
    156 double	erfc __P((double));
    157 double	gamma __P((double));
    158 double	hypot __P((double, double));
    159 int	isnan __P((double));
    160 int	finite __P((double));
    161 double	j0 __P((double));
    162 double	j1 __P((double));
    163 double	jn __P((int, double));
    164 double	lgamma __P((double));
    165 double	y0 __P((double));
    166 double	y1 __P((double));
    167 double	yn __P((int, double));
    168 
    169 #if (_XOPEN_SOURCE - 0) >= 500 || defined(_NETBSD_SOURCE)
    170 double	acosh __P((double));
    171 double	asinh __P((double));
    172 double	atanh __P((double));
    173 double	cbrt __P((double));
    174 double	expm1 __P((double));
    175 int	ilogb __P((double));
    176 double	log1p __P((double));
    177 double	logb __P((double));
    178 double	nextafter __P((double, double));
    179 double	remainder __P((double, double));
    180 double	rint __P((double));
    181 double	scalb __P((double, double));
    182 #endif /* (_XOPEN_SOURCE - 0) >= 500 || defined(_NETBSD_SOURCE)*/
    183 #endif /* _XOPEN_SOURCE || _NETBSD_SOURCE */
    184 
    185 #if defined(_NETBSD_SOURCE)
    186 #ifndef __cplusplus
    187 int	matherr __P((struct exception *));
    188 #endif
    189 
    190 /*
    191  * IEEE Test Vector
    192  */
    193 double	significand __P((double));
    194 
    195 /*
    196  * Functions callable from C, intended to support IEEE arithmetic.
    197  */
    198 double	copysign __P((double, double));
    199 double	scalbn __P((double, int));
    200 
    201 /*
    202  * BSD math library entry points
    203  */
    204 #ifndef __MATH_PRIVATE__
    205 double	cabs __P((/* struct complex { double r; double i; } */));
    206 #endif
    207 double	drem __P((double, double));
    208 
    209 #endif /* _NETBSD_SOURCE */
    210 
    211 #if defined(_NETBSD_SOURCE) || defined(_REENTRANT)
    212 /*
    213  * Reentrant version of gamma & lgamma; passes signgam back by reference
    214  * as the second argument; user must allocate space for signgam.
    215  */
    216 double	gamma_r __P((double, int *));
    217 double	lgamma_r __P((double, int *));
    218 #endif /* _NETBSD_SOURCE || _REENTRANT */
    219 
    220 
    221 #if defined(_NETBSD_SOURCE)
    222 int	isinf __P((double));
    223 
    224 /* float versions of ANSI/POSIX functions */
    225 float	acosf __P((float));
    226 float	asinf __P((float));
    227 float	atanf __P((float));
    228 float	atan2f __P((float, float));
    229 float	cosf __P((float));
    230 float	sinf __P((float));
    231 float	tanf __P((float));
    232 
    233 float	coshf __P((float));
    234 float	sinhf __P((float));
    235 float	tanhf __P((float));
    236 
    237 float	expf __P((float));
    238 float	frexpf __P((float, int *));
    239 float	ldexpf __P((float, int));
    240 float	logf __P((float));
    241 float	log10f __P((float));
    242 float	modff __P((float, float *));
    243 
    244 float	powf __P((float, float));
    245 float	sqrtf __P((float));
    246 
    247 float	ceilf __P((float));
    248 float	fabsf __P((float));
    249 float	floorf __P((float));
    250 float	fmodf __P((float, float));
    251 
    252 float	erff __P((float));
    253 float	erfcf __P((float));
    254 float	gammaf __P((float));
    255 float	hypotf __P((float, float));
    256 int	isinff __P((float));
    257 int	isnanf __P((float));
    258 int	finitef __P((float));
    259 float	j0f __P((float));
    260 float	j1f __P((float));
    261 float	jnf __P((int, float));
    262 float	lgammaf __P((float));
    263 float	y0f __P((float));
    264 float	y1f __P((float));
    265 float	ynf __P((int, float));
    266 
    267 float	acoshf __P((float));
    268 float	asinhf __P((float));
    269 float	atanhf __P((float));
    270 float	cbrtf __P((float));
    271 float	logbf __P((float));
    272 float	nextafterf __P((float, float));
    273 float	remainderf __P((float, float));
    274 float	scalbf __P((float, float));
    275 
    276 /*
    277  * float version of IEEE Test Vector
    278  */
    279 float	significandf __P((float));
    280 
    281 /*
    282  * Float versions of functions callable from C, intended to support
    283  * IEEE arithmetic.
    284  */
    285 float	copysignf __P((float, float));
    286 int	ilogbf __P((float));
    287 float	rintf __P((float));
    288 float	scalbnf __P((float, int));
    289 
    290 /*
    291  * float versions of BSD math library entry points
    292  */
    293 #ifndef __MATH_PRIVATE__
    294 float	cabsf __P((/* struct complex { float r; float i; } */));
    295 #endif
    296 float	dremf __P((float, float));
    297 float	expm1f __P((float));
    298 float	log1pf __P((float));
    299 #endif /* _NETBSD_SOURCE */
    300 
    301 #if defined(_NETBSD_SOURCE) || defined(_REENTRANT)
    302 /*
    303  * Float versions of reentrant version of gamma & lgamma; passes
    304  * signgam back by reference as the second argument; user must
    305  * allocate space for signgam.
    306  */
    307 float	gammaf_r __P((float, int *));
    308 float	lgammaf_r __P((float, int *));
    309 #endif /* !... || _REENTRANT */
    310 
    311 __END_DECLS
    312 
    313 #endif /* _MATH_H_ */
    314