Mercurial > hg > CbC > CbC_gcc

annotate libquadmath/math/floorq.c @ 158:494b0b89df80 default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
...
author Shinji KONO <kono@ie.u-ryukyu.ac.jp>
date Mon, 25 May 2020 18:13:55 +0900
parents 1830386684a0
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
1 /* s_floorl.c -- long double version of s_floor.c.
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
2 * Conversion to IEEE quad long double by Jakub Jelinek, jj@ultra.linux.cz.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
3 */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
4
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
5 /*
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
6 * ====================================================
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
8 *
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
10 * Permission to use, copy, modify, and distribute this
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
11 * software is freely granted, provided that this notice
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
12 * is preserved.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
13 * ====================================================
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
14 */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
15
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
16 #if defined(LIBM_SCCS) && !defined(lint)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
17 static char rcsid[] = "$NetBSD: $";
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
18 #endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
19
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
20 /*
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
21 * floorq(x)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
22 * Return x rounded toward -inf to integral value
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
23 * Method:
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
24 * Bit twiddling.
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
25 */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
26
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
27 #define NO_MATH_REDIRECT
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
28
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
29 #include "quadmath-imp.h"
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
30
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
31 __float128 floorq(__float128 x)
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
32 {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
33 int64_t i0,i1,j0;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
34 uint64_t i,j;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
35 GET_FLT128_WORDS64(i0,i1,x);
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
36 j0 = ((i0>>48)&0x7fff)-0x3fff;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
37 if(j0<48) {
111
kono
parents: 68
diff changeset
38 if(j0<0) {
kono
parents: 68
diff changeset
39 /* return 0*sign(x) if |x|<1 */
kono
parents: 68
diff changeset
40 if(i0>=0) {i0=i1=0;}
kono
parents: 68
diff changeset
41 else if(((i0&0x7fffffffffffffffLL)|i1)!=0)
kono
parents: 68
diff changeset
42 { i0=0xbfff000000000000ULL;i1=0;}
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
43 } else {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
44 i = (0x0000ffffffffffffULL)>>j0;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
45 if(((i0&i)|i1)==0) return x; /* x is integral */
111
kono
parents: 68
diff changeset
46 if(i0<0) i0 += (0x0001000000000000LL)>>j0;
kono
parents: 68
diff changeset
47 i0 &= (~i); i1=0;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
48 }
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
49 } else if (j0>111) {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
50 if(j0==0x4000) return x+x; /* inf or NaN */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
51 else return x; /* x is integral */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
52 } else {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
53 i = -1ULL>>(j0-48);
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
54 if((i1&i)==0) return x; /* x is integral */
111
kono
parents: 68
diff changeset
55 if(i0<0) {
kono
parents: 68
diff changeset
56 if(j0==48) i0+=1;
kono
parents: 68
diff changeset
57 else {
kono
parents: 68
diff changeset
58 j = i1+(1LL<<(112-j0));
kono
parents: 68
diff changeset
59 if(j<i1) i0 +=1 ; /* got a carry */
kono
parents: 68
diff changeset
60 i1=j;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
61 }
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
62 }
111
kono
parents: 68
diff changeset
63 i1 &= (~i);
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
64 }
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
65 SET_FLT128_WORDS64(x,i0,i1);
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
66 return x;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
67 }