/* @(#)w_exp.c 5.1 93/09/24 */ /* * ==================================================== * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. * * Developed at SunPro, a Sun Microsystems, Inc. business. * Permission to use, copy, modify, and distribute this * software is freely granted, provided that this notice * is preserved. * ==================================================== */ /* FUNCTION <>, <>---exponential INDEX exp INDEX expf ANSI_SYNOPSIS #include double exp(double <[x]>); float expf(float <[x]>); TRAD_SYNOPSIS #include double exp(<[x]>); double <[x]>; float expf(<[x]>); float <[x]>; DESCRIPTION <> and <> calculate the exponential of <[x]>, that is, @ifinfo e raised to the power <[x]> (where e @end ifinfo @tex $e^x$ (where $e$ @end tex is the base of the natural system of logarithms, approximately 2.71828). You can use the (non-ANSI) function <> to specify error handling for these functions. RETURNS On success, <> and <> return the calculated value. If the result underflows, the returned value is <<0>>. If the result overflows, the returned value is <>. In either case, <> is set to <>. PORTABILITY <> is ANSI C. <> is an extension. */ /* * wrapper exp(x) */ #include "fdlibm.h" #include #ifndef _DOUBLE_IS_32BITS #ifdef __STDC__ static const double #else static double #endif o_threshold= 7.09782712893383973096e+02, /* 0x40862E42, 0xFEFA39EF */ u_threshold= -7.45133219101941108420e+02; /* 0xc0874910, 0xD52D3051 */ #ifdef __STDC__ double exp(double x) /* wrapper exp */ #else double exp(x) /* wrapper exp */ double x; #endif { #ifdef _IEEE_LIBM return __ieee754_exp(x); #else double z; struct exception exc; z = __ieee754_exp(x); if(_LIB_VERSION == _IEEE_) return z; if(finite(x)) { if(x>o_threshold) { /* exp(finite) overflow */ #ifndef HUGE_VAL #define HUGE_VAL inf double inf = 0.0; SET_HIGH_WORD(inf,0x7ff00000); /* set inf to infinite */ #endif exc.type = OVERFLOW; exc.name = "exp"; exc.err = 0; exc.arg1 = exc.arg2 = x; if (_LIB_VERSION == _SVID_) exc.retval = HUGE; else exc.retval = HUGE_VAL; if (_LIB_VERSION == _POSIX_) errno = ERANGE; else if (!matherr(&exc)) { errno = ERANGE; } if (exc.err != 0) errno = exc.err; return exc.retval; } else if(x