inc/include/math.h

*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * CDDL HEADER START
*0Sstevel@tonic-gate *
*0Sstevel@tonic-gate * The contents of this file are subject to the terms of the
*0Sstevel@tonic-gate * Common Development and Distribution License, Version 1.0 only
*0Sstevel@tonic-gate * (the "License").  You may not use this file except in compliance
*0Sstevel@tonic-gate * with the License.
*0Sstevel@tonic-gate *
*0Sstevel@tonic-gate * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
*0Sstevel@tonic-gate * or http://www.opensolaris.org/os/licensing.
*0Sstevel@tonic-gate * See the License for the specific language governing permissions
*0Sstevel@tonic-gate * and limitations under the License.
*0Sstevel@tonic-gate *
*0Sstevel@tonic-gate * When distributing Covered Code, include this CDDL HEADER in each
*0Sstevel@tonic-gate * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
*0Sstevel@tonic-gate * If applicable, add the following below this CDDL HEADER, with the
*0Sstevel@tonic-gate * fields enclosed by brackets "[]" replaced with your own identifying
*0Sstevel@tonic-gate * information: Portions Copyright [yyyy] [name of copyright owner]
*0Sstevel@tonic-gate *
*0Sstevel@tonic-gate * CDDL HEADER END
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Copyright (c) 1988 by Sun Microsystems, Inc.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*      Copyright (c) 1984 AT&T */
*0Sstevel@tonic-gate/*        All Rights Reserved   */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#pragma ident	"%Z%%M%	%I%	%E% SMI"
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Math library definitions for all the public functions implemented in libm.a.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#ifndef	__math_h
*0Sstevel@tonic-gate#define	__math_h
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Posix (actually ansi C) section
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate#define	HUGE_VAL	(__infinity())	/* Produces IEEE Infinity. */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate
*0Sstevel@tonic-gateextern double	__infinity();
*0Sstevel@tonic-gateextern double	acos(/* double x */);
*0Sstevel@tonic-gateextern double	asin(/* double x */);
*0Sstevel@tonic-gateextern double	atan(/* double x */);
*0Sstevel@tonic-gateextern double	atan2(/* double y, double x */);
*0Sstevel@tonic-gateextern double	ceil(/* double x */);
*0Sstevel@tonic-gateextern double	cos(/* double x */);
*0Sstevel@tonic-gateextern double	cosh(/* double x */);
*0Sstevel@tonic-gateextern double	exp(/* double x */);
*0Sstevel@tonic-gateextern double	fabs(/* double x */);
*0Sstevel@tonic-gateextern double	floor(/* double x */);
*0Sstevel@tonic-gateextern double	fmod(/* double x, double y */);
*0Sstevel@tonic-gateextern double	frexp(/* double value, int *exp */);
*0Sstevel@tonic-gateextern double	ldexp(/* double value, int exp */);
*0Sstevel@tonic-gateextern double	log(/* double x */);
*0Sstevel@tonic-gateextern double	log10(/* double x */);
*0Sstevel@tonic-gateextern double	modf(/* double value, double *iptr */);
*0Sstevel@tonic-gateextern double	pow(/* double x, double y */);
*0Sstevel@tonic-gateextern double	sin(/* double x */);
*0Sstevel@tonic-gateextern double	sinh(/* double x */);
*0Sstevel@tonic-gateextern double	sqrt(/* double x */);
*0Sstevel@tonic-gateextern double	tan(/* double x */);
*0Sstevel@tonic-gateextern double	tanh(/* double x */);
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#ifndef	_POSIX_SOURCE			/* the rest of the file is !POSIX */
*0Sstevel@tonic-gate#include <floatingpoint.h>		/* Contains definitions for types and
*0Sstevel@tonic-gate					 * functions implemented in libc.a.
*0Sstevel@tonic-gate					 */
*0Sstevel@tonic-gateextern double	acosh();
*0Sstevel@tonic-gateextern double	asinh();
*0Sstevel@tonic-gateextern double	atanh();
*0Sstevel@tonic-gateextern double	cbrt();
*0Sstevel@tonic-gateextern double	copysign();
*0Sstevel@tonic-gateextern double	erf();
*0Sstevel@tonic-gateextern double	erfc();
*0Sstevel@tonic-gateextern double	expm1();
*0Sstevel@tonic-gateextern int	finite();
*0Sstevel@tonic-gateextern double	hypot();
*0Sstevel@tonic-gateextern double	j0();
*0Sstevel@tonic-gateextern double	j1();
*0Sstevel@tonic-gateextern double	jn();
*0Sstevel@tonic-gateextern double	lgamma();
*0Sstevel@tonic-gateextern double	log1p();
*0Sstevel@tonic-gateextern double	rint();
*0Sstevel@tonic-gateextern double	y0();
*0Sstevel@tonic-gateextern double	y1();
*0Sstevel@tonic-gateextern double	yn();
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Sun definitions.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Implemented precisions for trigonometric argument reduction.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gateenum fp_pi_type {
*0Sstevel@tonic-gate	fp_pi_infinite	= 0,	/* Infinite-precision approximation to pi. */
*0Sstevel@tonic-gate	fp_pi_66	= 1,	/* 66-bit approximation to pi. */
*0Sstevel@tonic-gate	fp_pi_53	= 2	/* 53-bit approximation to pi. */
*0Sstevel@tonic-gate};
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Pi precision to use for trigonometric argument reduction.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gateextern enum	fp_pi_type fp_pi;
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Functions callable from C, intended to support IEEE arithmetic.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gateextern enum	fp_class_type fp_class();
*0Sstevel@tonic-gateextern int	ieee_flags();
*0Sstevel@tonic-gateextern int	ieee_handler();
*0Sstevel@tonic-gateextern void	ieee_retrospective();
*0Sstevel@tonic-gateextern int	ilogb();
*0Sstevel@tonic-gateextern double	infinity();
*0Sstevel@tonic-gateextern int	irint();
*0Sstevel@tonic-gateextern int	isinf();
*0Sstevel@tonic-gateextern int	isnan();
*0Sstevel@tonic-gateextern int	isnormal();
*0Sstevel@tonic-gateextern int	issubnormal();
*0Sstevel@tonic-gateextern int	iszero();
*0Sstevel@tonic-gateextern double	logb();
*0Sstevel@tonic-gateextern double	max_normal();
*0Sstevel@tonic-gateextern double	max_subnormal();
*0Sstevel@tonic-gateextern double	min_normal();
*0Sstevel@tonic-gateextern double	min_subnormal();
*0Sstevel@tonic-gateextern double	nextafter();
*0Sstevel@tonic-gateextern void	nonstandard_arithmetic();
*0Sstevel@tonic-gateextern double	quiet_nan();
*0Sstevel@tonic-gateextern double	remainder();
*0Sstevel@tonic-gateextern double	scalb();
*0Sstevel@tonic-gateextern double	scalbn();
*0Sstevel@tonic-gateextern double	signaling_nan();
*0Sstevel@tonic-gateextern int	signbit();
*0Sstevel@tonic-gateextern double	significand();
*0Sstevel@tonic-gateextern void	standard_arithmetic();
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * Other functions for C programmers.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gateextern double	acospi();
*0Sstevel@tonic-gateextern double	aint();
*0Sstevel@tonic-gateextern double	anint();
*0Sstevel@tonic-gateextern double	annuity();
*0Sstevel@tonic-gateextern double	asinpi();
*0Sstevel@tonic-gateextern double	atan2pi();
*0Sstevel@tonic-gateextern double	atanpi();
*0Sstevel@tonic-gateextern double	compound();
*0Sstevel@tonic-gateextern double	cospi();
*0Sstevel@tonic-gateextern double	exp10();
*0Sstevel@tonic-gateextern double	exp2();
*0Sstevel@tonic-gateextern double	log2();
*0Sstevel@tonic-gateextern int	nint();
*0Sstevel@tonic-gateextern void	sincos();
*0Sstevel@tonic-gateextern void	sincospi();
*0Sstevel@tonic-gateextern double	sinpi();
*0Sstevel@tonic-gateextern double	tanpi();
*0Sstevel@tonic-gateextern int	matherr();
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate *	Single-precision functions callable from Fortran, Pascal, Modula-2, etc,
*0Sstevel@tonic-gate *	take float* arguments instead of double and
*0Sstevel@tonic-gate *	return FLOATFUNCTIONTYPE results instead of double.
*0Sstevel@tonic-gate *	RETURNFLOAT is used to return a float function value without conversion
*0Sstevel@tonic-gate *	to double.
*0Sstevel@tonic-gate *	ASSIGNFLOAT is used to get the float value out of a FLOATFUNCTIONTYPE
*0Sstevel@tonic-gate *	result.
*0Sstevel@tonic-gate *	We don't want you to have to think about -fsingle2.
*0Sstevel@tonic-gate *
*0Sstevel@tonic-gate *	Some internal library functions pass float parameters as 32-bit values,
*0Sstevel@tonic-gate *	disguised as FLOATPARAMETER.  FLOATPARAMETERVALUE(x) extracts the
*0Sstevel@tonic-gate *	float value from the FLOATPARAMETER.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*	mc68000 returns float results in d0, same as int	*/
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#ifdef	mc68000
*0Sstevel@tonic-gate#define	FLOATFUNCTIONTYPE	int
*0Sstevel@tonic-gate#define	RETURNFLOAT(x) 		return (*(int *)(&(x)))
*0Sstevel@tonic-gate#define	ASSIGNFLOAT(x,y)	*(int *)(&x) = y
*0Sstevel@tonic-gate#endif
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*	sparc returns float results in %f0, same as top half of double	*/
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#ifdef	sparc
*0Sstevel@tonic-gate#define	FLOATFUNCTIONTYPE	double
*0Sstevel@tonic-gate#define	RETURNFLOAT(x) 		{ union {double _d; float _f } _kluge; _kluge._f = (x); return _kluge._d; }
*0Sstevel@tonic-gate#define	ASSIGNFLOAT(x,y)	{ union {double _d; float _f } _kluge; _kluge._d = (y); x = _kluge._f; }
*0Sstevel@tonic-gate#endif
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*	i386 returns float results on stack as extendeds, same as double */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#ifdef	i386
*0Sstevel@tonic-gate#define	FLOATFUNCTIONTYPE	float
*0Sstevel@tonic-gate#define	RETURNFLOAT(x) 		return (x)
*0Sstevel@tonic-gate#define	ASSIGNFLOAT(x,y)	x = y
*0Sstevel@tonic-gate#endif
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/* So far everybody passes float parameters as 32 bits on stack, same as int. */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#define	FLOATPARAMETER		int
*0Sstevel@tonic-gate#define	FLOATPARAMETERVALUE(x)	(*(float *)(&(x)))
*0Sstevel@tonic-gate
*0Sstevel@tonic-gateextern int		 ir_finite_();
*0Sstevel@tonic-gateextern enum fp_class_type ir_fp_class_();
*0Sstevel@tonic-gateextern int		 ir_ilogb_();
*0Sstevel@tonic-gateextern int		 ir_irint_();
*0Sstevel@tonic-gateextern int		 ir_isinf_();
*0Sstevel@tonic-gateextern int		 ir_isnan_();
*0Sstevel@tonic-gateextern int		 ir_isnormal_();
*0Sstevel@tonic-gateextern int		 ir_issubnormal_();
*0Sstevel@tonic-gateextern int		 ir_iszero_();
*0Sstevel@tonic-gateextern int		 ir_nint_();
*0Sstevel@tonic-gateextern int		 ir_signbit_();
*0Sstevel@tonic-gateextern void		 r_sincos_();
*0Sstevel@tonic-gateextern void		 r_sincospi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_acos_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_acosh_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_acospi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_aint_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_anint_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_annuity_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_asin_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_asinh_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_asinpi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_atan2_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_atan2pi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_atan_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_atanh_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_atanpi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_cbrt_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_ceil_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_compound_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_copysign_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_cos_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_cosh_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_cospi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_erf_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_erfc_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_exp10_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_exp2_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_exp_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_expm1_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_fabs_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_floor_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_fmod_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_hypot_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_infinity_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_j0_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_j1_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_jn_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_lgamma_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_log10_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_log1p_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_log2_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_log_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_logb_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_max_normal_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_max_subnormal_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_min_normal_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_min_subnormal_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_nextafter_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_pow_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_quiet_nan_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_remainder_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_rint_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_scalb_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_scalbn_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_signaling_nan_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_significand_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_sin_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_sinh_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_sinpi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_sqrt_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_tan_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_tanh_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_tanpi_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_y0_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_y1_();
*0Sstevel@tonic-gateextern FLOATFUNCTIONTYPE r_yn_();
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/* 	Constants, variables, and functions from System V */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#define	_ABS(x) ((x) < 0 ? -(x) : (x))
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#define	HUGE		(infinity())	/* For historical compatibility. */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#define	DOMAIN		1
*0Sstevel@tonic-gate#define	SING		2
*0Sstevel@tonic-gate#define	OVERFLOW	3
*0Sstevel@tonic-gate#define	UNDERFLOW	4
*0Sstevel@tonic-gate#define	TLOSS		5
*0Sstevel@tonic-gate#define	PLOSS		6
*0Sstevel@tonic-gate
*0Sstevel@tonic-gatestruct exception {
*0Sstevel@tonic-gate	int type;
*0Sstevel@tonic-gate	char *name;
*0Sstevel@tonic-gate	double arg1;
*0Sstevel@tonic-gate	double arg2;
*0Sstevel@tonic-gate	double retval;
*0Sstevel@tonic-gate};
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate * First three have to be defined exactly as in values.h including spacing!
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate#define	M_LN2	0.69314718055994530942
*0Sstevel@tonic-gate#define	M_PI	3.14159265358979323846
*0Sstevel@tonic-gate#define	M_SQRT2	1.41421356237309504880
*0Sstevel@tonic-gate
*0Sstevel@tonic-gate#define	M_E		2.7182818284590452354
*0Sstevel@tonic-gate#define	M_LOG2E		1.4426950408889634074
*0Sstevel@tonic-gate#define	M_LOG10E	0.43429448190325182765
*0Sstevel@tonic-gate#define	M_LN10		2.30258509299404568402
*0Sstevel@tonic-gate#define	M_PI_2		1.57079632679489661923
*0Sstevel@tonic-gate#define	M_PI_4		0.78539816339744830962
*0Sstevel@tonic-gate#define	M_1_PI		0.31830988618379067154
*0Sstevel@tonic-gate#define	M_2_PI		0.63661977236758134308
*0Sstevel@tonic-gate#define	M_2_SQRTPI	1.12837916709551257390
*0Sstevel@tonic-gate#define	M_SQRT1_2	0.70710678118654752440
*0Sstevel@tonic-gate#define	_REDUCE(TYPE, X, XN, C1, C2)	{ \
*0Sstevel@tonic-gate	double x1 = (double)(TYPE)X, x2 = X - x1; \
*0Sstevel@tonic-gate	X = x1 - (XN) * (C1); X += x2; X -= (XN) * (C2); }
*0Sstevel@tonic-gate#define	_POLY1(x, c)    ((c)[0] * (x) + (c)[1])
*0Sstevel@tonic-gate#define	_POLY2(x, c)    (_POLY1((x), (c)) * (x) + (c)[2])
*0Sstevel@tonic-gate#define	_POLY3(x, c)    (_POLY2((x), (c)) * (x) + (c)[3])
*0Sstevel@tonic-gate#define	_POLY4(x, c)    (_POLY3((x), (c)) * (x) + (c)[4])
*0Sstevel@tonic-gate#define	_POLY5(x, c)    (_POLY4((x), (c)) * (x) + (c)[5])
*0Sstevel@tonic-gate#define	_POLY6(x, c)    (_POLY5((x), (c)) * (x) + (c)[6])
*0Sstevel@tonic-gate#define	_POLY7(x, c)    (_POLY6((x), (c)) * (x) + (c)[7])
*0Sstevel@tonic-gate#define	_POLY8(x, c)    (_POLY7((x), (c)) * (x) + (c)[8])
*0Sstevel@tonic-gate#define	_POLY9(x, c)    (_POLY8((x), (c)) * (x) + (c)[9])
*0Sstevel@tonic-gate
*0Sstevel@tonic-gateextern int	signgam;
*0Sstevel@tonic-gate/*
*0Sstevel@tonic-gate *	Deprecated functions for compatibility with past.
*0Sstevel@tonic-gate *	Changes planned for future.
*0Sstevel@tonic-gate */
*0Sstevel@tonic-gate
*0Sstevel@tonic-gateextern double cabs();	/* Use double hypot(x,y)
*0Sstevel@tonic-gate			 * Traditional cabs usage is confused -
*0Sstevel@tonic-gate			 * is its argument two doubles or one struct?
*0Sstevel@tonic-gate			 */
*0Sstevel@tonic-gateextern double drem();	/* Use double remainder(x,y)
*0Sstevel@tonic-gate			 * drem will disappear in a future release.
*0Sstevel@tonic-gate			 */
*0Sstevel@tonic-gateextern double gamma();	/* Use double lgamma(x)
*0Sstevel@tonic-gate			 * to compute log of gamma function.
*0Sstevel@tonic-gate			 * Name gamma is reserved for true gamma function
*0Sstevel@tonic-gate			 * to appear in a future release.
*0Sstevel@tonic-gate			 */
*0Sstevel@tonic-gate#endif	/* !_POSIX_SOURCE */
*0Sstevel@tonic-gate#endif	/* !__math_h */