1 /* Round to int floating-point values. PowerPC32 version.
2 Copyright (C) 2011-2016 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4 Contributed by Adhemerval Zanella <azanella@br.ibm.com>, 2011
6 The GNU C Library is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Lesser General Public
8 License as published by the Free Software Foundation; either
9 version 2.1 of the License, or (at your option) any later version.
11 The GNU C Library is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Lesser General Public License for more details.
16 You should have received a copy of the GNU Lesser General Public
17 License along with the GNU C Library; if not, see
18 <http://www.gnu.org/licenses/>. */
20 /* This has been coded in assembler because GCC makes such a mess of it
21 when it's coded in C. */
24 #include <math_ldbl_opt.h>
27 /* double [fp1] nearbyint(double [fp1] x) */
29 .section .rodata.cst4,"aM",@progbits,4
32 .long 0x59800000 /* TWO52: 2**52 */
41 addis r9,r9,.LC0-1b@ha
42 lfs fp13,.LC0-1b@l(r9)
50 fsub fp12,fp13,fp13 /* generate 0.0 */
51 fcmpu cr7,fp0,fp13 /* if (fabs(x) > TWO52 */
53 fcmpu cr7,fp1,fp12 /* if (x > 0.0 */
54 ble cr7,L(lessthanzero)
56 mtfsb0 4*cr7+lt /* Disable FE_INEXACT exception */
57 fadd fp1,fp1,fp13 /* x += TWO52 */
58 fsub fp1,fp1,fp13 /* x -= TWO52 */
59 fabs fp1,fp1 /* if (x == 0.0 */
60 mtfsf 0xff,fp11 /* Restore FE_INEXACT state. */
65 mtfsb0 4*cr7+lt /* Disable FE_INEXACT exception */
66 fsub fp1,fp1,fp13 /* x -= TWO52 */
67 fadd fp1,fp1,fp13 /* x += TWO52 */
68 fnabs fp1,fp1 /* if (x == 0.0) */
69 mtfsf 0xff,fp11 /* Restore FE_INEXACT state. */
72 /* Ensure sNaN input is converted to qNaN. */
79 weak_alias (__nearbyint, nearbyint)
82 weak_alias (__nearbyint, nearbyintl)
83 strong_alias (__nearbyint, __nearbyintl)
85 #if LONG_DOUBLE_COMPAT(libm, GLIBC_2_1)
86 compat_symbol (libm, __nearbyint, nearbyintl, GLIBC_2_1)