1b077aed3SPierre Pronchery /*
2*e0c4386eSCy Schubert * Copyright 2016-2024 The OpenSSL Project Authors. All Rights Reserved.
3b077aed3SPierre Pronchery *
4b077aed3SPierre Pronchery * Licensed under the Apache License 2.0 (the "License"). You may not use
5b077aed3SPierre Pronchery * this file except in compliance with the License. You can obtain a copy
6b077aed3SPierre Pronchery * in the file LICENSE in the source distribution or at
7b077aed3SPierre Pronchery * https://www.openssl.org/source/license.html
8b077aed3SPierre Pronchery */
9b077aed3SPierre Pronchery
10b077aed3SPierre Pronchery /*
11b077aed3SPierre Pronchery * Derived from the BLAKE2 reference implementation written by Samuel Neves.
12b077aed3SPierre Pronchery * Copyright 2012, Samuel Neves <sneves@dei.uc.pt>
13b077aed3SPierre Pronchery * More information about the BLAKE2 hash function and its implementations
14b077aed3SPierre Pronchery * can be found at https://blake2.net.
15b077aed3SPierre Pronchery */
16b077aed3SPierre Pronchery
17b077aed3SPierre Pronchery #include <assert.h>
18b077aed3SPierre Pronchery #include <string.h>
19b077aed3SPierre Pronchery #include <openssl/crypto.h>
20b077aed3SPierre Pronchery #include "blake2_impl.h"
21b077aed3SPierre Pronchery #include "prov/blake2.h"
22b077aed3SPierre Pronchery
23b077aed3SPierre Pronchery static const uint32_t blake2s_IV[8] =
24b077aed3SPierre Pronchery {
25b077aed3SPierre Pronchery 0x6A09E667U, 0xBB67AE85U, 0x3C6EF372U, 0xA54FF53AU,
26b077aed3SPierre Pronchery 0x510E527FU, 0x9B05688CU, 0x1F83D9ABU, 0x5BE0CD19U
27b077aed3SPierre Pronchery };
28b077aed3SPierre Pronchery
29b077aed3SPierre Pronchery static const uint8_t blake2s_sigma[10][16] =
30b077aed3SPierre Pronchery {
31b077aed3SPierre Pronchery { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 } ,
32b077aed3SPierre Pronchery { 14, 10, 4, 8, 9, 15, 13, 6, 1, 12, 0, 2, 11, 7, 5, 3 } ,
33b077aed3SPierre Pronchery { 11, 8, 12, 0, 5, 2, 15, 13, 10, 14, 3, 6, 7, 1, 9, 4 } ,
34b077aed3SPierre Pronchery { 7, 9, 3, 1, 13, 12, 11, 14, 2, 6, 5, 10, 4, 0, 15, 8 } ,
35b077aed3SPierre Pronchery { 9, 0, 5, 7, 2, 4, 10, 15, 14, 1, 11, 12, 6, 8, 3, 13 } ,
36b077aed3SPierre Pronchery { 2, 12, 6, 10, 0, 11, 8, 3, 4, 13, 7, 5, 15, 14, 1, 9 } ,
37b077aed3SPierre Pronchery { 12, 5, 1, 15, 14, 13, 4, 10, 0, 7, 6, 3, 9, 2, 8, 11 } ,
38b077aed3SPierre Pronchery { 13, 11, 7, 14, 12, 1, 3, 9, 5, 0, 15, 4, 8, 6, 2, 10 } ,
39b077aed3SPierre Pronchery { 6, 15, 14, 9, 11, 3, 0, 8, 12, 2, 13, 7, 1, 4, 10, 5 } ,
40b077aed3SPierre Pronchery { 10, 2, 8, 4, 7, 6, 1, 5, 15, 11, 9, 14, 3, 12, 13 , 0 } ,
41b077aed3SPierre Pronchery };
42b077aed3SPierre Pronchery
43b077aed3SPierre Pronchery /* Set that it's the last block we'll compress */
blake2s_set_lastblock(BLAKE2S_CTX * S)44b077aed3SPierre Pronchery static ossl_inline void blake2s_set_lastblock(BLAKE2S_CTX *S)
45b077aed3SPierre Pronchery {
46b077aed3SPierre Pronchery S->f[0] = -1;
47b077aed3SPierre Pronchery }
48b077aed3SPierre Pronchery
49b077aed3SPierre Pronchery /* Initialize the hashing state. */
blake2s_init0(BLAKE2S_CTX * S)50b077aed3SPierre Pronchery static ossl_inline void blake2s_init0(BLAKE2S_CTX *S)
51b077aed3SPierre Pronchery {
52b077aed3SPierre Pronchery int i;
53b077aed3SPierre Pronchery
54b077aed3SPierre Pronchery memset(S, 0, sizeof(BLAKE2S_CTX));
55b077aed3SPierre Pronchery for (i = 0; i < 8; ++i) {
56b077aed3SPierre Pronchery S->h[i] = blake2s_IV[i];
57b077aed3SPierre Pronchery }
58b077aed3SPierre Pronchery }
59b077aed3SPierre Pronchery
60b077aed3SPierre Pronchery /* init xors IV with input parameter block and sets the output length */
blake2s_init_param(BLAKE2S_CTX * S,const BLAKE2S_PARAM * P)61b077aed3SPierre Pronchery static void blake2s_init_param(BLAKE2S_CTX *S, const BLAKE2S_PARAM *P)
62b077aed3SPierre Pronchery {
63b077aed3SPierre Pronchery size_t i;
64b077aed3SPierre Pronchery const uint8_t *p = (const uint8_t *)(P);
65b077aed3SPierre Pronchery
66b077aed3SPierre Pronchery blake2s_init0(S);
67b077aed3SPierre Pronchery S->outlen = P->digest_length;
68b077aed3SPierre Pronchery
69b077aed3SPierre Pronchery /* The param struct is carefully hand packed, and should be 32 bytes on
70b077aed3SPierre Pronchery * every platform. */
71b077aed3SPierre Pronchery assert(sizeof(BLAKE2S_PARAM) == 32);
72b077aed3SPierre Pronchery /* IV XOR ParamBlock */
73b077aed3SPierre Pronchery for (i = 0; i < 8; ++i) {
74b077aed3SPierre Pronchery S->h[i] ^= load32(&p[i*4]);
75b077aed3SPierre Pronchery }
76b077aed3SPierre Pronchery }
77b077aed3SPierre Pronchery
ossl_blake2s_param_init(BLAKE2S_PARAM * P)78b077aed3SPierre Pronchery void ossl_blake2s_param_init(BLAKE2S_PARAM *P)
79b077aed3SPierre Pronchery {
80b077aed3SPierre Pronchery P->digest_length = BLAKE2S_DIGEST_LENGTH;
81b077aed3SPierre Pronchery P->key_length = 0;
82b077aed3SPierre Pronchery P->fanout = 1;
83b077aed3SPierre Pronchery P->depth = 1;
84b077aed3SPierre Pronchery store32(P->leaf_length, 0);
85b077aed3SPierre Pronchery store48(P->node_offset, 0);
86b077aed3SPierre Pronchery P->node_depth = 0;
87b077aed3SPierre Pronchery P->inner_length = 0;
88b077aed3SPierre Pronchery memset(P->salt, 0, sizeof(P->salt));
89b077aed3SPierre Pronchery memset(P->personal, 0, sizeof(P->personal));
90b077aed3SPierre Pronchery }
91b077aed3SPierre Pronchery
ossl_blake2s_param_set_digest_length(BLAKE2S_PARAM * P,uint8_t outlen)92b077aed3SPierre Pronchery void ossl_blake2s_param_set_digest_length(BLAKE2S_PARAM *P, uint8_t outlen)
93b077aed3SPierre Pronchery {
94b077aed3SPierre Pronchery P->digest_length = outlen;
95b077aed3SPierre Pronchery }
96b077aed3SPierre Pronchery
ossl_blake2s_param_set_key_length(BLAKE2S_PARAM * P,uint8_t keylen)97b077aed3SPierre Pronchery void ossl_blake2s_param_set_key_length(BLAKE2S_PARAM *P, uint8_t keylen)
98b077aed3SPierre Pronchery {
99b077aed3SPierre Pronchery P->key_length = keylen;
100b077aed3SPierre Pronchery }
101b077aed3SPierre Pronchery
ossl_blake2s_param_set_personal(BLAKE2S_PARAM * P,const uint8_t * personal,size_t len)102b077aed3SPierre Pronchery void ossl_blake2s_param_set_personal(BLAKE2S_PARAM *P, const uint8_t *personal,
103b077aed3SPierre Pronchery size_t len)
104b077aed3SPierre Pronchery {
105b077aed3SPierre Pronchery memcpy(P->personal, personal, len);
106b077aed3SPierre Pronchery memset(P->personal + len, 0, BLAKE2S_PERSONALBYTES - len);
107b077aed3SPierre Pronchery }
108b077aed3SPierre Pronchery
ossl_blake2s_param_set_salt(BLAKE2S_PARAM * P,const uint8_t * salt,size_t len)109b077aed3SPierre Pronchery void ossl_blake2s_param_set_salt(BLAKE2S_PARAM *P, const uint8_t *salt,
110b077aed3SPierre Pronchery size_t len)
111b077aed3SPierre Pronchery {
112b077aed3SPierre Pronchery memcpy(P->salt, salt, len);
113b077aed3SPierre Pronchery memset(P->salt + len, 0, BLAKE2S_SALTBYTES - len);}
114b077aed3SPierre Pronchery
115b077aed3SPierre Pronchery /*
116b077aed3SPierre Pronchery * Initialize the hashing context with the given parameter block.
117b077aed3SPierre Pronchery * Always returns 1.
118b077aed3SPierre Pronchery */
ossl_blake2s_init(BLAKE2S_CTX * c,const BLAKE2S_PARAM * P)119b077aed3SPierre Pronchery int ossl_blake2s_init(BLAKE2S_CTX *c, const BLAKE2S_PARAM *P)
120b077aed3SPierre Pronchery {
121b077aed3SPierre Pronchery blake2s_init_param(c, P);
122b077aed3SPierre Pronchery return 1;
123b077aed3SPierre Pronchery }
124b077aed3SPierre Pronchery
125b077aed3SPierre Pronchery /*
126b077aed3SPierre Pronchery * Initialize the hashing context with the given parameter block and key.
127b077aed3SPierre Pronchery * Always returns 1.
128b077aed3SPierre Pronchery */
ossl_blake2s_init_key(BLAKE2S_CTX * c,const BLAKE2S_PARAM * P,const void * key)129b077aed3SPierre Pronchery int ossl_blake2s_init_key(BLAKE2S_CTX *c, const BLAKE2S_PARAM *P,
130b077aed3SPierre Pronchery const void *key)
131b077aed3SPierre Pronchery {
132b077aed3SPierre Pronchery blake2s_init_param(c, P);
133b077aed3SPierre Pronchery
134b077aed3SPierre Pronchery /* Pad the key to form first data block */
135b077aed3SPierre Pronchery {
136b077aed3SPierre Pronchery uint8_t block[BLAKE2S_BLOCKBYTES] = {0};
137b077aed3SPierre Pronchery
138b077aed3SPierre Pronchery memcpy(block, key, P->key_length);
139b077aed3SPierre Pronchery ossl_blake2s_update(c, block, BLAKE2S_BLOCKBYTES);
140b077aed3SPierre Pronchery OPENSSL_cleanse(block, BLAKE2S_BLOCKBYTES);
141b077aed3SPierre Pronchery }
142b077aed3SPierre Pronchery
143b077aed3SPierre Pronchery return 1;
144b077aed3SPierre Pronchery }
145b077aed3SPierre Pronchery
146b077aed3SPierre Pronchery /* Permute the state while xoring in the block of data. */
blake2s_compress(BLAKE2S_CTX * S,const uint8_t * blocks,size_t len)147b077aed3SPierre Pronchery static void blake2s_compress(BLAKE2S_CTX *S,
148b077aed3SPierre Pronchery const uint8_t *blocks,
149b077aed3SPierre Pronchery size_t len)
150b077aed3SPierre Pronchery {
151b077aed3SPierre Pronchery uint32_t m[16];
152b077aed3SPierre Pronchery uint32_t v[16];
153b077aed3SPierre Pronchery size_t i;
154b077aed3SPierre Pronchery size_t increment;
155b077aed3SPierre Pronchery
156b077aed3SPierre Pronchery /*
157b077aed3SPierre Pronchery * There are two distinct usage vectors for this function:
158b077aed3SPierre Pronchery *
159b077aed3SPierre Pronchery * a) BLAKE2s_Update uses it to process complete blocks,
160b077aed3SPierre Pronchery * possibly more than one at a time;
161b077aed3SPierre Pronchery *
162b077aed3SPierre Pronchery * b) BLAK2s_Final uses it to process last block, always
163b077aed3SPierre Pronchery * single but possibly incomplete, in which case caller
164b077aed3SPierre Pronchery * pads input with zeros.
165b077aed3SPierre Pronchery */
166b077aed3SPierre Pronchery assert(len < BLAKE2S_BLOCKBYTES || len % BLAKE2S_BLOCKBYTES == 0);
167b077aed3SPierre Pronchery
168b077aed3SPierre Pronchery /*
169b077aed3SPierre Pronchery * Since last block is always processed with separate call,
170b077aed3SPierre Pronchery * |len| not being multiple of complete blocks can be observed
171b077aed3SPierre Pronchery * only with |len| being less than BLAKE2S_BLOCKBYTES ("less"
172b077aed3SPierre Pronchery * including even zero), which is why following assignment doesn't
173b077aed3SPierre Pronchery * have to reside inside the main loop below.
174b077aed3SPierre Pronchery */
175b077aed3SPierre Pronchery increment = len < BLAKE2S_BLOCKBYTES ? len : BLAKE2S_BLOCKBYTES;
176b077aed3SPierre Pronchery
177b077aed3SPierre Pronchery for (i = 0; i < 8; ++i) {
178b077aed3SPierre Pronchery v[i] = S->h[i];
179b077aed3SPierre Pronchery }
180b077aed3SPierre Pronchery
181b077aed3SPierre Pronchery do {
182b077aed3SPierre Pronchery for (i = 0; i < 16; ++i) {
183b077aed3SPierre Pronchery m[i] = load32(blocks + i * sizeof(m[i]));
184b077aed3SPierre Pronchery }
185b077aed3SPierre Pronchery
186b077aed3SPierre Pronchery /* blake2s_increment_counter */
187b077aed3SPierre Pronchery S->t[0] += increment;
188b077aed3SPierre Pronchery S->t[1] += (S->t[0] < increment);
189b077aed3SPierre Pronchery
190b077aed3SPierre Pronchery v[ 8] = blake2s_IV[0];
191b077aed3SPierre Pronchery v[ 9] = blake2s_IV[1];
192b077aed3SPierre Pronchery v[10] = blake2s_IV[2];
193b077aed3SPierre Pronchery v[11] = blake2s_IV[3];
194b077aed3SPierre Pronchery v[12] = S->t[0] ^ blake2s_IV[4];
195b077aed3SPierre Pronchery v[13] = S->t[1] ^ blake2s_IV[5];
196b077aed3SPierre Pronchery v[14] = S->f[0] ^ blake2s_IV[6];
197b077aed3SPierre Pronchery v[15] = S->f[1] ^ blake2s_IV[7];
198b077aed3SPierre Pronchery #define G(r,i,a,b,c,d) \
199b077aed3SPierre Pronchery do { \
200b077aed3SPierre Pronchery a = a + b + m[blake2s_sigma[r][2*i+0]]; \
201b077aed3SPierre Pronchery d = rotr32(d ^ a, 16); \
202b077aed3SPierre Pronchery c = c + d; \
203b077aed3SPierre Pronchery b = rotr32(b ^ c, 12); \
204b077aed3SPierre Pronchery a = a + b + m[blake2s_sigma[r][2*i+1]]; \
205b077aed3SPierre Pronchery d = rotr32(d ^ a, 8); \
206b077aed3SPierre Pronchery c = c + d; \
207b077aed3SPierre Pronchery b = rotr32(b ^ c, 7); \
208b077aed3SPierre Pronchery } while (0)
209b077aed3SPierre Pronchery #define ROUND(r) \
210b077aed3SPierre Pronchery do { \
211b077aed3SPierre Pronchery G(r,0,v[ 0],v[ 4],v[ 8],v[12]); \
212b077aed3SPierre Pronchery G(r,1,v[ 1],v[ 5],v[ 9],v[13]); \
213b077aed3SPierre Pronchery G(r,2,v[ 2],v[ 6],v[10],v[14]); \
214b077aed3SPierre Pronchery G(r,3,v[ 3],v[ 7],v[11],v[15]); \
215b077aed3SPierre Pronchery G(r,4,v[ 0],v[ 5],v[10],v[15]); \
216b077aed3SPierre Pronchery G(r,5,v[ 1],v[ 6],v[11],v[12]); \
217b077aed3SPierre Pronchery G(r,6,v[ 2],v[ 7],v[ 8],v[13]); \
218b077aed3SPierre Pronchery G(r,7,v[ 3],v[ 4],v[ 9],v[14]); \
219b077aed3SPierre Pronchery } while (0)
220b077aed3SPierre Pronchery #if defined(OPENSSL_SMALL_FOOTPRINT)
221b077aed3SPierre Pronchery /* almost 3x reduction on x86_64, 4.5x on ARMv8, 4x on ARMv4 */
222b077aed3SPierre Pronchery for (i = 0; i < 10; i++) {
223b077aed3SPierre Pronchery ROUND(i);
224b077aed3SPierre Pronchery }
225b077aed3SPierre Pronchery #else
226b077aed3SPierre Pronchery ROUND(0);
227b077aed3SPierre Pronchery ROUND(1);
228b077aed3SPierre Pronchery ROUND(2);
229b077aed3SPierre Pronchery ROUND(3);
230b077aed3SPierre Pronchery ROUND(4);
231b077aed3SPierre Pronchery ROUND(5);
232b077aed3SPierre Pronchery ROUND(6);
233b077aed3SPierre Pronchery ROUND(7);
234b077aed3SPierre Pronchery ROUND(8);
235b077aed3SPierre Pronchery ROUND(9);
236b077aed3SPierre Pronchery #endif
237b077aed3SPierre Pronchery
238b077aed3SPierre Pronchery for (i = 0; i < 8; ++i) {
239b077aed3SPierre Pronchery S->h[i] = v[i] ^= v[i + 8] ^ S->h[i];
240b077aed3SPierre Pronchery }
241b077aed3SPierre Pronchery #undef G
242b077aed3SPierre Pronchery #undef ROUND
243b077aed3SPierre Pronchery blocks += increment;
244b077aed3SPierre Pronchery len -= increment;
245b077aed3SPierre Pronchery } while (len);
246b077aed3SPierre Pronchery }
247b077aed3SPierre Pronchery
248b077aed3SPierre Pronchery /* Absorb the input data into the hash state. Always returns 1. */
ossl_blake2s_update(BLAKE2S_CTX * c,const void * data,size_t datalen)249b077aed3SPierre Pronchery int ossl_blake2s_update(BLAKE2S_CTX *c, const void *data, size_t datalen)
250b077aed3SPierre Pronchery {
251b077aed3SPierre Pronchery const uint8_t *in = data;
252b077aed3SPierre Pronchery size_t fill;
253b077aed3SPierre Pronchery
254b077aed3SPierre Pronchery /*
255b077aed3SPierre Pronchery * Intuitively one would expect intermediate buffer, c->buf, to
256b077aed3SPierre Pronchery * store incomplete blocks. But in this case we are interested to
257b077aed3SPierre Pronchery * temporarily stash even complete blocks, because last one in the
258b077aed3SPierre Pronchery * stream has to be treated in special way, and at this point we
259b077aed3SPierre Pronchery * don't know if last block in *this* call is last one "ever". This
260b077aed3SPierre Pronchery * is the reason for why |datalen| is compared as >, and not >=.
261b077aed3SPierre Pronchery */
262b077aed3SPierre Pronchery fill = sizeof(c->buf) - c->buflen;
263b077aed3SPierre Pronchery if (datalen > fill) {
264b077aed3SPierre Pronchery if (c->buflen) {
265b077aed3SPierre Pronchery memcpy(c->buf + c->buflen, in, fill); /* Fill buffer */
266b077aed3SPierre Pronchery blake2s_compress(c, c->buf, BLAKE2S_BLOCKBYTES);
267b077aed3SPierre Pronchery c->buflen = 0;
268b077aed3SPierre Pronchery in += fill;
269b077aed3SPierre Pronchery datalen -= fill;
270b077aed3SPierre Pronchery }
271b077aed3SPierre Pronchery if (datalen > BLAKE2S_BLOCKBYTES) {
272b077aed3SPierre Pronchery size_t stashlen = datalen % BLAKE2S_BLOCKBYTES;
273b077aed3SPierre Pronchery /*
274b077aed3SPierre Pronchery * If |datalen| is a multiple of the blocksize, stash
275b077aed3SPierre Pronchery * last complete block, it can be final one...
276b077aed3SPierre Pronchery */
277b077aed3SPierre Pronchery stashlen = stashlen ? stashlen : BLAKE2S_BLOCKBYTES;
278b077aed3SPierre Pronchery datalen -= stashlen;
279b077aed3SPierre Pronchery blake2s_compress(c, in, datalen);
280b077aed3SPierre Pronchery in += datalen;
281b077aed3SPierre Pronchery datalen = stashlen;
282b077aed3SPierre Pronchery }
283b077aed3SPierre Pronchery }
284b077aed3SPierre Pronchery
285b077aed3SPierre Pronchery assert(datalen <= BLAKE2S_BLOCKBYTES);
286b077aed3SPierre Pronchery
287b077aed3SPierre Pronchery memcpy(c->buf + c->buflen, in, datalen);
288b077aed3SPierre Pronchery c->buflen += datalen; /* Be lazy, do not compress */
289b077aed3SPierre Pronchery
290b077aed3SPierre Pronchery return 1;
291b077aed3SPierre Pronchery }
292b077aed3SPierre Pronchery
293b077aed3SPierre Pronchery /*
294b077aed3SPierre Pronchery * Calculate the final hash and save it in md.
295b077aed3SPierre Pronchery * Always returns 1.
296b077aed3SPierre Pronchery */
ossl_blake2s_final(unsigned char * md,BLAKE2S_CTX * c)297b077aed3SPierre Pronchery int ossl_blake2s_final(unsigned char *md, BLAKE2S_CTX *c)
298b077aed3SPierre Pronchery {
299b077aed3SPierre Pronchery uint8_t outbuffer[BLAKE2S_OUTBYTES] = {0};
300b077aed3SPierre Pronchery uint8_t *target = outbuffer;
301b077aed3SPierre Pronchery int iter = (c->outlen + 3) / 4;
302b077aed3SPierre Pronchery int i;
303b077aed3SPierre Pronchery
304b077aed3SPierre Pronchery /* Avoid writing to the temporary buffer if possible */
305b077aed3SPierre Pronchery if ((c->outlen % sizeof(c->h[0])) == 0)
306b077aed3SPierre Pronchery target = md;
307b077aed3SPierre Pronchery
308b077aed3SPierre Pronchery blake2s_set_lastblock(c);
309b077aed3SPierre Pronchery /* Padding */
310b077aed3SPierre Pronchery memset(c->buf + c->buflen, 0, sizeof(c->buf) - c->buflen);
311b077aed3SPierre Pronchery blake2s_compress(c, c->buf, c->buflen);
312b077aed3SPierre Pronchery
313b077aed3SPierre Pronchery /* Output full hash to buffer */
314b077aed3SPierre Pronchery for (i = 0; i < iter; ++i)
315b077aed3SPierre Pronchery store32(target + sizeof(c->h[i]) * i, c->h[i]);
316b077aed3SPierre Pronchery
317*e0c4386eSCy Schubert if (target != md) {
318b077aed3SPierre Pronchery memcpy(md, target, c->outlen);
319*e0c4386eSCy Schubert OPENSSL_cleanse(target, sizeof(outbuffer));
320*e0c4386eSCy Schubert }
321b077aed3SPierre Pronchery
322b077aed3SPierre Pronchery OPENSSL_cleanse(c, sizeof(BLAKE2S_CTX));
323b077aed3SPierre Pronchery return 1;
324b077aed3SPierre Pronchery }
325