projects / thirdparty / glibc.git / blame
| Commit | Line | Data |
|---|---|---|
| d4697bc9 | 1 | /* Copyright (C) 2007-2014 Free Software Foundation, Inc. |
| f878f0bb RH |
2 | This file is part of the GNU C Library. |
| 3 | Contributed by Richard Henderson. | |
| 4 | ||
| 5 | The GNU C Library is free software; you can redistribute it and/or | |
| 6 | modify it under the terms of the GNU Lesser General Public | |
| 7 | License as published by the Free Software Foundation; either | |
| 8 | version 2.1 of the License, or (at your option) any later version. | |
| 9 | ||
| 10 | The GNU C Library is distributed in the hope that it will be useful, | |
| 11 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 12 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
| 13 | Lesser General Public License for more details. | |
| 14 | ||
| 15 | You should have received a copy of the GNU Lesser General Public | |
| ab84e3ff PE |
16 | License along with the GNU C Library. If not, see |
| 17 | <http://www.gnu.org/licenses/>. */ | |
| f878f0bb RH |
18 | |
| 19 | #include <math.h> | |
| 20 | ||
| 21 | ||
| 22 | float | |
| 23 | __roundf (float x) | |
| 24 | { | |
| 25 | const float almost_half = 0x1.fffffep-2; | |
| 26 | const float two23 = 0x1.0p23; | |
| 27 | float r, tmp; | |
| 5556231d | 28 | |
| f878f0bb RH |
29 | __asm ( |
| 30 | #ifdef _IEEE_FP_INEXACT | |
| 31 | "adds/suic %2, %3, %1\n\tsubs/suic %1, %3, %0" | |
| 32 | #else | |
| 33 | "adds/suc %2, %3, %1\n\tsubs/suc %1, %3, %0" | |
| 34 | #endif | |
| 35 | : "=&f"(r), "=&f"(tmp) | |
| 36 | : "f"(fabsf (x) + almost_half), "f"(two23)); | |
| 37 | ||
| 38 | /* round(-0) == -0, and in general we'll always have the same | |
| 39 | sign as our input. */ | |
| 40 | return copysignf (r, x); | |
| 41 | } | |
| 42 | ||
| 43 | weak_alias (__roundf, roundf) |