diff options
author | bde <bde@FreeBSD.org> | 1997-02-20 13:35:15 +0000 |
---|---|---|
committer | bde <bde@FreeBSD.org> | 1997-02-20 13:35:15 +0000 |
commit | a9570d9f942c778c099c84d08f34773b14f0d9c8 (patch) | |
tree | 968171ff9ded709c4e63152a09c3e290b805168b /lib/msun/src/w_ynf.c | |
parent | 0071d9dabbe23115ee60715fbeba793f63803d49 (diff) | |
download | FreeBSD-src-a9570d9f942c778c099c84d08f34773b14f0d9c8.zip FreeBSD-src-a9570d9f942c778c099c84d08f34773b14f0d9c8.tar.gz |
Split up the Bessel function wrapper files so that most wrapper functions
are in their own file.
Diffstat (limited to 'lib/msun/src/w_ynf.c')
-rw-r--r-- | lib/msun/src/w_ynf.c | 50 |
1 files changed, 50 insertions, 0 deletions
diff --git a/lib/msun/src/w_ynf.c b/lib/msun/src/w_ynf.c new file mode 100644 index 0000000..f7a2498 --- /dev/null +++ b/lib/msun/src/w_ynf.c @@ -0,0 +1,50 @@ +/* w_ynf.c -- float version of w_yn.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 + +#include "math.h" +#include "math_private.h" + +#ifdef __STDC__ + float ynf(int n, float x) /* wrapper ynf */ +#else + float ynf(n,x) /* wrapper ynf */ + float x; int n; +#endif +{ +#ifdef _IEEE_LIBM + return __ieee754_ynf(n,x); +#else + float z; + z = __ieee754_ynf(n,x); + if(_LIB_VERSION == _IEEE_ || isnanf(x) ) return z; + if(x <= (float)0.0){ + if(x==(float)0.0) + /* d= -one/(x-x); */ + return (float)__kernel_standard((double)n,(double)x,112); + else + /* d = zero/(x-x); */ + return (float)__kernel_standard((double)n,(double)x,113); + } + if(x>(float)X_TLOSS) { + /* yn(x>X_TLOSS,n) */ + return (float)__kernel_standard((double)n,(double)x,139); + } else + return z; +#endif +} |