void Math_y1(void *fp) { F_Math_y1 *f; f = fp; *f->ret = __ieee754_y1(f->x); }
/* * @implemented */ double _y1(double num) { double retval; if (!_finite(num)) *_errno() = EDOM; retval = __ieee754_y1(num); if (_fpclass(retval) == _FPCLASS_NINF) { *_errno() = EDOM; retval = sqrt(-1); } return retval; }
/* wrapper y1 */ double y1 (double x) { if (__builtin_expect (x <= 0.0 || x > X_TLOSS, 0) && _LIB_VERSION != _IEEE_) { if (x < 0.0) { /* d = zero/(x-x) */ feraiseexcept (FE_INVALID); return __kernel_standard (x, x, 11); } else if (x == 0.0) /* d = -one/(x-x) */ return __kernel_standard (x, x, 10); else /* y1(x>X_TLOSS) */ return __kernel_standard (x, x, 37); } return __ieee754_y1 (x); }