sysdeps/powerpc/power5+/fpu/s_modf.c

   1 /* Copyright (C) 2013-2018 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 <http://www.gnu.org/licenses/>.  */
  17
  18 #include <math.h>
  19 #include <math_private.h>
  20 #include <math_ldbl_opt.h>
  21 #include <libm-alias-double.h>
  22
  23 double
  24 __modf (double x, double *iptr)
  25 {
  26   if (__builtin_isinf (x))
  27     {
  28       *iptr = x;
  29       return __copysign (0.0, x);
  30     }
  31   else if (__builtin_isnan (x))
  32     {
  33       *iptr = NAN;
  34       return NAN;
  35     }
  36
  37   if (x >= 0.0)
  38     {
  39       *iptr = __floor (x);
  40       return __copysign (x - *iptr, x);
  41     }
  42   else
  43     {
  44       *iptr = __ceil (x);
  45       return __copysign (x - *iptr, x);
  46     }
  47 }
  48 libm_alias_double (__modf, modf)
  49 #if LONG_DOUBLE_COMPAT (libc, GLIBC_2_0)
  50 compat_symbol (libc, __modf, modfl, GLIBC_2_0);
  51 #endif