Absolute File Name: | /home/opencoverage/opencoverage/guest-scripts/openssl/src/crypto/bn/bn_div.c |
Source code | Switch to Preprocessed file |
Line | Source | Count | ||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | /* | - | ||||||||||||||||||
2 | * Copyright 1995-2018 The OpenSSL Project Authors. All Rights Reserved. | - | ||||||||||||||||||
3 | * | - | ||||||||||||||||||
4 | * Licensed under the OpenSSL license (the "License"). You may not use | - | ||||||||||||||||||
5 | * this file except in compliance with the License. You can obtain a copy | - | ||||||||||||||||||
6 | * in the file LICENSE in the source distribution or at | - | ||||||||||||||||||
7 | * https://www.openssl.org/source/license.html | - | ||||||||||||||||||
8 | */ | - | ||||||||||||||||||
9 | - | |||||||||||||||||||
10 | #include <openssl/bn.h> | - | ||||||||||||||||||
11 | #include "internal/cryptlib.h" | - | ||||||||||||||||||
12 | #include "bn_lcl.h" | - | ||||||||||||||||||
13 | - | |||||||||||||||||||
14 | /* The old slow way */ | - | ||||||||||||||||||
15 | #if 0 | - | ||||||||||||||||||
16 | int BN_div(BIGNUM *dv, BIGNUM *rem, const BIGNUM *m, const BIGNUM *d, | - | ||||||||||||||||||
17 | BN_CTX *ctx) | - | ||||||||||||||||||
18 | { | - | ||||||||||||||||||
19 | int i, nm, nd; | - | ||||||||||||||||||
20 | int ret = 0; | - | ||||||||||||||||||
21 | BIGNUM *D; | - | ||||||||||||||||||
22 | - | |||||||||||||||||||
23 | bn_check_top(m); | - | ||||||||||||||||||
24 | bn_check_top(d); | - | ||||||||||||||||||
25 | if (BN_is_zero(d)) { | - | ||||||||||||||||||
26 | BNerr(BN_F_BN_DIV, BN_R_DIV_BY_ZERO); | - | ||||||||||||||||||
27 | return 0; | - | ||||||||||||||||||
28 | } | - | ||||||||||||||||||
29 | - | |||||||||||||||||||
30 | if (BN_ucmp(m, d) < 0) { | - | ||||||||||||||||||
31 | if (rem != NULL) { | - | ||||||||||||||||||
32 | if (BN_copy(rem, m) == NULL) | - | ||||||||||||||||||
33 | return 0; | - | ||||||||||||||||||
34 | } | - | ||||||||||||||||||
35 | if (dv != NULL) | - | ||||||||||||||||||
36 | BN_zero(dv); | - | ||||||||||||||||||
37 | return 1; | - | ||||||||||||||||||
38 | } | - | ||||||||||||||||||
39 | - | |||||||||||||||||||
40 | BN_CTX_start(ctx); | - | ||||||||||||||||||
41 | D = BN_CTX_get(ctx); | - | ||||||||||||||||||
42 | if (dv == NULL) | - | ||||||||||||||||||
43 | dv = BN_CTX_get(ctx); | - | ||||||||||||||||||
44 | if (rem == NULL) | - | ||||||||||||||||||
45 | rem = BN_CTX_get(ctx); | - | ||||||||||||||||||
46 | if (D == NULL || dv == NULL || rem == NULL) | - | ||||||||||||||||||
47 | goto end; | - | ||||||||||||||||||
48 | - | |||||||||||||||||||
49 | nd = BN_num_bits(d); | - | ||||||||||||||||||
50 | nm = BN_num_bits(m); | - | ||||||||||||||||||
51 | if (BN_copy(D, d) == NULL) | - | ||||||||||||||||||
52 | goto end; | - | ||||||||||||||||||
53 | if (BN_copy(rem, m) == NULL) | - | ||||||||||||||||||
54 | goto end; | - | ||||||||||||||||||
55 | - | |||||||||||||||||||
56 | /* | - | ||||||||||||||||||
57 | * The next 2 are needed so we can do a dv->d[0]|=1 later since | - | ||||||||||||||||||
58 | * BN_lshift1 will only work once there is a value :-) | - | ||||||||||||||||||
59 | */ | - | ||||||||||||||||||
60 | BN_zero(dv); | - | ||||||||||||||||||
61 | if (bn_wexpand(dv, 1) == NULL) | - | ||||||||||||||||||
62 | goto end; | - | ||||||||||||||||||
63 | dv->top = 1; | - | ||||||||||||||||||
64 | - | |||||||||||||||||||
65 | if (!BN_lshift(D, D, nm - nd)) | - | ||||||||||||||||||
66 | goto end; | - | ||||||||||||||||||
67 | for (i = nm - nd; i >= 0; i--) { | - | ||||||||||||||||||
68 | if (!BN_lshift1(dv, dv)) | - | ||||||||||||||||||
69 | goto end; | - | ||||||||||||||||||
70 | if (BN_ucmp(rem, D) >= 0) { | - | ||||||||||||||||||
71 | dv->d[0] |= 1; | - | ||||||||||||||||||
72 | if (!BN_usub(rem, rem, D)) | - | ||||||||||||||||||
73 | goto end; | - | ||||||||||||||||||
74 | } | - | ||||||||||||||||||
75 | /* CAN IMPROVE (and have now :=) */ | - | ||||||||||||||||||
76 | if (!BN_rshift1(D, D)) | - | ||||||||||||||||||
77 | goto end; | - | ||||||||||||||||||
78 | } | - | ||||||||||||||||||
79 | rem->neg = BN_is_zero(rem) ? 0 : m->neg; | - | ||||||||||||||||||
80 | dv->neg = m->neg ^ d->neg; | - | ||||||||||||||||||
81 | ret = 1; | - | ||||||||||||||||||
82 | end: | - | ||||||||||||||||||
83 | BN_CTX_end(ctx); | - | ||||||||||||||||||
84 | return ret; | - | ||||||||||||||||||
85 | } | - | ||||||||||||||||||
86 | - | |||||||||||||||||||
87 | #else | - | ||||||||||||||||||
88 | - | |||||||||||||||||||
89 | # if !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM) \ | - | ||||||||||||||||||
90 | && !defined(PEDANTIC) && !defined(BN_DIV3W) | - | ||||||||||||||||||
91 | # if defined(__GNUC__) && __GNUC__>=2 | - | ||||||||||||||||||
92 | # if defined(__i386) || defined (__i386__) | - | ||||||||||||||||||
93 | /*- | - | ||||||||||||||||||
94 | * There were two reasons for implementing this template: | - | ||||||||||||||||||
95 | * - GNU C generates a call to a function (__udivdi3 to be exact) | - | ||||||||||||||||||
96 | * in reply to ((((BN_ULLONG)n0)<<BN_BITS2)|n1)/d0 (I fail to | - | ||||||||||||||||||
97 | * understand why...); | - | ||||||||||||||||||
98 | * - divl doesn't only calculate quotient, but also leaves | - | ||||||||||||||||||
99 | * remainder in %edx which we can definitely use here:-) | - | ||||||||||||||||||
100 | */ | - | ||||||||||||||||||
101 | # undef bn_div_words | - | ||||||||||||||||||
102 | # define bn_div_words(n0,n1,d0) \ | - | ||||||||||||||||||
103 | ({ asm volatile ( \ | - | ||||||||||||||||||
104 | "divl %4" \ | - | ||||||||||||||||||
105 | : "=a"(q), "=d"(rem) \ | - | ||||||||||||||||||
106 | : "a"(n1), "d"(n0), "r"(d0) \ | - | ||||||||||||||||||
107 | : "cc"); \ | - | ||||||||||||||||||
108 | q; \ | - | ||||||||||||||||||
109 | }) | - | ||||||||||||||||||
110 | # define REMAINDER_IS_ALREADY_CALCULATED | - | ||||||||||||||||||
111 | # elif defined(__x86_64) && defined(SIXTY_FOUR_BIT_LONG) | - | ||||||||||||||||||
112 | /* | - | ||||||||||||||||||
113 | * Same story here, but it's 128-bit by 64-bit division. Wow! | - | ||||||||||||||||||
114 | */ | - | ||||||||||||||||||
115 | # undef bn_div_words | - | ||||||||||||||||||
116 | # define bn_div_words(n0,n1,d0) \ | - | ||||||||||||||||||
117 | ({ asm volatile ( \ | - | ||||||||||||||||||
118 | "divq %4" \ | - | ||||||||||||||||||
119 | : "=a"(q), "=d"(rem) \ | - | ||||||||||||||||||
120 | : "a"(n1), "d"(n0), "r"(d0) \ | - | ||||||||||||||||||
121 | : "cc"); \ | - | ||||||||||||||||||
122 | q; \ | - | ||||||||||||||||||
123 | }) | - | ||||||||||||||||||
124 | # define REMAINDER_IS_ALREADY_CALCULATED | - | ||||||||||||||||||
125 | # endif /* __<cpu> */ | - | ||||||||||||||||||
126 | # endif /* __GNUC__ */ | - | ||||||||||||||||||
127 | # endif /* OPENSSL_NO_ASM */ | - | ||||||||||||||||||
128 | - | |||||||||||||||||||
129 | /*- | - | ||||||||||||||||||
130 | * BN_div computes dv := num / divisor, rounding towards | - | ||||||||||||||||||
131 | * zero, and sets up rm such that dv*divisor + rm = num holds. | - | ||||||||||||||||||
132 | * Thus: | - | ||||||||||||||||||
133 | * dv->neg == num->neg ^ divisor->neg (unless the result is zero) | - | ||||||||||||||||||
134 | * rm->neg == num->neg (unless the remainder is zero) | - | ||||||||||||||||||
135 | * If 'dv' or 'rm' is NULL, the respective value is not returned. | - | ||||||||||||||||||
136 | */ | - | ||||||||||||||||||
137 | int BN_div(BIGNUM *dv, BIGNUM *rm, const BIGNUM *num, const BIGNUM *divisor, | - | ||||||||||||||||||
138 | BN_CTX *ctx) | - | ||||||||||||||||||
139 | { | - | ||||||||||||||||||
140 | int norm_shift, i, loop; | - | ||||||||||||||||||
141 | BIGNUM *tmp, wnum, *snum, *sdiv, *res; | - | ||||||||||||||||||
142 | BN_ULONG *resp, *wnump; | - | ||||||||||||||||||
143 | BN_ULONG d0, d1; | - | ||||||||||||||||||
144 | int num_n, div_n; | - | ||||||||||||||||||
145 | int no_branch = 0; | - | ||||||||||||||||||
146 | - | |||||||||||||||||||
147 | /* | - | ||||||||||||||||||
148 | * Invalid zero-padding would have particularly bad consequences so don't | - | ||||||||||||||||||
149 | * just rely on bn_check_top() here (bn_check_top() works only for | - | ||||||||||||||||||
150 | * BN_DEBUG builds) | - | ||||||||||||||||||
151 | */ | - | ||||||||||||||||||
152 | if ((num->top > 0 && num->d[num->top - 1] == 0) ||
| 0-7089159 | ||||||||||||||||||
153 | (divisor->top > 0 && divisor->d[divisor->top - 1] == 0)) {
| 0-7098495 | ||||||||||||||||||
154 | BNerr(BN_F_BN_DIV, BN_R_NOT_INITIALIZED); | - | ||||||||||||||||||
155 | return 0; never executed: return 0; | 0 | ||||||||||||||||||
156 | } | - | ||||||||||||||||||
157 | - | |||||||||||||||||||
158 | bn_check_top(num); | - | ||||||||||||||||||
159 | bn_check_top(divisor); | - | ||||||||||||||||||
160 | - | |||||||||||||||||||
161 | if ((BN_get_flags(num, BN_FLG_CONSTTIME) != 0)
| 2422072-4676426 | ||||||||||||||||||
162 | || (BN_get_flags(divisor, BN_FLG_CONSTTIME) != 0)) {
| 24686-4651740 | ||||||||||||||||||
163 | no_branch = 1; | - | ||||||||||||||||||
164 | } executed 2446758 times by 2 tests: end of block Executed by:
| 2446758 | ||||||||||||||||||
165 | - | |||||||||||||||||||
166 | bn_check_top(dv); | - | ||||||||||||||||||
167 | bn_check_top(rm); | - | ||||||||||||||||||
168 | /*- bn_check_top(num); *//* | - | ||||||||||||||||||
169 | * 'num' has been checked already | - | ||||||||||||||||||
170 | */ | - | ||||||||||||||||||
171 | /*- bn_check_top(divisor); *//* | - | ||||||||||||||||||
172 | * 'divisor' has been checked already | - | ||||||||||||||||||
173 | */ | - | ||||||||||||||||||
174 | - | |||||||||||||||||||
175 | if (BN_is_zero(divisor)) {
| 3-7098495 | ||||||||||||||||||
176 | BNerr(BN_F_BN_DIV, BN_R_DIV_BY_ZERO); | - | ||||||||||||||||||
177 | return 0; executed 3 times by 1 test: return 0; Executed by:
| 3 | ||||||||||||||||||
178 | } | - | ||||||||||||||||||
179 | - | |||||||||||||||||||
180 | if (!no_branch && BN_ucmp(num, divisor) < 0) {
| 230043-4651737 | ||||||||||||||||||
181 | if (rm != NULL) {
| 1-230042 | ||||||||||||||||||
182 | if (BN_copy(rm, num) == NULL)
| 0-230042 | ||||||||||||||||||
183 | return 0; never executed: return 0; | 0 | ||||||||||||||||||
184 | } executed 230042 times by 2 tests: end of block Executed by:
| 230042 | ||||||||||||||||||
185 | if (dv != NULL)
| 96-229947 | ||||||||||||||||||
186 | BN_zero(dv); executed 96 times by 1 test: (BN_set_word((dv),0)); Executed by:
| 96 | ||||||||||||||||||
187 | return 1; executed 230043 times by 2 tests: return 1; Executed by:
| 230043 | ||||||||||||||||||
188 | } | - | ||||||||||||||||||
189 | - | |||||||||||||||||||
190 | BN_CTX_start(ctx); | - | ||||||||||||||||||
191 | res = (dv == NULL) ? BN_CTX_get(ctx) : dv;
| 2507737-4360715 | ||||||||||||||||||
192 | tmp = BN_CTX_get(ctx); | - | ||||||||||||||||||
193 | snum = BN_CTX_get(ctx); | - | ||||||||||||||||||
194 | sdiv = BN_CTX_get(ctx); | - | ||||||||||||||||||
195 | if (sdiv == NULL)
| 0-6868452 | ||||||||||||||||||
196 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
197 | - | |||||||||||||||||||
198 | /* First we normalise the numbers */ | - | ||||||||||||||||||
199 | norm_shift = BN_BITS2 - ((BN_num_bits(divisor)) % BN_BITS2); | - | ||||||||||||||||||
200 | if (!(BN_lshift(sdiv, divisor, norm_shift)))
| 0-6868452 | ||||||||||||||||||
201 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
202 | sdiv->neg = 0; | - | ||||||||||||||||||
203 | norm_shift += BN_BITS2; | - | ||||||||||||||||||
204 | if (!(BN_lshift(snum, num, norm_shift)))
| 0-6868452 | ||||||||||||||||||
205 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
206 | snum->neg = 0; | - | ||||||||||||||||||
207 | - | |||||||||||||||||||
208 | if (no_branch) {
| 2446758-4421694 | ||||||||||||||||||
209 | /* | - | ||||||||||||||||||
210 | * Since we don't know whether snum is larger than sdiv, we pad snum | - | ||||||||||||||||||
211 | * with enough zeroes without changing its value. | - | ||||||||||||||||||
212 | */ | - | ||||||||||||||||||
213 | if (snum->top <= sdiv->top + 1) {
| 532746-1914012 | ||||||||||||||||||
214 | if (bn_wexpand(snum, sdiv->top + 2) == NULL)
| 0-532746 | ||||||||||||||||||
215 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
216 | for (i = snum->top; i < sdiv->top + 2; i++)
| 532746-532814 | ||||||||||||||||||
217 | snum->d[i] = 0; executed 532814 times by 1 test: snum->d[i] = 0; Executed by:
| 532814 | ||||||||||||||||||
218 | snum->top = sdiv->top + 2; | - | ||||||||||||||||||
219 | } else { executed 532746 times by 1 test: end of block Executed by:
| 532746 | ||||||||||||||||||
220 | if (bn_wexpand(snum, snum->top + 1) == NULL)
| 0-1914012 | ||||||||||||||||||
221 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
222 | snum->d[snum->top] = 0; | - | ||||||||||||||||||
223 | snum->top++; | - | ||||||||||||||||||
224 | } executed 1914012 times by 2 tests: end of block Executed by:
| 1914012 | ||||||||||||||||||
225 | } | - | ||||||||||||||||||
226 | - | |||||||||||||||||||
227 | div_n = sdiv->top; | - | ||||||||||||||||||
228 | num_n = snum->top; | - | ||||||||||||||||||
229 | loop = num_n - div_n; | - | ||||||||||||||||||
230 | /* | - | ||||||||||||||||||
231 | * Lets setup a 'window' into snum This is the part that corresponds to | - | ||||||||||||||||||
232 | * the current 'area' being divided | - | ||||||||||||||||||
233 | */ | - | ||||||||||||||||||
234 | wnum.neg = 0; | - | ||||||||||||||||||
235 | wnum.d = &(snum->d[loop]); | - | ||||||||||||||||||
236 | wnum.top = div_n; | - | ||||||||||||||||||
237 | wnum.flags = BN_FLG_STATIC_DATA; | - | ||||||||||||||||||
238 | /* | - | ||||||||||||||||||
239 | * only needed when BN_ucmp messes up the values between top and max | - | ||||||||||||||||||
240 | */ | - | ||||||||||||||||||
241 | wnum.dmax = snum->dmax - loop; /* so we don't step out of bounds */ | - | ||||||||||||||||||
242 | - | |||||||||||||||||||
243 | /* Get the top 2 words of sdiv */ | - | ||||||||||||||||||
244 | /* div_n=sdiv->top; */ | - | ||||||||||||||||||
245 | d0 = sdiv->d[div_n - 1]; | - | ||||||||||||||||||
246 | d1 = (div_n == 1) ? 0 : sdiv->d[div_n - 2];
| 534748-6333704 | ||||||||||||||||||
247 | - | |||||||||||||||||||
248 | /* pointer to the 'top' of snum */ | - | ||||||||||||||||||
249 | wnump = &(snum->d[num_n - 1]); | - | ||||||||||||||||||
250 | - | |||||||||||||||||||
251 | /* Setup to 'res' */ | - | ||||||||||||||||||
252 | if (!bn_wexpand(res, (loop + 1)))
| 0-6868452 | ||||||||||||||||||
253 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
254 | res->neg = (num->neg ^ divisor->neg); | - | ||||||||||||||||||
255 | res->top = loop - no_branch; | - | ||||||||||||||||||
256 | resp = &(res->d[loop - 1]); | - | ||||||||||||||||||
257 | - | |||||||||||||||||||
258 | /* space for temp */ | - | ||||||||||||||||||
259 | if (!bn_wexpand(tmp, (div_n + 1)))
| 0-6868452 | ||||||||||||||||||
260 | goto err; never executed: goto err; | 0 | ||||||||||||||||||
261 | - | |||||||||||||||||||
262 | if (!no_branch) {
| 2446758-4421694 | ||||||||||||||||||
263 | if (BN_ucmp(&wnum, sdiv) >= 0) {
| 32266-4389428 | ||||||||||||||||||
264 | /* | - | ||||||||||||||||||
265 | * If BN_DEBUG_RAND is defined BN_ucmp changes (via bn_pollute) | - | ||||||||||||||||||
266 | * the const bignum arguments => clean the values between top and | - | ||||||||||||||||||
267 | * max again | - | ||||||||||||||||||
268 | */ | - | ||||||||||||||||||
269 | bn_clear_top2max(&wnum); | - | ||||||||||||||||||
270 | bn_sub_words(wnum.d, wnum.d, sdiv->d, div_n); | - | ||||||||||||||||||
271 | *resp = 1; | - | ||||||||||||||||||
272 | } else executed 32266 times by 2 tests: end of block Executed by:
| 32266 | ||||||||||||||||||
273 | res->top--; executed 4389428 times by 2 tests: res->top--; Executed by:
| 4389428 | ||||||||||||||||||
274 | } | - | ||||||||||||||||||
275 | - | |||||||||||||||||||
276 | /* Increase the resp pointer so that we never create an invalid pointer. */ | - | ||||||||||||||||||
277 | resp++; | - | ||||||||||||||||||
278 | - | |||||||||||||||||||
279 | /* | - | ||||||||||||||||||
280 | * if res->top == 0 then clear the neg value otherwise decrease the resp | - | ||||||||||||||||||
281 | * pointer | - | ||||||||||||||||||
282 | */ | - | ||||||||||||||||||
283 | if (res->top == 0)
| 0-6868452 | ||||||||||||||||||
284 | res->neg = 0; never executed: res->neg = 0; | 0 | ||||||||||||||||||
285 | else | - | ||||||||||||||||||
286 | resp--; executed 6868452 times by 2 tests: resp--; Executed by:
| 6868452 | ||||||||||||||||||
287 | - | |||||||||||||||||||
288 | for (i = 0; i < loop - 1; i++, wnump--) {
| 6868452-22741823 | ||||||||||||||||||
289 | BN_ULONG q, l0; | - | ||||||||||||||||||
290 | /* | - | ||||||||||||||||||
291 | * the first part of the loop uses the top two words of snum and sdiv | - | ||||||||||||||||||
292 | * to calculate a BN_ULONG q such that | wnum - sdiv * q | < sdiv | - | ||||||||||||||||||
293 | */ | - | ||||||||||||||||||
294 | # if defined(BN_DIV3W) && !defined(OPENSSL_NO_ASM) | - | ||||||||||||||||||
295 | BN_ULONG bn_div_3_words(BN_ULONG *, BN_ULONG, BN_ULONG); | - | ||||||||||||||||||
296 | q = bn_div_3_words(wnump, d1, d0); | - | ||||||||||||||||||
297 | # else | - | ||||||||||||||||||
298 | BN_ULONG n0, n1, rem = 0; | - | ||||||||||||||||||
299 | - | |||||||||||||||||||
300 | n0 = wnump[0]; | - | ||||||||||||||||||
301 | n1 = wnump[-1]; | - | ||||||||||||||||||
302 | if (n0 == d0)
| 90174-22651649 | ||||||||||||||||||
303 | q = BN_MASK2; executed 90174 times by 1 test: q = (0xffffffffffffffffL); Executed by:
| 90174 | ||||||||||||||||||
304 | else { /* n0 < d0 */ | - | ||||||||||||||||||
305 | - | |||||||||||||||||||
306 | # ifdef BN_LLONG | - | ||||||||||||||||||
307 | BN_ULLONG t2; | - | ||||||||||||||||||
308 | - | |||||||||||||||||||
309 | # if defined(BN_LLONG) && defined(BN_DIV2W) && !defined(bn_div_words) | - | ||||||||||||||||||
310 | q = (BN_ULONG)(((((BN_ULLONG) n0) << BN_BITS2) | n1) / d0); | - | ||||||||||||||||||
311 | # else | - | ||||||||||||||||||
312 | q = bn_div_words(n0, n1, d0); | - | ||||||||||||||||||
313 | # endif | - | ||||||||||||||||||
314 | - | |||||||||||||||||||
315 | # ifndef REMAINDER_IS_ALREADY_CALCULATED | - | ||||||||||||||||||
316 | /* | - | ||||||||||||||||||
317 | * rem doesn't have to be BN_ULLONG. The least we | - | ||||||||||||||||||
318 | * know it's less that d0, isn't it? | - | ||||||||||||||||||
319 | */ | - | ||||||||||||||||||
320 | rem = (n1 - q * d0) & BN_MASK2; | - | ||||||||||||||||||
321 | # endif | - | ||||||||||||||||||
322 | t2 = (BN_ULLONG) d1 *q; | - | ||||||||||||||||||
323 | - | |||||||||||||||||||
324 | for (;;) { | - | ||||||||||||||||||
325 | if (t2 <= ((((BN_ULLONG) rem) << BN_BITS2) | wnump[-2])) | - | ||||||||||||||||||
326 | break; | - | ||||||||||||||||||
327 | q--; | - | ||||||||||||||||||
328 | rem += d0; | - | ||||||||||||||||||
329 | if (rem < d0) | - | ||||||||||||||||||
330 | break; /* don't let rem overflow */ | - | ||||||||||||||||||
331 | t2 -= d1; | - | ||||||||||||||||||
332 | } | - | ||||||||||||||||||
333 | # else /* !BN_LLONG */ | - | ||||||||||||||||||
334 | BN_ULONG t2l, t2h; | - | ||||||||||||||||||
335 | - | |||||||||||||||||||
336 | q = bn_div_words(n0, n1, d0); | - | ||||||||||||||||||
337 | # ifndef REMAINDER_IS_ALREADY_CALCULATED | - | ||||||||||||||||||
338 | rem = (n1 - q * d0) & BN_MASK2; | - | ||||||||||||||||||
339 | # endif | - | ||||||||||||||||||
340 | - | |||||||||||||||||||
341 | # if defined(BN_UMULT_LOHI) | - | ||||||||||||||||||
342 | BN_UMULT_LOHI(t2l, t2h, d1, q); | - | ||||||||||||||||||
343 | # elif defined(BN_UMULT_HIGH) | - | ||||||||||||||||||
344 | t2l = d1 * q; | - | ||||||||||||||||||
345 | t2h = BN_UMULT_HIGH(d1, q); | - | ||||||||||||||||||
346 | # else | - | ||||||||||||||||||
347 | { | - | ||||||||||||||||||
348 | BN_ULONG ql, qh; | - | ||||||||||||||||||
349 | t2l = LBITS(d1); | - | ||||||||||||||||||
350 | t2h = HBITS(d1); | - | ||||||||||||||||||
351 | ql = LBITS(q); | - | ||||||||||||||||||
352 | qh = HBITS(q); | - | ||||||||||||||||||
353 | mul64(t2l, t2h, ql, qh); /* t2=(BN_ULLONG)d1*q; */ | - | ||||||||||||||||||
354 | } | - | ||||||||||||||||||
355 | # endif | - | ||||||||||||||||||
356 | - | |||||||||||||||||||
357 | for (;;) { | - | ||||||||||||||||||
358 | if ((t2h < rem) || ((t2h == rem) && (t2l <= wnump[-2])))
| 78657-16235308 | ||||||||||||||||||
359 | break; executed 16568745 times by 2 tests: break; Executed by:
| 16568745 | ||||||||||||||||||
360 | q--; | - | ||||||||||||||||||
361 | rem += d0; | - | ||||||||||||||||||
362 | if (rem < d0)
| 470532-6082904 | ||||||||||||||||||
363 | break; /* don't let rem overflow */ executed 6082904 times by 2 tests: break; Executed by:
| 6082904 | ||||||||||||||||||
364 | if (t2l < d1)
| 168517-302015 | ||||||||||||||||||
365 | t2h--; executed 302015 times by 2 tests: t2h--; Executed by:
| 302015 | ||||||||||||||||||
366 | t2l -= d1; | - | ||||||||||||||||||
367 | } executed 470532 times by 2 tests: end of block Executed by:
| 470532 | ||||||||||||||||||
368 | # endif /* !BN_LLONG */ | - | ||||||||||||||||||
369 | } executed 22651649 times by 2 tests: end of block Executed by:
| 22651649 | ||||||||||||||||||
370 | # endif /* !BN_DIV3W */ | - | ||||||||||||||||||
371 | - | |||||||||||||||||||
372 | l0 = bn_mul_words(tmp->d, sdiv->d, div_n, q); | - | ||||||||||||||||||
373 | tmp->d[div_n] = l0; | - | ||||||||||||||||||
374 | wnum.d--; | - | ||||||||||||||||||
375 | /* | - | ||||||||||||||||||
376 | * ingore top values of the bignums just sub the two BN_ULONG arrays | - | ||||||||||||||||||
377 | * with bn_sub_words | - | ||||||||||||||||||
378 | */ | - | ||||||||||||||||||
379 | if (bn_sub_words(wnum.d, wnum.d, tmp->d, div_n + 1)) {
| 78007-22663816 | ||||||||||||||||||
380 | /* | - | ||||||||||||||||||
381 | * Note: As we have considered only the leading two BN_ULONGs in | - | ||||||||||||||||||
382 | * the calculation of q, sdiv * q might be greater than wnum (but | - | ||||||||||||||||||
383 | * then (q-1) * sdiv is less or equal than wnum) | - | ||||||||||||||||||
384 | */ | - | ||||||||||||||||||
385 | q--; | - | ||||||||||||||||||
386 | if (bn_add_words(wnum.d, wnum.d, sdiv->d, div_n))
| 0-78007 | ||||||||||||||||||
387 | /* | - | ||||||||||||||||||
388 | * we can't have an overflow here (assuming that q != 0, but | - | ||||||||||||||||||
389 | * if q == 0 then tmp is zero anyway) | - | ||||||||||||||||||
390 | */ | - | ||||||||||||||||||
391 | (*wnump)++; executed 78007 times by 1 test: (*wnump)++; Executed by:
| 78007 | ||||||||||||||||||
392 | } executed 78007 times by 1 test: end of block Executed by:
| 78007 | ||||||||||||||||||
393 | /* store part of the result */ | - | ||||||||||||||||||
394 | resp--; | - | ||||||||||||||||||
395 | *resp = q; | - | ||||||||||||||||||
396 | } executed 22741823 times by 2 tests: end of block Executed by:
| 22741823 | ||||||||||||||||||
397 | bn_correct_top(snum); | - | ||||||||||||||||||
398 | if (rm != NULL) {
| 84287-6784165 | ||||||||||||||||||
399 | /* | - | ||||||||||||||||||
400 | * Keep a copy of the neg flag in num because if rm==num BN_rshift() | - | ||||||||||||||||||
401 | * will overwrite it. | - | ||||||||||||||||||
402 | */ | - | ||||||||||||||||||
403 | int neg = num->neg; | - | ||||||||||||||||||
404 | BN_rshift(rm, snum, norm_shift); | - | ||||||||||||||||||
405 | if (!BN_is_zero(rm))
| 32906-6751259 | ||||||||||||||||||
406 | rm->neg = neg; executed 6751259 times by 2 tests: rm->neg = neg; Executed by:
| 6751259 | ||||||||||||||||||
407 | bn_check_top(rm); | - | ||||||||||||||||||
408 | } executed 6784165 times by 2 tests: end of block Executed by:
| 6784165 | ||||||||||||||||||
409 | if (no_branch)
| 2446758-4421694 | ||||||||||||||||||
410 | bn_correct_top(res); executed 2446758 times by 2 tests: bn_correct_top(res); Executed by:
| 2446758 | ||||||||||||||||||
411 | BN_CTX_end(ctx); | - | ||||||||||||||||||
412 | return 1; executed 6868452 times by 2 tests: return 1; Executed by:
| 6868452 | ||||||||||||||||||
413 | err: | - | ||||||||||||||||||
414 | bn_check_top(rm); | - | ||||||||||||||||||
415 | BN_CTX_end(ctx); | - | ||||||||||||||||||
416 | return 0; never executed: return 0; | 0 | ||||||||||||||||||
417 | } | - | ||||||||||||||||||
418 | #endif | - | ||||||||||||||||||
Source code | Switch to Preprocessed file |