Mercurial > hg > CbC > CbC_gcc

annotate libquadmath/math/remainderq.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 /* e_fmodl.c -- long double version of e_fmod.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 * 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
7 *
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
8 * 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
9 * Permission to use, copy, modify, and distribute this
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
10 * software is freely granted, provided that this notice
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
11 * is preserved.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
12 * ====================================================
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 /* remainderq(x,p)
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
16 * Return :
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
17 * returns x REM p = x - [x/p]*p as if in infinite
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
18 * precise arithmetic, where [x/p] is the (infinite bit)
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
19 * integer nearest x/p (in half way case choose the even one).
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
20 * Method :
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
21 * Based on fmodl() return x-[x/p]chopped*p exactlp.
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
22 */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
23
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
24 #include "quadmath-imp.h"
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
25
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
26 static const __float128 zero = 0;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
27
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
28
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
29 __float128
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
30 remainderq(__float128 x, __float128 p)
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
31 {
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
32 int64_t hx,hp;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
33 uint64_t sx,lx,lp;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
34 __float128 p_half;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
35
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
36 GET_FLT128_WORDS64(hx,lx,x);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
37 GET_FLT128_WORDS64(hp,lp,p);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
38 sx = hx&0x8000000000000000ULL;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
39 hp &= 0x7fffffffffffffffLL;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
40 hx &= 0x7fffffffffffffffLL;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
41
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
42 /* purge off exception values */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
43 if((hp|lp)==0) return (x*p)/(x*p); /* p = 0 */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
44 if((hx>=0x7fff000000000000LL)|| /* x not finite */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
45 ((hp>=0x7fff000000000000LL)&& /* p is NaN */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
46 (((hp-0x7fff000000000000LL)|lp)!=0)))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
47 return (x*p)/(x*p);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
48
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
49
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
50 if (hp<=0x7ffdffffffffffffLL) x = fmodq(x,p+p); /* now x < 2p */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
51 if (((hx-hp)|(lx-lp))==0) return zero*x;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
52 x = fabsq(x);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
53 p = fabsq(p);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
54 if (hp<0x0002000000000000LL) {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
55 if(x+x>p) {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
56 x-=p;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
57 if(x+x>=p) x -= p;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
58 }
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
59 } else {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
60 p_half = 0.5Q*p;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
61 if(x>p_half) {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
62 x-=p;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
63 if(x>=p_half) x -= p;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
64 }
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
65 }
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
66 GET_FLT128_MSW64(hx,x);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
67 SET_FLT128_MSW64(x,hx^sx);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
68 return x;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
69 }