Back to index

glibc  2.9
s_rintf.c
Go to the documentation of this file.
00001 /* Round a 32-bit floating point value to the nearest integer.
00002    Copyright (C) 1997, 2000 Free Software Foundation, Inc.
00003    This file is part of the GNU C Library.
00004 
00005    The GNU C Library is free software; you can redistribute it and/or
00006    modify it under the terms of the GNU Lesser General Public
00007    License as published by the Free Software Foundation; either
00008    version 2.1 of the License, or (at your option) any later version.
00009 
00010    The GNU C Library is distributed in the hope that it will be useful,
00011    but WITHOUT ANY WARRANTY; without even the implied warranty of
00012    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00013    Lesser General Public License for more details.
00014 
00015    You should have received a copy of the GNU Lesser General Public
00016    License along with the GNU C Library; if not, write to the Free
00017    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
00018    02111-1307 USA.  */
00019 
00020 #include "math.h"
00021 
00022 float
00023 __rintf (float x)
00024 {
00025   static const float TWO23 = 8388608.0;
00026 
00027   if (fabsf (x) < TWO23)
00028     {
00029       if (x > 0.0)
00030        {
00031          x += TWO23;
00032          x -= TWO23;
00033        }
00034       else if (x < 0.0)
00035        {
00036          x = TWO23 - x;
00037          x = -(x - TWO23);
00038        }
00039     }
00040   
00041   return x;
00042 }
00043 weak_alias (__rintf, rintf)