diff options
Diffstat (limited to 'libm/e_jn.c')
| -rw-r--r-- | libm/e_jn.c | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/libm/e_jn.c b/libm/e_jn.c index b12a7f33a..3825bead8 100644 --- a/libm/e_jn.c +++ b/libm/e_jn.c @@ -200,6 +200,23 @@ double attribute_hidden __ieee754_jn(int n, double x) if(sgn==1) return -b; else return b; } +/* + * wrapper jn(int n, double x) + */ +#ifndef _IEEE_LIBM +double jn(int n, double x) +{ + double z = __ieee754_jn(n, x); + if (_LIB_VERSION == _IEEE_ || isnan(x)) + return z; + if (fabs(x) > X_TLOSS) + return __kernel_standard((double)n, x, 38); /* jn(|x|>X_TLOSS,n) */ + return z; +} +#else +strong_alias(__ieee754_jn, jn) +#endif + double attribute_hidden __ieee754_yn(int n, double x) { int32_t i,hx,ix,lx; @@ -256,3 +273,26 @@ double attribute_hidden __ieee754_yn(int n, double x) } if(sign>0) return b; else return -b; } + +/* + * wrapper yn(int n, double x) + */ +#ifndef _IEEE_LIBM +double yn(int n, double x) /* wrapper yn */ +{ + double z = __ieee754_yn(n, x); + if (_LIB_VERSION == _IEEE_ || isnan(x)) + return z; + if (x <= 0.0) { + if(x == 0.0) /* d= -one/(x-x); */ + return __kernel_standard((double)n, x, 12); + /* d = zero/(x-x); */ + return __kernel_standard((double)n, x, 13); + } + if (x > X_TLOSS) + return __kernel_standard((double)n, x, 39); /* yn(x>X_TLOSS,n) */ + return z; +} +#else +strong_alias(__ieee754_yn, yn) +#endif |