Esempio n. 1
0
void
Math_y1(void *fp)
{
	F_Math_y1 *f;

	f = fp;

	*f->ret = __ieee754_y1(f->x);
}
Esempio n. 2
0
/*
 * @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;
}
Esempio n. 3
0
/* 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);
}