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

   1 /* Copyright (C) 2013-2014 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
  22 double
  23 __modf (double x, double *iptr)
  24 {
  25   if (__builtin_isinf (x))
  26     {
  27       *iptr = x;
  28       return __copysign (0.0, x);
  29     }
  30   else if (__builtin_isnan (x))
  31     {
  32       *iptr = NAN;
  33       return NAN;
  34     }
  35
  36   if (x >= 0.0)
  37     {
  38       *iptr = __floor (x);
  39       return __copysign (x - *iptr, x);
  40     }
  41   else
  42     {
  43       *iptr = __ceil (x);
  44       return __copysign (x - *iptr, x);
  45     }
  46 }
  47 weak_alias (__modf, modf)
  48 #ifdef NO_LONG_DOUBLE
  49 strong_alias (__modf, __modfl)
  50 weak_alias (__modf, modfl)
  51 #endif
  52 #if IS_IN (libm)
  53 # if LONG_DOUBLE_COMPAT(libm, GLIBC_2_0)
  54 compat_symbol (libm, __modf, modfl, GLIBC_2_0);
  55 # endif
  56 #elif LONG_DOUBLE_COMPAT(libc, GLIBC_2_0)
  57 compat_symbol (libc, __modf, modfl, GLIBC_2_0);
  58 #endif