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