1 /*
2    BLAKE2 reference source code package - optimized C implementations
3 
4    Copyright 2012, Samuel Neves <sneves@dei.uc.pt>.  You may use this under the
5    terms of the CC0, the OpenSSL Licence, or the Apache Public License 2.0, at
6    your option.  The terms of these licenses can be found at:
7 
8    - CC0 1.0 Universal : http://creativecommons.org/publicdomain/zero/1.0
9    - OpenSSL license   : https://www.openssl.org/source/license.html
10    - Apache 2.0        : http://www.apache.org/licenses/LICENSE-2.0
11 
12    More information about the BLAKE2 hash function can be found at
13    https://blake2.net.
14 */
15 
16 #include <stdint.h>
17 #include <string.h>
18 #include <stdio.h>
19 
20 #include "blake2.h"
21 #include "blake2-impl.h"
22 
23 #include "blake2-config.h"
24 
25 #ifdef _MSC_VER
26 #include <intrin.h> /* for _mm_set_epi64x */
27 #endif
28 #include <emmintrin.h>
29 #if defined(HAVE_SSSE3)
30 #include <tmmintrin.h>
31 #endif
32 #if defined(HAVE_SSE41)
33 #include <smmintrin.h>
34 #endif
35 #if defined(HAVE_AVX)
36 #include <immintrin.h>
37 #endif
38 #if defined(HAVE_XOP)
39 #include <x86intrin.h>
40 #endif
41 
42 #include "blake2b-round.h"
43 
44 static const uint64_t blake2b_IV[8] =
45 {
46   0x6a09e667f3bcc908ULL, 0xbb67ae8584caa73bULL,
47   0x3c6ef372fe94f82bULL, 0xa54ff53a5f1d36f1ULL,
48   0x510e527fade682d1ULL, 0x9b05688c2b3e6c1fULL,
49   0x1f83d9abfb41bd6bULL, 0x5be0cd19137e2179ULL
50 };
51 
52 static const uint8_t blake2b_sigma[12][16] =
53 {
54   {  0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15 } ,
55   { 14, 10,  4,  8,  9, 15, 13,  6,  1, 12,  0,  2, 11,  7,  5,  3 } ,
56   { 11,  8, 12,  0,  5,  2, 15, 13, 10, 14,  3,  6,  7,  1,  9,  4 } ,
57   {  7,  9,  3,  1, 13, 12, 11, 14,  2,  6,  5, 10,  4,  0, 15,  8 } ,
58   {  9,  0,  5,  7,  2,  4, 10, 15, 14,  1, 11, 12,  6,  8,  3, 13 } ,
59   {  2, 12,  6, 10,  0, 11,  8,  3,  4, 13,  7,  5, 15, 14,  1,  9 } ,
60   { 12,  5,  1, 15, 14, 13,  4, 10,  0,  7,  6,  3,  9,  2,  8, 11 } ,
61   { 13, 11,  7, 14, 12,  1,  3,  9,  5,  0, 15,  4,  8,  6,  2, 10 } ,
62   {  6, 15, 14,  9, 11,  3,  0,  8, 12,  2, 13,  7,  1,  4, 10,  5 } ,
63   { 10,  2,  8,  4,  7,  6,  1,  5, 15, 11,  9, 14,  3, 12, 13 , 0 } ,
64   {  0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15 } ,
65   { 14, 10,  4,  8,  9, 15, 13,  6,  1, 12,  0,  2, 11,  7,  5,  3 }
66 };
67 
68 
69 /* Some helper functions, not necessarily useful */
blake2b_set_lastnode(blake2b_state * S)70 BLAKE2_LOCAL_INLINE(int) blake2b_set_lastnode( blake2b_state *S )
71 {
72   S->f[1] = -1;
73   return 0;
74 }
75 
blake2b_clear_lastnode(blake2b_state * S)76 BLAKE2_LOCAL_INLINE(int) blake2b_clear_lastnode( blake2b_state *S )
77 {
78   S->f[1] = 0;
79   return 0;
80 }
81 
blake2b_is_lastblock(const blake2b_state * S)82 BLAKE2_LOCAL_INLINE(int) blake2b_is_lastblock( const blake2b_state *S )
83 {
84   return S->f[0] != 0;
85 }
86 
blake2b_set_lastblock(blake2b_state * S)87 BLAKE2_LOCAL_INLINE(int) blake2b_set_lastblock( blake2b_state *S )
88 {
89   if( S->last_node ) blake2b_set_lastnode( S );
90 
91   S->f[0] = -1;
92   return 0;
93 }
94 
blake2b_clear_lastblock(blake2b_state * S)95 BLAKE2_LOCAL_INLINE(int) blake2b_clear_lastblock( blake2b_state *S )
96 {
97   if( S->last_node ) blake2b_clear_lastnode( S );
98 
99   S->f[0] = 0;
100   return 0;
101 }
102 
103 
blake2b_increment_counter(blake2b_state * S,const uint64_t inc)104 BLAKE2_LOCAL_INLINE(int) blake2b_increment_counter( blake2b_state *S, const uint64_t inc )
105 {
106 #if __x86_64__
107   /* ADD/ADC chain */
108   __uint128_t t = ( ( __uint128_t )S->t[1] << 64 ) | S->t[0];
109   t += inc;
110   S->t[0] = ( uint64_t )( t >>  0 );
111   S->t[1] = ( uint64_t )( t >> 64 );
112 #else
113   S->t[0] += inc;
114   S->t[1] += ( S->t[0] < inc );
115 #endif
116   return 0;
117 }
118 
119 
120 /* Parameter-related functions */
blake2b_param_set_digest_length(blake2b_param * P,const uint8_t digest_length)121 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_digest_length( blake2b_param *P, const uint8_t digest_length )
122 {
123   P->digest_length = digest_length;
124   return 0;
125 }
126 
blake2b_param_set_fanout(blake2b_param * P,const uint8_t fanout)127 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_fanout( blake2b_param *P, const uint8_t fanout )
128 {
129   P->fanout = fanout;
130   return 0;
131 }
132 
blake2b_param_set_max_depth(blake2b_param * P,const uint8_t depth)133 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_max_depth( blake2b_param *P, const uint8_t depth )
134 {
135   P->depth = depth;
136   return 0;
137 }
138 
blake2b_param_set_leaf_length(blake2b_param * P,const uint32_t leaf_length)139 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_leaf_length( blake2b_param *P, const uint32_t leaf_length )
140 {
141   P->leaf_length = leaf_length;
142   return 0;
143 }
144 
blake2b_param_set_node_offset(blake2b_param * P,const uint64_t node_offset)145 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_node_offset( blake2b_param *P, const uint64_t node_offset )
146 {
147   P->node_offset = node_offset;
148   return 0;
149 }
150 
blake2b_param_set_node_depth(blake2b_param * P,const uint8_t node_depth)151 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_node_depth( blake2b_param *P, const uint8_t node_depth )
152 {
153   P->node_depth = node_depth;
154   return 0;
155 }
156 
blake2b_param_set_inner_length(blake2b_param * P,const uint8_t inner_length)157 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_inner_length( blake2b_param *P, const uint8_t inner_length )
158 {
159   P->inner_length = inner_length;
160   return 0;
161 }
162 
blake2b_param_set_salt(blake2b_param * P,const uint8_t salt[BLAKE2B_SALTBYTES])163 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_salt( blake2b_param *P, const uint8_t salt[BLAKE2B_SALTBYTES] )
164 {
165   memcpy( P->salt, salt, BLAKE2B_SALTBYTES );
166   return 0;
167 }
168 
blake2b_param_set_personal(blake2b_param * P,const uint8_t personal[BLAKE2B_PERSONALBYTES])169 BLAKE2_LOCAL_INLINE(int) blake2b_param_set_personal( blake2b_param *P, const uint8_t personal[BLAKE2B_PERSONALBYTES] )
170 {
171   memcpy( P->personal, personal, BLAKE2B_PERSONALBYTES );
172   return 0;
173 }
174 
blake2b_init0(blake2b_state * S)175 BLAKE2_LOCAL_INLINE(int) blake2b_init0( blake2b_state *S )
176 {
177   int i;
178   memset( S, 0, sizeof( blake2b_state ) );
179 
180   for( i = 0; i < 8; ++i ) S->h[i] = blake2b_IV[i];
181 
182   return 0;
183 }
184 
185 /* init xors IV with input parameter block */
blake2b_init_param(blake2b_state * S,const blake2b_param * P)186 int blake2b_init_param( blake2b_state *S, const blake2b_param *P )
187 {
188   /*blake2b_init0( S ); */
189   const uint8_t * v = ( const uint8_t * )( blake2b_IV );
190   const uint8_t * p = ( const uint8_t * )( P );
191   uint8_t * h = ( uint8_t * )( S->h );
192   int i;
193   /* IV XOR ParamBlock */
194   memset( S, 0, sizeof( blake2b_state ) );
195 
196   for( i = 0; i < BLAKE2B_OUTBYTES; ++i ) h[i] = v[i] ^ p[i];
197 
198   return 0;
199 }
200 
201 
202 /* Some sort of default parameter block initialization, for sequential blake2b */
blake2b_init(blake2b_state * S,const uint8_t outlen)203 int blake2b_init( blake2b_state *S, const uint8_t outlen )
204 {
205   const blake2b_param P =
206   {
207     outlen,
208     0,
209     1,
210     1,
211     0,
212     0,
213     0,
214     0,
215     {0},
216     {0},
217     {0}
218   };
219 
220   if ( ( !outlen ) || ( outlen > BLAKE2B_OUTBYTES ) ) return -1;
221 
222   return blake2b_init_param( S, &P );
223 }
224 
blake2b_init_key(blake2b_state * S,const uint8_t outlen,const void * key,const uint8_t keylen)225 int blake2b_init_key( blake2b_state *S, const uint8_t outlen, const void *key, const uint8_t keylen )
226 {
227   const blake2b_param P =
228   {
229     outlen,
230     keylen,
231     1,
232     1,
233     0,
234     0,
235     0,
236     0,
237     {0},
238     {0},
239     {0}
240   };
241 
242   if ( ( !outlen ) || ( outlen > BLAKE2B_OUTBYTES ) ) return -1;
243 
244   if ( ( !keylen ) || keylen > BLAKE2B_KEYBYTES ) return -1;
245 
246   if( blake2b_init_param( S, &P ) < 0 )
247     return 0;
248 
249   {
250     uint8_t block[BLAKE2B_BLOCKBYTES];
251     memset( block, 0, BLAKE2B_BLOCKBYTES );
252     memcpy( block, key, keylen );
253     blake2b_update( S, block, BLAKE2B_BLOCKBYTES );
254     secure_zero_memory( block, BLAKE2B_BLOCKBYTES ); /* Burn the key from stack */
255   }
256   return 0;
257 }
258 
blake2b_compress(blake2b_state * S,const uint8_t block[BLAKE2B_BLOCKBYTES])259 BLAKE2_LOCAL_INLINE(int) blake2b_compress( blake2b_state *S, const uint8_t block[BLAKE2B_BLOCKBYTES] )
260 {
261   __m128i row1l, row1h;
262   __m128i row2l, row2h;
263   __m128i row3l, row3h;
264   __m128i row4l, row4h;
265   __m128i b0, b1;
266   __m128i t0, t1;
267 #if defined(HAVE_SSSE3) && !defined(HAVE_XOP)
268   const __m128i r16 = _mm_setr_epi8( 2, 3, 4, 5, 6, 7, 0, 1, 10, 11, 12, 13, 14, 15, 8, 9 );
269   const __m128i r24 = _mm_setr_epi8( 3, 4, 5, 6, 7, 0, 1, 2, 11, 12, 13, 14, 15, 8, 9, 10 );
270 #endif
271 #if defined(HAVE_SSE41)
272   const __m128i m0 = LOADU( block + 00 );
273   const __m128i m1 = LOADU( block + 16 );
274   const __m128i m2 = LOADU( block + 32 );
275   const __m128i m3 = LOADU( block + 48 );
276   const __m128i m4 = LOADU( block + 64 );
277   const __m128i m5 = LOADU( block + 80 );
278   const __m128i m6 = LOADU( block + 96 );
279   const __m128i m7 = LOADU( block + 112 );
280 #else
281   const uint64_t  m0 = ( ( uint64_t * )block )[ 0];
282   const uint64_t  m1 = ( ( uint64_t * )block )[ 1];
283   const uint64_t  m2 = ( ( uint64_t * )block )[ 2];
284   const uint64_t  m3 = ( ( uint64_t * )block )[ 3];
285   const uint64_t  m4 = ( ( uint64_t * )block )[ 4];
286   const uint64_t  m5 = ( ( uint64_t * )block )[ 5];
287   const uint64_t  m6 = ( ( uint64_t * )block )[ 6];
288   const uint64_t  m7 = ( ( uint64_t * )block )[ 7];
289   const uint64_t  m8 = ( ( uint64_t * )block )[ 8];
290   const uint64_t  m9 = ( ( uint64_t * )block )[ 9];
291   const uint64_t m10 = ( ( uint64_t * )block )[10];
292   const uint64_t m11 = ( ( uint64_t * )block )[11];
293   const uint64_t m12 = ( ( uint64_t * )block )[12];
294   const uint64_t m13 = ( ( uint64_t * )block )[13];
295   const uint64_t m14 = ( ( uint64_t * )block )[14];
296   const uint64_t m15 = ( ( uint64_t * )block )[15];
297 #endif
298   row1l = LOADU( &S->h[0] );
299   row1h = LOADU( &S->h[2] );
300   row2l = LOADU( &S->h[4] );
301   row2h = LOADU( &S->h[6] );
302   row3l = LOADU( &blake2b_IV[0] );
303   row3h = LOADU( &blake2b_IV[2] );
304   row4l = _mm_xor_si128( LOADU( &blake2b_IV[4] ), LOADU( &S->t[0] ) );
305   row4h = _mm_xor_si128( LOADU( &blake2b_IV[6] ), LOADU( &S->f[0] ) );
306   ROUND( 0 );
307   ROUND( 1 );
308   ROUND( 2 );
309   ROUND( 3 );
310   ROUND( 4 );
311   ROUND( 5 );
312   ROUND( 6 );
313   ROUND( 7 );
314   ROUND( 8 );
315   ROUND( 9 );
316   ROUND( 10 );
317   ROUND( 11 );
318   row1l = _mm_xor_si128( row3l, row1l );
319   row1h = _mm_xor_si128( row3h, row1h );
320   STOREU( &S->h[0], _mm_xor_si128( LOADU( &S->h[0] ), row1l ) );
321   STOREU( &S->h[2], _mm_xor_si128( LOADU( &S->h[2] ), row1h ) );
322   row2l = _mm_xor_si128( row4l, row2l );
323   row2h = _mm_xor_si128( row4h, row2h );
324   STOREU( &S->h[4], _mm_xor_si128( LOADU( &S->h[4] ), row2l ) );
325   STOREU( &S->h[6], _mm_xor_si128( LOADU( &S->h[6] ), row2h ) );
326   return 0;
327 }
328 
329 
blake2b_update(blake2b_state * S,const uint8_t * in,uint64_t inlen)330 int blake2b_update( blake2b_state *S, const uint8_t *in, uint64_t inlen )
331 {
332   while( inlen > 0 )
333   {
334     size_t left = S->buflen;
335     size_t fill = 2 * BLAKE2B_BLOCKBYTES - left;
336 
337     if( inlen > fill )
338     {
339       memcpy( S->buf + left, in, fill ); /* Fill buffer */
340       S->buflen += fill;
341       blake2b_increment_counter( S, BLAKE2B_BLOCKBYTES );
342       blake2b_compress( S, S->buf ); /* Compress */
343       memcpy( S->buf, S->buf + BLAKE2B_BLOCKBYTES, BLAKE2B_BLOCKBYTES ); /* Shift buffer left */
344       S->buflen -= BLAKE2B_BLOCKBYTES;
345       in += fill;
346       inlen -= fill;
347     }
348     else /* inlen <= fill */
349     {
350       memcpy( S->buf + left, in, inlen );
351       S->buflen += inlen; /* Be lazy, do not compress */
352       in += inlen;
353       inlen -= inlen;
354     }
355   }
356 
357   return 0;
358 }
359 
360 
blake2b_final(blake2b_state * S,uint8_t * out,uint8_t outlen)361 int blake2b_final( blake2b_state *S, uint8_t *out, uint8_t outlen )
362 {
363   if( outlen > BLAKE2B_OUTBYTES )
364     return -1;
365 
366   if( blake2b_is_lastblock( S ) )
367     return -1;
368 
369   if( S->buflen > BLAKE2B_BLOCKBYTES )
370   {
371     blake2b_increment_counter( S, BLAKE2B_BLOCKBYTES );
372     blake2b_compress( S, S->buf );
373     S->buflen -= BLAKE2B_BLOCKBYTES;
374     memmove( S->buf, S->buf + BLAKE2B_BLOCKBYTES, S->buflen );
375   }
376 
377   blake2b_increment_counter( S, S->buflen );
378   blake2b_set_lastblock( S );
379   memset( S->buf + S->buflen, 0, 2 * BLAKE2B_BLOCKBYTES - S->buflen ); /* Padding */
380   blake2b_compress( S, S->buf );
381   memcpy( out, &S->h[0], outlen );
382   return 0;
383 }
384 
385 
blake2b(uint8_t * out,const void * in,const void * key,const uint8_t outlen,const uint64_t inlen,uint8_t keylen)386 int blake2b( uint8_t *out, const void *in, const void *key, const uint8_t outlen, const uint64_t inlen, uint8_t keylen )
387 {
388   blake2b_state S[1];
389 
390   /* Verify parameters */
391   if ( NULL == in && inlen > 0 ) return -1;
392 
393   if ( NULL == out ) return -1;
394 
395   if( NULL == key && keylen > 0 ) return -1;
396 
397   if( !outlen || outlen > BLAKE2B_OUTBYTES ) return -1;
398 
399   if( keylen > BLAKE2B_KEYBYTES ) return -1;
400 
401   if( keylen )
402   {
403     if( blake2b_init_key( S, outlen, key, keylen ) < 0 ) return -1;
404   }
405   else
406   {
407     if( blake2b_init( S, outlen ) < 0 ) return -1;
408   }
409 
410   blake2b_update( S, ( const uint8_t * )in, inlen );
411   blake2b_final( S, out, outlen );
412   return 0;
413 }
414 
415 #if defined(SUPERCOP)
crypto_hash(unsigned char * out,unsigned char * in,unsigned long long inlen)416 int crypto_hash( unsigned char *out, unsigned char *in, unsigned long long inlen )
417 {
418   return blake2b( out, in, NULL, BLAKE2B_OUTBYTES, inlen, 0 );
419 }
420 #endif
421 
422 #if defined(BLAKE2B_SELFTEST)
423 #include <string.h>
424 #include "blake2-kat.h"
main(int argc,char ** argv)425 int main( int argc, char **argv )
426 {
427   uint8_t key[BLAKE2B_KEYBYTES];
428   uint8_t buf[KAT_LENGTH];
429   size_t i;
430 
431   for( i = 0; i < BLAKE2B_KEYBYTES; ++i )
432     key[i] = ( uint8_t )i;
433 
434   for( i = 0; i < KAT_LENGTH; ++i )
435     buf[i] = ( uint8_t )i;
436 
437   for( i = 0; i < KAT_LENGTH; ++i )
438   {
439     uint8_t hash[BLAKE2B_OUTBYTES];
440     blake2b( hash, buf, key, BLAKE2B_OUTBYTES, i, BLAKE2B_KEYBYTES );
441 
442     if( 0 != memcmp( hash, blake2b_keyed_kat[i], BLAKE2B_OUTBYTES ) )
443     {
444       puts( "error" );
445       return -1;
446     }
447   }
448 
449   puts( "ok" );
450   return 0;
451 }
452 #endif
453 
454