hjl / source-git / glibc

Forked from source-git/glibc 3 years ago
Clone

Blame sysdeps/powerpc/fpu/s_rintf.c

Packit 6c4009
/* Round a 32-bit floating point value to the nearest integer.
Packit 6c4009
   Copyright (C) 1997-2018 Free Software Foundation, Inc.
Packit 6c4009
   This file is part of the GNU C Library.
Packit 6c4009
Packit 6c4009
   The GNU C Library is free software; you can redistribute it and/or
Packit 6c4009
   modify it under the terms of the GNU Lesser General Public
Packit 6c4009
   License as published by the Free Software Foundation; either
Packit 6c4009
   version 2.1 of the License, or (at your option) any later version.
Packit 6c4009
Packit 6c4009
   The GNU C Library is distributed in the hope that it will be useful,
Packit 6c4009
   but WITHOUT ANY WARRANTY; without even the implied warranty of
Packit 6c4009
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
Packit 6c4009
   Lesser General Public License for more details.
Packit 6c4009
Packit 6c4009
   You should have received a copy of the GNU Lesser General Public
Packit 6c4009
   License along with the GNU C Library; if not, see
Packit 6c4009
   <http://www.gnu.org/licenses/>.  */
Packit 6c4009
Packit 6c4009
#include <math.h>
Packit 6c4009
#include <libm-alias-float.h>
Packit 6c4009
Packit 6c4009
float
Packit 6c4009
__rintf (float x)
Packit 6c4009
{
Packit 6c4009
  static const float TWO23 = 8388608.0;
Packit 6c4009
Packit 6c4009
  if (fabsf (x) < TWO23)
Packit 6c4009
    {
Packit 6c4009
      if (x > 0.0)
Packit 6c4009
	{
Packit 6c4009
	  x += TWO23;
Packit 6c4009
	  x -= TWO23;
Packit 6c4009
	}
Packit 6c4009
      else if (x < 0.0)
Packit 6c4009
	{
Packit 6c4009
	  x = TWO23 - x;
Packit 6c4009
	  x = -(x - TWO23);
Packit 6c4009
	}
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  return x;
Packit 6c4009
}
Packit 6c4009
libm_alias_float (__rint, rint)