Line data Source code
1 : #line 2 "../src/kernel/none/divll_pre.h"
2 : /* Copyright (C) 2014 The PARI group.
3 :
4 : This file is part of the PARI/GP package.
5 :
6 : PARI/GP is free software; you can redistribute it and/or modify it under the
7 : terms of the GNU General Public License as published by the Free Software
8 : Foundation; either version 2 of the License, or (at your option) any later
9 : version. It is distributed in the hope that it will be useful, but WITHOUT
10 : ANY WARRANTY WHATSOEVER.
11 :
12 : Check the License for details. You should have received a copy of it, along
13 : with the package; see the file 'COPYING'. If not, write to the Free Software
14 : Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. */
15 :
16 : #undef LOCAL_HIREMAINDER
17 : extern ulong hiremainder;
18 : #if defined(INLINE) && defined(__GNUC__) && !defined(DISABLE_INLINE)
19 : #define LOCAL_HIREMAINDER ulong hiremainder
20 : #else
21 : #define LOCAL_HIREMAINDER
22 : #endif
23 :
24 : #if defined(INLINE) && defined(__GNUC__) && !defined(DISABLE_INLINE)
25 : INLINE ulong /* precompute inverse of n */
26 1389473620 : get_Fl_red(ulong n)
27 : {
28 : LOCAL_HIREMAINDER;
29 1389473620 : n <<= bfffo(n);
30 1389473620 : hiremainder = ~n;
31 1389473620 : return divll(~0UL, n);
32 : }
33 : #else
34 : INLINE ulong /* precompute inverse of n */
35 414892434 : get_Fl_red(ulong n)
36 : {
37 414892434 : ulong q, oldhi = hiremainder;
38 414892434 : n <<= bfffo(n);
39 414892434 : hiremainder = ~n;
40 414892434 : q = divll(~0UL, n);
41 414892434 : hiremainder = oldhi;
42 414892434 : return q;
43 : }
44 : #endif
45 :
46 : INLINE ulong /* requires u1 <= n, n normalised */
47 10663680684 : divll_pre_normalized(ulong u1, ulong u0, ulong n, ulong ninv, ulong *pt_r)
48 : {
49 : ulong q0, q1, r;
50 : LOCAL_HIREMAINDER;
51 : LOCAL_OVERFLOW;
52 10663680684 : q0 = mulll(ninv, u1); q1 = hiremainder;
53 10663680684 : q0 = addll(q0, u0);
54 10663680684 : q1 = addllx(q1+1, u1);
55 10663680684 : r = u0 - q1 * n;
56 10663680684 : if (r > q0)
57 : {
58 7275160374 : r += n; q1--;
59 : }
60 10663680684 : if (r >= n)
61 : {
62 22098969 : r -= n; q1++;
63 : }
64 10663680684 : *pt_r = r; return q1;
65 : }
66 :
67 : INLINE ulong /* requires u1 <= n, n normalised */
68 14293376209 : remll_pre_normalized(ulong u1, ulong u0, ulong n, ulong ninv)
69 : {
70 : ulong q0, q1, r;
71 : LOCAL_HIREMAINDER;
72 : LOCAL_OVERFLOW;
73 14293376209 : q0 = mulll(ninv, u1); q1 = hiremainder;
74 14293376209 : q0 = addll(q0, u0);
75 14293376209 : q1 = addllx(q1, u1);
76 14293376209 : r = u0 - (q1 + 1) * n;
77 14293376209 : if (r >= q0)
78 10116492658 : r += n;
79 14293376209 : return r < n ? r : r - n;
80 : }
81 :
82 : INLINE ulong /* reduce <a_hi, a_lo> mod n */
83 14220843245 : remll_pre(ulong a_hi, ulong a_lo, ulong n, ulong ninv)
84 : {
85 14220843245 : int norm = bfffo(n);
86 14220843245 : int bits = BITS_IN_LONG - norm;
87 14220843245 : ulong sn = n << norm;
88 14220843245 : if (a_hi >= n) /* reduce a_hi first */
89 : {
90 71068273 : const ulong u1 = norm ? a_hi >> bits : 0;
91 71068273 : const ulong u0 = a_hi << norm;
92 71068273 : a_hi = remll_pre_normalized(u1, u0, sn, ninv) >> norm;
93 : }
94 : /* now reduce <a_hi, a_lo> */
95 : {
96 14220864649 : const ulong u1 = ((a_hi << norm) | (norm ? a_lo >> bits: 0));
97 14220864649 : const ulong u0 = a_lo << norm;
98 14220864649 : return remll_pre_normalized(u1, u0, sn, ninv) >> norm;
99 : }
100 : }
101 :
102 : #if !defined(INLINE)
103 : extern ulong divll_pre(ulong a_lo, ulong n, ulong ninv);
104 : #else
105 :
106 : #if defined(__GNUC__) && !defined(DISABLE_INLINE)
107 : #define divll_pre(a, n, ninv) \
108 : __extension__ ({ \
109 : ulong __a = (a); \
110 : ulong __n = (n); \
111 : int norm = bfffo(__n); \
112 : int bits = BITS_IN_LONG - norm; \
113 : ulong r, sn = __n << norm; \
114 : const ulong u1 = ((hiremainder << norm) | (norm ? __a >> bits: 0)); \
115 : const ulong u0 = __a << norm; \
116 : const ulong q = divll_pre_normalized(u1, u0, sn, ninv, &r); \
117 : hiremainder = r>>norm; q; \
118 : })
119 :
120 : #else /* __GNUC__ */
121 : INLINE ulong
122 2086262775 : divll_pre(ulong a_lo, ulong n, ulong ninv)
123 : {
124 2086262775 : int norm = bfffo(n);
125 2086262775 : int bits = BITS_IN_LONG - norm;
126 2086262775 : ulong r, sn = n << norm;
127 2086262775 : const ulong u1 = ((hiremainder << norm) | (norm ? a_lo >> bits: 0));
128 2086262775 : const ulong u0 = a_lo << norm;
129 2086262775 : const ulong q = divll_pre_normalized(u1, u0, sn, ninv, &r);
130 2086262775 : hiremainder = r>>norm; return q;
131 : }
132 : #endif /* __GNUC__ */
133 :
134 : #endif
|