xref: /freebsd/sys/contrib/openzfs/module/icp/algs/blake3/blake3_impl.h (revision 058ac3e8063366dafa634d9107642e12b038bf09)
1 /*
2  * CDDL HEADER START
3  *
4  * The contents of this file are subject to the terms of the
5  * Common Development and Distribution License (the "License").
6  * You may not use this file except in compliance with the License.
7  *
8  * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9  * or https://opensource.org/licenses/CDDL-1.0.
10  * See the License for the specific language governing permissions
11  * and limitations under the License.
12  *
13  * When distributing Covered Code, include this CDDL HEADER in each
14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15  * If applicable, add the following below this CDDL HEADER, with the
16  * fields enclosed by brackets "[]" replaced with your own identifying
17  * information: Portions Copyright [yyyy] [name of copyright owner]
18  *
19  * CDDL HEADER END
20  */
21 
22 /*
23  * Based on BLAKE3 v1.3.1, https://github.com/BLAKE3-team/BLAKE3
24  * Copyright (c) 2019-2020 Samuel Neves and Jack O'Connor
25  * Copyright (c) 2021-2022 Tino Reichardt <milky-zfs@mcmilk.de>
26  */
27 
28 #ifndef BLAKE3_IMPL_H
29 #define	BLAKE3_IMPL_H
30 
31 #ifdef	__cplusplus
32 extern "C" {
33 #endif
34 
35 #include <sys/types.h>
36 #include <sys/blake3.h>
37 #include <sys/simd.h>
38 
39 /*
40  * Methods used to define BLAKE3 assembler implementations
41  */
42 typedef void (*blake3_compress_in_place_f)(uint32_t cv[8],
43     const uint8_t block[BLAKE3_BLOCK_LEN],
44     uint8_t block_len, uint64_t counter,
45     uint8_t flags);
46 
47 typedef void (*blake3_compress_xof_f)(const uint32_t cv[8],
48     const uint8_t block[BLAKE3_BLOCK_LEN], uint8_t block_len,
49     uint64_t counter, uint8_t flags, uint8_t out[64]);
50 
51 typedef void (*blake3_hash_many_f)(const uint8_t * const *inputs,
52     size_t num_inputs, size_t blocks, const uint32_t key[8],
53     uint64_t counter, boolean_t increment_counter, uint8_t flags,
54     uint8_t flags_start, uint8_t flags_end, uint8_t *out);
55 
56 typedef boolean_t (*blake3_is_supported_f)(void);
57 
58 typedef struct blake3_impl_ops {
59 	blake3_compress_in_place_f compress_in_place;
60 	blake3_compress_xof_f compress_xof;
61 	blake3_hash_many_f hash_many;
62 	blake3_is_supported_f is_supported;
63 	int degree;
64 	const char *name;
65 } blake3_ops_t;
66 
67 /* Return selected BLAKE3 implementation ops */
68 extern const blake3_ops_t *blake3_impl_get_ops(void);
69 
70 extern const blake3_ops_t blake3_generic_impl;
71 
72 #if defined(__aarch64__) || \
73 	(defined(__x86_64) && defined(HAVE_SSE2)) || \
74 	(defined(__PPC64__) && defined(__LITTLE_ENDIAN__))
75 extern const blake3_ops_t blake3_sse2_impl;
76 #endif
77 
78 #if defined(__aarch64__) || \
79 	(defined(__x86_64) && defined(HAVE_SSE4_1)) || \
80 	(defined(__PPC64__) && defined(__LITTLE_ENDIAN__))
81 extern const blake3_ops_t blake3_sse41_impl;
82 #endif
83 
84 #if defined(__x86_64) && defined(HAVE_SSE4_1) && defined(HAVE_AVX2)
85 extern const blake3_ops_t blake3_avx2_impl;
86 #endif
87 
88 #if defined(__x86_64) && defined(HAVE_AVX512F) && defined(HAVE_AVX512VL)
89 extern const blake3_ops_t blake3_avx512_impl;
90 #endif
91 
92 #if defined(__x86_64)
93 #define	MAX_SIMD_DEGREE 16
94 #else
95 #define	MAX_SIMD_DEGREE 4
96 #endif
97 
98 #define	MAX_SIMD_DEGREE_OR_2	(MAX_SIMD_DEGREE > 2 ? MAX_SIMD_DEGREE : 2)
99 
100 static const uint32_t BLAKE3_IV[8] = {
101 	0x6A09E667UL, 0xBB67AE85UL, 0x3C6EF372UL, 0xA54FF53AUL,
102 	0x510E527FUL, 0x9B05688CUL, 0x1F83D9ABUL, 0x5BE0CD19UL};
103 
104 static const uint8_t BLAKE3_MSG_SCHEDULE[7][16] = {
105 	{0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15},
106 	{2, 6, 3, 10, 7, 0, 4, 13, 1, 11, 12, 5, 9, 14, 15, 8},
107 	{3, 4, 10, 12, 13, 2, 7, 14, 6, 5, 9, 0, 11, 15, 8, 1},
108 	{10, 7, 12, 9, 14, 3, 13, 15, 4, 0, 11, 2, 5, 8, 1, 6},
109 	{12, 13, 9, 11, 15, 10, 14, 8, 7, 2, 5, 3, 0, 1, 6, 4},
110 	{9, 14, 11, 5, 8, 12, 15, 1, 13, 3, 0, 10, 2, 6, 4, 7},
111 	{11, 15, 5, 0, 1, 9, 8, 6, 14, 10, 2, 12, 3, 4, 7, 13},
112 };
113 
114 /* Find index of the highest set bit */
115 static inline unsigned int highest_one(uint64_t x) {
116 #if defined(__GNUC__) || defined(__clang__)
117 	return (63 ^ __builtin_clzll(x));
118 #elif defined(_MSC_VER) && defined(IS_X86_64)
119 	unsigned long index;
120 	_BitScanReverse64(&index, x);
121 	return (index);
122 #elif defined(_MSC_VER) && defined(IS_X86_32)
123 	if (x >> 32) {
124 		unsigned long index;
125 		_BitScanReverse(&index, x >> 32);
126 		return (32 + index);
127 	} else {
128 		unsigned long index;
129 		_BitScanReverse(&index, x);
130 		return (index);
131 	}
132 #else
133 	unsigned int c = 0;
134 	if (x & 0xffffffff00000000ULL) { x >>= 32; c += 32; }
135 	if (x & 0x00000000ffff0000ULL) { x >>= 16; c += 16; }
136 	if (x & 0x000000000000ff00ULL) { x >>=  8; c +=  8; }
137 	if (x & 0x00000000000000f0ULL) { x >>=  4; c +=  4; }
138 	if (x & 0x000000000000000cULL) { x >>=  2; c +=  2; }
139 	if (x & 0x0000000000000002ULL) { c +=  1; }
140 	return (c);
141 #endif
142 }
143 
144 /* Count the number of 1 bits. */
145 static inline unsigned int popcnt(uint64_t x) {
146 	unsigned int count = 0;
147 
148 	while (x != 0) {
149 		count += 1;
150 		x &= x - 1;
151 	}
152 
153 	return (count);
154 }
155 
156 /*
157  * Largest power of two less than or equal to x.
158  * As a special case, returns 1 when x is 0.
159  */
160 static inline uint64_t round_down_to_power_of_2(uint64_t x) {
161 	return (1ULL << highest_one(x | 1));
162 }
163 
164 static inline uint32_t counter_low(uint64_t counter) {
165 	return ((uint32_t)counter);
166 }
167 
168 static inline uint32_t counter_high(uint64_t counter) {
169 	return ((uint32_t)(counter >> 32));
170 }
171 
172 static inline uint32_t load32(const void *src) {
173 	const uint8_t *p = (const uint8_t *)src;
174 	return ((uint32_t)(p[0]) << 0) | ((uint32_t)(p[1]) << 8) |
175 	    ((uint32_t)(p[2]) << 16) | ((uint32_t)(p[3]) << 24);
176 }
177 
178 static inline void load_key_words(const uint8_t key[BLAKE3_KEY_LEN],
179     uint32_t key_words[8]) {
180 	key_words[0] = load32(&key[0 * 4]);
181 	key_words[1] = load32(&key[1 * 4]);
182 	key_words[2] = load32(&key[2 * 4]);
183 	key_words[3] = load32(&key[3 * 4]);
184 	key_words[4] = load32(&key[4 * 4]);
185 	key_words[5] = load32(&key[5 * 4]);
186 	key_words[6] = load32(&key[6 * 4]);
187 	key_words[7] = load32(&key[7 * 4]);
188 }
189 
190 static inline void store32(void *dst, uint32_t w) {
191 	uint8_t *p = (uint8_t *)dst;
192 	p[0] = (uint8_t)(w >> 0);
193 	p[1] = (uint8_t)(w >> 8);
194 	p[2] = (uint8_t)(w >> 16);
195 	p[3] = (uint8_t)(w >> 24);
196 }
197 
198 static inline void store_cv_words(uint8_t bytes_out[32], uint32_t cv_words[8]) {
199 	store32(&bytes_out[0 * 4], cv_words[0]);
200 	store32(&bytes_out[1 * 4], cv_words[1]);
201 	store32(&bytes_out[2 * 4], cv_words[2]);
202 	store32(&bytes_out[3 * 4], cv_words[3]);
203 	store32(&bytes_out[4 * 4], cv_words[4]);
204 	store32(&bytes_out[5 * 4], cv_words[5]);
205 	store32(&bytes_out[6 * 4], cv_words[6]);
206 	store32(&bytes_out[7 * 4], cv_words[7]);
207 }
208 
209 #ifdef	__cplusplus
210 }
211 #endif
212 
213 #endif				/* BLAKE3_IMPL_H */
214