1  /* Copyright (C) 2011-2023 Free Software Foundation, Inc.
       2     This file is part of the GNU C Library.
       3  
       4     The GNU C Library is free software; you can redistribute it and/or
       5     modify it under the terms of the GNU Lesser General Public
       6     License as published by the Free Software Foundation; either
       7     version 2.1 of the License, or (at your option) any later version.
       8  
       9     The GNU C Library is distributed in the hope that it will be useful,
      10     but WITHOUT ANY WARRANTY; without even the implied warranty of
      11     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
      12     Lesser General Public License for more details.
      13  
      14     You should have received a copy of the GNU Lesser General Public
      15     License along with the GNU C Library; if not, see
      16     <https://www.gnu.org/licenses/>.  */
      17  
      18  #include <fenv.h>
      19  #include <math.h>
      20  #include <math_private.h>
      21  #include <math-svid-compat.h>
      22  #include <libm-alias-float.h>
      23  
      24  
      25  #if LIBM_SVID_COMPAT
      26  /* wrapper j0f */
      27  float
      28  __j0f (float x)
      29  {
      30    if (__builtin_expect (isgreater (fabsf (x),
      31  				   AS_FLOAT_CONSTANT (X_TLOSS)), 0)
      32        && _LIB_VERSION != _IEEE_ && _LIB_VERSION != _POSIX_)
      33      /* j0(|x|>X_TLOSS) */
      34      return __kernel_standard_f (x, x, 134);
      35  
      36    return __ieee754_j0f (x);
      37  }
      38  libm_alias_float (__j0, j0)
      39  
      40  
      41  /* wrapper y0f */
      42  float
      43  __y0f (float x)
      44  {
      45    if (__builtin_expect (islessequal (x, 0.0f)
      46                          || isgreater (x, AS_FLOAT_CONSTANT (X_TLOSS)), 0)
      47        && _LIB_VERSION != _IEEE_)
      48      {
      49        if (x < 0.0f)
      50  	{
      51  	  /* d = zero/(x-x) */
      52  	  feraiseexcept (FE_INVALID);
      53  	  return __kernel_standard_f (x, x, 109);
      54  	}
      55        else if (x == 0.0f)
      56  	{
      57  	  /* d = -one/(x-x) */
      58  	  feraiseexcept (FE_DIVBYZERO);
      59  	  return __kernel_standard_f (x, x, 108);
      60  	}
      61        else if (_LIB_VERSION != _POSIX_)
      62  	/* y0(x>X_TLOSS) */
      63  	return __kernel_standard_f (x, x, 135);
      64      }
      65  
      66    return __ieee754_y0f (x);
      67  }
      68  libm_alias_float (__y0, y0)
      69  #endif