1 /* Round to int floating-point values. RISC-V version. 2 Copyright (C) 2017-2021 Free Software Foundation, Inc. 3 This file is part of the GNU C Library. 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 16 License along with the GNU C Library. If not, see 17 <https://www.gnu.org/licenses/>. */ 18 19 #include <math.h> 20 #include <fenv_private.h> 21 #include <libm-alias-float.h> 22 #include <stdint.h> 23 24 float __nearbyintf(float x)25__nearbyintf (float x) 26 { 27 int flags = riscv_getflags (); 28 bool nan = isnan (x); 29 float mag = fabsf (x); 30 31 if (nan) 32 return x + x; 33 34 if (mag < (1 << __FLT_MANT_DIG__)) 35 { 36 int32_t i; 37 float new_x; 38 39 asm volatile ("fcvt.w.s %0, %1" : "=r" (i) : "f" (x)); 40 asm volatile ("fcvt.s.w %0, %1" : "=f" (new_x) : "r" (i)); 41 42 /* nearbyint(-0) == -0, and in general we'll always have the same 43 sign as our input. */ 44 x = copysignf (new_x, x); 45 46 riscv_setflags (flags); 47 } 48 49 return x; 50 } 51 52 libm_alias_float (__nearbyint, nearbyint) 53