1*b077aed3SPierre Pronchery /*
2*b077aed3SPierre Pronchery * Copyright 2019-2022 The OpenSSL Project Authors. All Rights Reserved.
3*b077aed3SPierre Pronchery * Copyright (c) 2019, Oracle and/or its affiliates. All rights reserved.
4*b077aed3SPierre Pronchery *
5*b077aed3SPierre Pronchery * Licensed under the Apache License 2.0 (the "License"). You may not use
6*b077aed3SPierre Pronchery * this file except in compliance with the License. You can obtain a copy
7*b077aed3SPierre Pronchery * in the file LICENSE in the source distribution or at
8*b077aed3SPierre Pronchery * https://www.openssl.org/source/license.html
9*b077aed3SPierre Pronchery */
10*b077aed3SPierre Pronchery
11*b077aed3SPierre Pronchery #include <openssl/crypto.h>
12*b077aed3SPierre Pronchery #include <openssl/bn.h>
13*b077aed3SPierre Pronchery #include "crypto/sparse_array.h"
14*b077aed3SPierre Pronchery
15*b077aed3SPierre Pronchery /*
16*b077aed3SPierre Pronchery * How many bits are used to index each level in the tree structure?
17*b077aed3SPierre Pronchery * This setting determines the number of pointers stored in each node of the
18*b077aed3SPierre Pronchery * tree used to represent the sparse array. Having more pointers reduces the
19*b077aed3SPierre Pronchery * depth of the tree but potentially wastes more memory. That is, this is a
20*b077aed3SPierre Pronchery * direct space versus time tradeoff.
21*b077aed3SPierre Pronchery *
22*b077aed3SPierre Pronchery * The default is to use four bits which means that the are 16
23*b077aed3SPierre Pronchery * pointers in each tree node.
24*b077aed3SPierre Pronchery *
25*b077aed3SPierre Pronchery * The library builder is also permitted to define other sizes in the closed
26*b077aed3SPierre Pronchery * interval [2, sizeof(ossl_uintmax_t) * 8]. Space use generally scales
27*b077aed3SPierre Pronchery * exponentially with the block size, although the implementation only
28*b077aed3SPierre Pronchery * creates enough blocks to support the largest used index. The depth is:
29*b077aed3SPierre Pronchery * ceil(log_2(largest index) / 2^{block size})
30*b077aed3SPierre Pronchery * E.g. with a block size of 4, and a largest index of 1000, the depth
31*b077aed3SPierre Pronchery * will be three.
32*b077aed3SPierre Pronchery */
33*b077aed3SPierre Pronchery #ifndef OPENSSL_SA_BLOCK_BITS
34*b077aed3SPierre Pronchery # define OPENSSL_SA_BLOCK_BITS 4
35*b077aed3SPierre Pronchery #elif OPENSSL_SA_BLOCK_BITS < 2 || OPENSSL_SA_BLOCK_BITS > (BN_BITS2 - 1)
36*b077aed3SPierre Pronchery # error OPENSSL_SA_BLOCK_BITS is out of range
37*b077aed3SPierre Pronchery #endif
38*b077aed3SPierre Pronchery
39*b077aed3SPierre Pronchery /*
40*b077aed3SPierre Pronchery * From the number of bits, work out:
41*b077aed3SPierre Pronchery * the number of pointers in a tree node;
42*b077aed3SPierre Pronchery * a bit mask to quickly extract an index and
43*b077aed3SPierre Pronchery * the maximum depth of the tree structure.
44*b077aed3SPierre Pronchery */
45*b077aed3SPierre Pronchery #define SA_BLOCK_MAX (1 << OPENSSL_SA_BLOCK_BITS)
46*b077aed3SPierre Pronchery #define SA_BLOCK_MASK (SA_BLOCK_MAX - 1)
47*b077aed3SPierre Pronchery #define SA_BLOCK_MAX_LEVELS (((int)sizeof(ossl_uintmax_t) * 8 \
48*b077aed3SPierre Pronchery + OPENSSL_SA_BLOCK_BITS - 1) \
49*b077aed3SPierre Pronchery / OPENSSL_SA_BLOCK_BITS)
50*b077aed3SPierre Pronchery
51*b077aed3SPierre Pronchery struct sparse_array_st {
52*b077aed3SPierre Pronchery int levels;
53*b077aed3SPierre Pronchery ossl_uintmax_t top;
54*b077aed3SPierre Pronchery size_t nelem;
55*b077aed3SPierre Pronchery void **nodes;
56*b077aed3SPierre Pronchery };
57*b077aed3SPierre Pronchery
ossl_sa_new(void)58*b077aed3SPierre Pronchery OPENSSL_SA *ossl_sa_new(void)
59*b077aed3SPierre Pronchery {
60*b077aed3SPierre Pronchery OPENSSL_SA *res = OPENSSL_zalloc(sizeof(*res));
61*b077aed3SPierre Pronchery
62*b077aed3SPierre Pronchery return res;
63*b077aed3SPierre Pronchery }
64*b077aed3SPierre Pronchery
sa_doall(const OPENSSL_SA * sa,void (* node)(void **),void (* leaf)(ossl_uintmax_t,void *,void *),void * arg)65*b077aed3SPierre Pronchery static void sa_doall(const OPENSSL_SA *sa, void (*node)(void **),
66*b077aed3SPierre Pronchery void (*leaf)(ossl_uintmax_t, void *, void *), void *arg)
67*b077aed3SPierre Pronchery {
68*b077aed3SPierre Pronchery int i[SA_BLOCK_MAX_LEVELS];
69*b077aed3SPierre Pronchery void *nodes[SA_BLOCK_MAX_LEVELS];
70*b077aed3SPierre Pronchery ossl_uintmax_t idx = 0;
71*b077aed3SPierre Pronchery int l = 0;
72*b077aed3SPierre Pronchery
73*b077aed3SPierre Pronchery i[0] = 0;
74*b077aed3SPierre Pronchery nodes[0] = sa->nodes;
75*b077aed3SPierre Pronchery while (l >= 0) {
76*b077aed3SPierre Pronchery const int n = i[l];
77*b077aed3SPierre Pronchery void ** const p = nodes[l];
78*b077aed3SPierre Pronchery
79*b077aed3SPierre Pronchery if (n >= SA_BLOCK_MAX) {
80*b077aed3SPierre Pronchery if (p != NULL && node != NULL)
81*b077aed3SPierre Pronchery (*node)(p);
82*b077aed3SPierre Pronchery l--;
83*b077aed3SPierre Pronchery idx >>= OPENSSL_SA_BLOCK_BITS;
84*b077aed3SPierre Pronchery } else {
85*b077aed3SPierre Pronchery i[l] = n + 1;
86*b077aed3SPierre Pronchery if (p != NULL && p[n] != NULL) {
87*b077aed3SPierre Pronchery idx = (idx & ~SA_BLOCK_MASK) | n;
88*b077aed3SPierre Pronchery if (l < sa->levels - 1) {
89*b077aed3SPierre Pronchery i[++l] = 0;
90*b077aed3SPierre Pronchery nodes[l] = p[n];
91*b077aed3SPierre Pronchery idx <<= OPENSSL_SA_BLOCK_BITS;
92*b077aed3SPierre Pronchery } else if (leaf != NULL) {
93*b077aed3SPierre Pronchery (*leaf)(idx, p[n], arg);
94*b077aed3SPierre Pronchery }
95*b077aed3SPierre Pronchery }
96*b077aed3SPierre Pronchery }
97*b077aed3SPierre Pronchery }
98*b077aed3SPierre Pronchery }
99*b077aed3SPierre Pronchery
sa_free_node(void ** p)100*b077aed3SPierre Pronchery static void sa_free_node(void **p)
101*b077aed3SPierre Pronchery {
102*b077aed3SPierre Pronchery OPENSSL_free(p);
103*b077aed3SPierre Pronchery }
104*b077aed3SPierre Pronchery
sa_free_leaf(ossl_uintmax_t n,void * p,void * arg)105*b077aed3SPierre Pronchery static void sa_free_leaf(ossl_uintmax_t n, void *p, void *arg)
106*b077aed3SPierre Pronchery {
107*b077aed3SPierre Pronchery OPENSSL_free(p);
108*b077aed3SPierre Pronchery }
109*b077aed3SPierre Pronchery
ossl_sa_free(OPENSSL_SA * sa)110*b077aed3SPierre Pronchery void ossl_sa_free(OPENSSL_SA *sa)
111*b077aed3SPierre Pronchery {
112*b077aed3SPierre Pronchery if (sa != NULL) {
113*b077aed3SPierre Pronchery sa_doall(sa, &sa_free_node, NULL, NULL);
114*b077aed3SPierre Pronchery OPENSSL_free(sa);
115*b077aed3SPierre Pronchery }
116*b077aed3SPierre Pronchery }
117*b077aed3SPierre Pronchery
ossl_sa_free_leaves(OPENSSL_SA * sa)118*b077aed3SPierre Pronchery void ossl_sa_free_leaves(OPENSSL_SA *sa)
119*b077aed3SPierre Pronchery {
120*b077aed3SPierre Pronchery sa_doall(sa, &sa_free_node, &sa_free_leaf, NULL);
121*b077aed3SPierre Pronchery OPENSSL_free(sa);
122*b077aed3SPierre Pronchery }
123*b077aed3SPierre Pronchery
124*b077aed3SPierre Pronchery /* Wrap this in a structure to avoid compiler warnings */
125*b077aed3SPierre Pronchery struct trampoline_st {
126*b077aed3SPierre Pronchery void (*func)(ossl_uintmax_t, void *);
127*b077aed3SPierre Pronchery };
128*b077aed3SPierre Pronchery
trampoline(ossl_uintmax_t n,void * l,void * arg)129*b077aed3SPierre Pronchery static void trampoline(ossl_uintmax_t n, void *l, void *arg)
130*b077aed3SPierre Pronchery {
131*b077aed3SPierre Pronchery ((const struct trampoline_st *)arg)->func(n, l);
132*b077aed3SPierre Pronchery }
133*b077aed3SPierre Pronchery
ossl_sa_doall(const OPENSSL_SA * sa,void (* leaf)(ossl_uintmax_t,void *))134*b077aed3SPierre Pronchery void ossl_sa_doall(const OPENSSL_SA *sa, void (*leaf)(ossl_uintmax_t, void *))
135*b077aed3SPierre Pronchery {
136*b077aed3SPierre Pronchery struct trampoline_st tramp;
137*b077aed3SPierre Pronchery
138*b077aed3SPierre Pronchery tramp.func = leaf;
139*b077aed3SPierre Pronchery if (sa != NULL)
140*b077aed3SPierre Pronchery sa_doall(sa, NULL, &trampoline, &tramp);
141*b077aed3SPierre Pronchery }
142*b077aed3SPierre Pronchery
ossl_sa_doall_arg(const OPENSSL_SA * sa,void (* leaf)(ossl_uintmax_t,void *,void *),void * arg)143*b077aed3SPierre Pronchery void ossl_sa_doall_arg(const OPENSSL_SA *sa,
144*b077aed3SPierre Pronchery void (*leaf)(ossl_uintmax_t, void *, void *),
145*b077aed3SPierre Pronchery void *arg)
146*b077aed3SPierre Pronchery {
147*b077aed3SPierre Pronchery if (sa != NULL)
148*b077aed3SPierre Pronchery sa_doall(sa, NULL, leaf, arg);
149*b077aed3SPierre Pronchery }
150*b077aed3SPierre Pronchery
ossl_sa_num(const OPENSSL_SA * sa)151*b077aed3SPierre Pronchery size_t ossl_sa_num(const OPENSSL_SA *sa)
152*b077aed3SPierre Pronchery {
153*b077aed3SPierre Pronchery return sa == NULL ? 0 : sa->nelem;
154*b077aed3SPierre Pronchery }
155*b077aed3SPierre Pronchery
ossl_sa_get(const OPENSSL_SA * sa,ossl_uintmax_t n)156*b077aed3SPierre Pronchery void *ossl_sa_get(const OPENSSL_SA *sa, ossl_uintmax_t n)
157*b077aed3SPierre Pronchery {
158*b077aed3SPierre Pronchery int level;
159*b077aed3SPierre Pronchery void **p, *r = NULL;
160*b077aed3SPierre Pronchery
161*b077aed3SPierre Pronchery if (sa == NULL || sa->nelem == 0)
162*b077aed3SPierre Pronchery return NULL;
163*b077aed3SPierre Pronchery
164*b077aed3SPierre Pronchery if (n <= sa->top) {
165*b077aed3SPierre Pronchery p = sa->nodes;
166*b077aed3SPierre Pronchery for (level = sa->levels - 1; p != NULL && level > 0; level--)
167*b077aed3SPierre Pronchery p = (void **)p[(n >> (OPENSSL_SA_BLOCK_BITS * level))
168*b077aed3SPierre Pronchery & SA_BLOCK_MASK];
169*b077aed3SPierre Pronchery r = p == NULL ? NULL : p[n & SA_BLOCK_MASK];
170*b077aed3SPierre Pronchery }
171*b077aed3SPierre Pronchery return r;
172*b077aed3SPierre Pronchery }
173*b077aed3SPierre Pronchery
alloc_node(void)174*b077aed3SPierre Pronchery static ossl_inline void **alloc_node(void)
175*b077aed3SPierre Pronchery {
176*b077aed3SPierre Pronchery return OPENSSL_zalloc(SA_BLOCK_MAX * sizeof(void *));
177*b077aed3SPierre Pronchery }
178*b077aed3SPierre Pronchery
ossl_sa_set(OPENSSL_SA * sa,ossl_uintmax_t posn,void * val)179*b077aed3SPierre Pronchery int ossl_sa_set(OPENSSL_SA *sa, ossl_uintmax_t posn, void *val)
180*b077aed3SPierre Pronchery {
181*b077aed3SPierre Pronchery int i, level = 1;
182*b077aed3SPierre Pronchery ossl_uintmax_t n = posn;
183*b077aed3SPierre Pronchery void **p;
184*b077aed3SPierre Pronchery
185*b077aed3SPierre Pronchery if (sa == NULL)
186*b077aed3SPierre Pronchery return 0;
187*b077aed3SPierre Pronchery
188*b077aed3SPierre Pronchery for (level = 1; level < SA_BLOCK_MAX_LEVELS; level++)
189*b077aed3SPierre Pronchery if ((n >>= OPENSSL_SA_BLOCK_BITS) == 0)
190*b077aed3SPierre Pronchery break;
191*b077aed3SPierre Pronchery
192*b077aed3SPierre Pronchery for (;sa->levels < level; sa->levels++) {
193*b077aed3SPierre Pronchery p = alloc_node();
194*b077aed3SPierre Pronchery if (p == NULL)
195*b077aed3SPierre Pronchery return 0;
196*b077aed3SPierre Pronchery p[0] = sa->nodes;
197*b077aed3SPierre Pronchery sa->nodes = p;
198*b077aed3SPierre Pronchery }
199*b077aed3SPierre Pronchery if (sa->top < posn)
200*b077aed3SPierre Pronchery sa->top = posn;
201*b077aed3SPierre Pronchery
202*b077aed3SPierre Pronchery p = sa->nodes;
203*b077aed3SPierre Pronchery for (level = sa->levels - 1; level > 0; level--) {
204*b077aed3SPierre Pronchery i = (posn >> (OPENSSL_SA_BLOCK_BITS * level)) & SA_BLOCK_MASK;
205*b077aed3SPierre Pronchery if (p[i] == NULL && (p[i] = alloc_node()) == NULL)
206*b077aed3SPierre Pronchery return 0;
207*b077aed3SPierre Pronchery p = p[i];
208*b077aed3SPierre Pronchery }
209*b077aed3SPierre Pronchery p += posn & SA_BLOCK_MASK;
210*b077aed3SPierre Pronchery if (val == NULL && *p != NULL)
211*b077aed3SPierre Pronchery sa->nelem--;
212*b077aed3SPierre Pronchery else if (val != NULL && *p == NULL)
213*b077aed3SPierre Pronchery sa->nelem++;
214*b077aed3SPierre Pronchery *p = val;
215*b077aed3SPierre Pronchery return 1;
216*b077aed3SPierre Pronchery }
217