xref: /freebsd-src/crypto/openssl/crypto/sparse_array.c (revision b077aed33b7b6aefca7b17ddb250cf521f938613)
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