xref: /openbsd-src/lib/libcrypto/bn/bn_primitives.c (revision 99bf4692bf9b3bd9f10766e8554a4dfc4f9ebe9b)
1*99bf4692Sjsing /* $OpenBSD: bn_primitives.c,v 1.2 2023/06/21 07:48:41 jsing Exp $ */
2b5aa971dSjsing /*
3b5aa971dSjsing  * Copyright (c) 2023 Joel Sing <jsing@openbsd.org>
4b5aa971dSjsing  *
5b5aa971dSjsing  * Permission to use, copy, modify, and distribute this software for any
6b5aa971dSjsing  * purpose with or without fee is hereby granted, provided that the above
7b5aa971dSjsing  * copyright notice and this permission notice appear in all copies.
8b5aa971dSjsing  *
9b5aa971dSjsing  * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
10b5aa971dSjsing  * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
11b5aa971dSjsing  * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
12b5aa971dSjsing  * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
13b5aa971dSjsing  * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
14b5aa971dSjsing  * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
15b5aa971dSjsing  * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
16b5aa971dSjsing  */
17b5aa971dSjsing 
18b5aa971dSjsing #include <openssl/bn.h>
19b5aa971dSjsing 
20b5aa971dSjsing #include "bn_arch.h"
21b5aa971dSjsing #include "bn_internal.h"
22b5aa971dSjsing #include "bn_local.h"
23b5aa971dSjsing 
24*99bf4692Sjsing #ifndef HAVE_BN_CLZW
25b5aa971dSjsing #ifndef HAVE_BN_WORD_CLZ
26b5aa971dSjsing int
bn_word_clz(BN_ULONG w)27b5aa971dSjsing bn_word_clz(BN_ULONG w)
28b5aa971dSjsing {
29b5aa971dSjsing 	BN_ULONG bits, mask, shift;
30b5aa971dSjsing 
31b5aa971dSjsing 	bits = shift = BN_BITS2;
32b5aa971dSjsing 	mask = 0;
33b5aa971dSjsing 
34b5aa971dSjsing 	while ((shift >>= 1) != 0) {
35b5aa971dSjsing 		bits += (shift & mask) - (shift & ~mask);
36b5aa971dSjsing 		mask = bn_ct_ne_zero_mask(w >> bits);
37b5aa971dSjsing 	}
38b5aa971dSjsing 	bits += 1 & mask;
39b5aa971dSjsing 
40b5aa971dSjsing 	bits -= bn_ct_eq_zero(w);
41b5aa971dSjsing 
42b5aa971dSjsing 	return BN_BITS2 - bits;
43b5aa971dSjsing }
44b5aa971dSjsing #endif
45*99bf4692Sjsing #endif
46b5aa971dSjsing 
47b5aa971dSjsing #ifndef HAVE_BN_BITSIZE
48b5aa971dSjsing int
bn_bitsize(const BIGNUM * bn)49b5aa971dSjsing bn_bitsize(const BIGNUM *bn)
50b5aa971dSjsing {
51b5aa971dSjsing 	BN_ULONG n = 0, x = 0;
52b5aa971dSjsing 	BN_ULONG mask, w;
53b5aa971dSjsing 	int i = 0;
54b5aa971dSjsing 
55b5aa971dSjsing 	while (i < bn->top) {
56b5aa971dSjsing 		w = bn->d[i];
57b5aa971dSjsing 		mask = bn_ct_ne_zero_mask(w);
58b5aa971dSjsing 		n = ((BN_ULONG)i & mask) | (n & ~mask);
59b5aa971dSjsing 		x = (w & mask) | (x & ~mask);
60b5aa971dSjsing 		i++;
61b5aa971dSjsing 	}
62b5aa971dSjsing 
63*99bf4692Sjsing 	return (n + 1) * BN_BITS2 - bn_clzw(x);
64b5aa971dSjsing }
65b5aa971dSjsing #endif
66