double jn(int n, double x) /* wrapper jn */ { #ifdef _IEEE_LIBM return __ieee754_jn(n,x); #else double z; 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) */ } else return z; #endif }
void Math_jn(void *fp) { F_Math_jn *f; f = fp; *f->ret = __ieee754_jn(f->n, f->x); }
/* wrapper jn */ double jn (int n, double x) { if (__builtin_expect (isgreater (fabs (x), X_TLOSS), 0) && _LIB_VERSION != _IEEE_ && _LIB_VERSION != _POSIX_) /* jn(n,|x|>X_TLOSS) */ return __kernel_standard (n, x, 38); return __ieee754_jn (n, x); }