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  #define NO_MATH_REDIRECT
      19  #define dsqrtl __hide_dsqrtl
      20  #define f32xsqrtf64 __hide_f32xsqrtf64
      21  #include <math.h>
      22  #undef dsqrtl
      23  #undef f32xsqrtf64
      24  #include <math_private.h>
      25  #include <math-svid-compat.h>
      26  #include <math-narrow-alias.h>
      27  #include <libm-alias-double.h>
      28  
      29  
      30  #if LIBM_SVID_COMPAT
      31  /* wrapper sqrt */
      32  double
      33  __sqrt (double x)
      34  {
      35    if (__builtin_expect (isless (x, 0.0), 0) && _LIB_VERSION != _IEEE_)
      36      return __kernel_standard (x, x, 26); /* sqrt(negative) */
      37  
      38    return __ieee754_sqrt (x);
      39  }
      40  libm_alias_double (__sqrt, sqrt)
      41  libm_alias_double_narrow (__sqrt, sqrt)
      42  #endif