diff options
Diffstat (limited to 'lib/msun/src/w_j0f.c')
-rw-r--r-- | lib/msun/src/w_j0f.c | 41 |
1 files changed, 41 insertions, 0 deletions
diff --git a/lib/msun/src/w_j0f.c b/lib/msun/src/w_j0f.c new file mode 100644 index 0000000..8530626 --- /dev/null +++ b/lib/msun/src/w_j0f.c @@ -0,0 +1,41 @@ +/* w_j0f.c -- float version of w_j0.c. + * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. + */ + +/* + * ==================================================== + * 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. + * ==================================================== + */ + +#ifndef lint +static char rcsid[] = "$FreeBSD$"; +#endif + +/* + * wrapper j0f(float x) + */ + +#include "math.h" +#include "math_private.h" + +float +j0f(float x) /* wrapper j0f */ +{ +#ifdef _IEEE_LIBM + return __ieee754_j0f(x); +#else + float z = __ieee754_j0f(x); + if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z; + if(fabsf(x)>(float)X_TLOSS) { + /* j0f(|x|>X_TLOSS) */ + return (float)__kernel_standard((double)x,(double)x,134); + } else + return z; +#endif +} |