1  /* Copyright (C) 2013-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 Library General Public License as
       6     published by the Free Software Foundation; either version 2 of the
       7     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     Library General Public License for more details.
      13  
      14     You should have received a copy of the GNU Library General Public
      15     License along with the GNU C Library; see the file COPYING.LIB.  If
      16     not, see <https://www.gnu.org/licenses/>.  */
      17  
      18  /* ISA 2.07 provides fast GPR to FP instruction (mfvsr{d,wz}) which make
      19     generic implementation faster.  Also disables for old ISAs that do not
      20     have ceil/floor instructions.  */
      21  #if defined(_ARCH_PWR8) || !defined(_ARCH_PWR5X)
      22  # include <sysdeps/ieee754/flt-32/s_modff.c>
      23  #else
      24  # include <math.h>
      25  # include <libm-alias-float.h>
      26  
      27  float
      28  __modff (float x, float *iptr)
      29  {
      30    if (__builtin_isinff (x))
      31      {
      32        *iptr = x;
      33        return copysignf (0.0, x);
      34      }
      35    else if (__builtin_isnanf (x))
      36      {
      37        *iptr = NAN;
      38        return NAN;
      39      }
      40  
      41    if (x >= 0.0)
      42      {
      43        *iptr = floorf (x);
      44        return copysignf (x - *iptr, x);
      45      }
      46    else
      47      {
      48        *iptr = ceilf (x);
      49        return copysignf (x - *iptr, x);
      50      }
      51  }
      52  # ifndef __modff
      53  libm_alias_float (__modf, modf)
      54  # endif
      55  #endif