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