]>
git.ipfire.org Git - thirdparty/glibc.git/blob - sysdeps/powerpc/powerpc32/fpu/s_llrint.c
1 /* Round a double value to a long long in the current rounding mode.
2 Copyright (C) 1997-2021 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
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.
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.
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/>. */
21 #include <math_ldbl_opt.h>
22 #include <math_private.h>
24 #include <libm-alias-double.h>
30 /* Assume powerpc64 instructions availability. */
32 __asm__ ("fctid %0, %1" : "=d" (ret
) : "d" (x
));
36 if (HAVE_PPC_FCTIDZ
|| rx
!= x
)
37 return (long long int) rx
;
40 /* Avoid incorrect exceptions from libgcc conversions (as of GCC
41 5): <https://gcc.gnu.org/bugzilla/show_bug.cgi?id=59412>. */
42 if (fabs (rx
) < 0x1p
31)
43 return (long long int) (long int) rx
;
45 EXTRACT_WORDS64 (i0
, rx
);
46 int exponent
= ((i0
>> 52) & 0x7ff) - 0x3ff;
49 unsigned long long int mant
50 = (i0
& ((1ULL << 52) - 1)) | (1ULL << 52);
52 mant
>>= 52 - exponent
;
54 mant
<<= exponent
- 52;
55 return (long long int) ((i0
& (1ULL << 63)) != 0 ? -mant
: mant
);
57 else if (rx
== (double) LLONG_MIN
)
60 return (long long int) (long int) rx
<< 32;
65 libm_alias_double (__llrint
, llrint
)