xref: /freebsd/sys/contrib/openzfs/module/icp/algs/blake3/blake3_impl.h (revision 61145dc2b94f12f6a47344fb9aac702321880e43)
1*61145dc2SMartin Matuska // SPDX-License-Identifier: CDDL-1.0
21f1e2261SMartin Matuska /*
31f1e2261SMartin Matuska  * CDDL HEADER START
41f1e2261SMartin Matuska  *
51f1e2261SMartin Matuska  * The contents of this file are subject to the terms of the
61f1e2261SMartin Matuska  * Common Development and Distribution License (the "License").
71f1e2261SMartin Matuska  * You may not use this file except in compliance with the License.
81f1e2261SMartin Matuska  *
91f1e2261SMartin Matuska  * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
10271171e0SMartin Matuska  * or https://opensource.org/licenses/CDDL-1.0.
111f1e2261SMartin Matuska  * See the License for the specific language governing permissions
121f1e2261SMartin Matuska  * and limitations under the License.
131f1e2261SMartin Matuska  *
141f1e2261SMartin Matuska  * When distributing Covered Code, include this CDDL HEADER in each
151f1e2261SMartin Matuska  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
161f1e2261SMartin Matuska  * If applicable, add the following below this CDDL HEADER, with the
171f1e2261SMartin Matuska  * fields enclosed by brackets "[]" replaced with your own identifying
181f1e2261SMartin Matuska  * information: Portions Copyright [yyyy] [name of copyright owner]
191f1e2261SMartin Matuska  *
201f1e2261SMartin Matuska  * CDDL HEADER END
211f1e2261SMartin Matuska  */
221f1e2261SMartin Matuska 
231f1e2261SMartin Matuska /*
241f1e2261SMartin Matuska  * Based on BLAKE3 v1.3.1, https://github.com/BLAKE3-team/BLAKE3
251f1e2261SMartin Matuska  * Copyright (c) 2019-2020 Samuel Neves and Jack O'Connor
261f1e2261SMartin Matuska  * Copyright (c) 2021-2022 Tino Reichardt <milky-zfs@mcmilk.de>
271f1e2261SMartin Matuska  */
281f1e2261SMartin Matuska 
291f1e2261SMartin Matuska #ifndef	BLAKE3_IMPL_H
301f1e2261SMartin Matuska #define	BLAKE3_IMPL_H
311f1e2261SMartin Matuska 
321f1e2261SMartin Matuska #ifdef	__cplusplus
331f1e2261SMartin Matuska extern "C" {
341f1e2261SMartin Matuska #endif
351f1e2261SMartin Matuska 
361f1e2261SMartin Matuska #include <sys/blake3.h>
371f1e2261SMartin Matuska #include <sys/simd.h>
3815f0b8c3SMartin Matuska #include <sys/asm_linkage.h>
391f1e2261SMartin Matuska 
401f1e2261SMartin Matuska /*
411f1e2261SMartin Matuska  * Methods used to define BLAKE3 assembler implementations
421f1e2261SMartin Matuska  */
431f1e2261SMartin Matuska typedef void (*blake3_compress_in_place_f)(uint32_t cv[8],
441f1e2261SMartin Matuska     const uint8_t block[BLAKE3_BLOCK_LEN],
451f1e2261SMartin Matuska     uint8_t block_len, uint64_t counter,
461f1e2261SMartin Matuska     uint8_t flags);
471f1e2261SMartin Matuska 
481f1e2261SMartin Matuska typedef void (*blake3_compress_xof_f)(const uint32_t cv[8],
491f1e2261SMartin Matuska     const uint8_t block[BLAKE3_BLOCK_LEN], uint8_t block_len,
501f1e2261SMartin Matuska     uint64_t counter, uint8_t flags, uint8_t out[64]);
511f1e2261SMartin Matuska 
521f1e2261SMartin Matuska typedef void (*blake3_hash_many_f)(const uint8_t * const *inputs,
531f1e2261SMartin Matuska     size_t num_inputs, size_t blocks, const uint32_t key[8],
541f1e2261SMartin Matuska     uint64_t counter, boolean_t increment_counter, uint8_t flags,
551f1e2261SMartin Matuska     uint8_t flags_start, uint8_t flags_end, uint8_t *out);
561f1e2261SMartin Matuska 
571f1e2261SMartin Matuska typedef boolean_t (*blake3_is_supported_f)(void);
581f1e2261SMartin Matuska 
592a58b312SMartin Matuska typedef struct {
601f1e2261SMartin Matuska 	blake3_compress_in_place_f compress_in_place;
611f1e2261SMartin Matuska 	blake3_compress_xof_f compress_xof;
621f1e2261SMartin Matuska 	blake3_hash_many_f hash_many;
631f1e2261SMartin Matuska 	blake3_is_supported_f is_supported;
641f1e2261SMartin Matuska 	int degree;
651f1e2261SMartin Matuska 	const char *name;
66c7046f76SMartin Matuska } blake3_ops_t;
671f1e2261SMartin Matuska 
682a58b312SMartin Matuska /* return selected BLAKE3 implementation ops */
692a58b312SMartin Matuska extern const blake3_ops_t *blake3_get_ops(void);
701f1e2261SMartin Matuska 
711f1e2261SMartin Matuska #if defined(__x86_64)
721f1e2261SMartin Matuska #define	MAX_SIMD_DEGREE 16
731f1e2261SMartin Matuska #else
741f1e2261SMartin Matuska #define	MAX_SIMD_DEGREE 4
751f1e2261SMartin Matuska #endif
761f1e2261SMartin Matuska 
771f1e2261SMartin Matuska #define	MAX_SIMD_DEGREE_OR_2	(MAX_SIMD_DEGREE > 2 ? MAX_SIMD_DEGREE : 2)
781f1e2261SMartin Matuska 
791f1e2261SMartin Matuska static const uint32_t BLAKE3_IV[8] = {
801f1e2261SMartin Matuska 	0x6A09E667UL, 0xBB67AE85UL, 0x3C6EF372UL, 0xA54FF53AUL,
811f1e2261SMartin Matuska 	0x510E527FUL, 0x9B05688CUL, 0x1F83D9ABUL, 0x5BE0CD19UL};
821f1e2261SMartin Matuska 
831f1e2261SMartin Matuska static const uint8_t BLAKE3_MSG_SCHEDULE[7][16] = {
841f1e2261SMartin Matuska 	{0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15},
851f1e2261SMartin Matuska 	{2, 6, 3, 10, 7, 0, 4, 13, 1, 11, 12, 5, 9, 14, 15, 8},
861f1e2261SMartin Matuska 	{3, 4, 10, 12, 13, 2, 7, 14, 6, 5, 9, 0, 11, 15, 8, 1},
871f1e2261SMartin Matuska 	{10, 7, 12, 9, 14, 3, 13, 15, 4, 0, 11, 2, 5, 8, 1, 6},
881f1e2261SMartin Matuska 	{12, 13, 9, 11, 15, 10, 14, 8, 7, 2, 5, 3, 0, 1, 6, 4},
891f1e2261SMartin Matuska 	{9, 14, 11, 5, 8, 12, 15, 1, 13, 3, 0, 10, 2, 6, 4, 7},
901f1e2261SMartin Matuska 	{11, 15, 5, 0, 1, 9, 8, 6, 14, 10, 2, 12, 3, 4, 7, 13},
911f1e2261SMartin Matuska };
921f1e2261SMartin Matuska 
931f1e2261SMartin Matuska /* Find index of the highest set bit */
highest_one(uint64_t x)941f1e2261SMartin Matuska static inline unsigned int highest_one(uint64_t x) {
951f1e2261SMartin Matuska #if defined(__GNUC__) || defined(__clang__)
961f1e2261SMartin Matuska 	return (63 ^ __builtin_clzll(x));
971f1e2261SMartin Matuska #elif defined(_MSC_VER) && defined(IS_X86_64)
981f1e2261SMartin Matuska 	unsigned long index;
991f1e2261SMartin Matuska 	_BitScanReverse64(&index, x);
1001f1e2261SMartin Matuska 	return (index);
1011f1e2261SMartin Matuska #elif defined(_MSC_VER) && defined(IS_X86_32)
1021f1e2261SMartin Matuska 	if (x >> 32) {
1031f1e2261SMartin Matuska 		unsigned long index;
1041f1e2261SMartin Matuska 		_BitScanReverse(&index, x >> 32);
1051f1e2261SMartin Matuska 		return (32 + index);
1061f1e2261SMartin Matuska 	} else {
1071f1e2261SMartin Matuska 		unsigned long index;
1081f1e2261SMartin Matuska 		_BitScanReverse(&index, x);
1091f1e2261SMartin Matuska 		return (index);
1101f1e2261SMartin Matuska 	}
1111f1e2261SMartin Matuska #else
1121f1e2261SMartin Matuska 	unsigned int c = 0;
1131f1e2261SMartin Matuska 	if (x & 0xffffffff00000000ULL) { x >>= 32; c += 32; }
1141f1e2261SMartin Matuska 	if (x & 0x00000000ffff0000ULL) { x >>= 16; c += 16; }
1151f1e2261SMartin Matuska 	if (x & 0x000000000000ff00ULL) { x >>=  8; c +=  8; }
1161f1e2261SMartin Matuska 	if (x & 0x00000000000000f0ULL) { x >>=  4; c +=  4; }
1171f1e2261SMartin Matuska 	if (x & 0x000000000000000cULL) { x >>=  2; c +=  2; }
1181f1e2261SMartin Matuska 	if (x & 0x0000000000000002ULL) { c +=  1; }
1191f1e2261SMartin Matuska 	return (c);
1201f1e2261SMartin Matuska #endif
1211f1e2261SMartin Matuska }
1221f1e2261SMartin Matuska 
1231f1e2261SMartin Matuska /* Count the number of 1 bits. */
popcnt(uint64_t x)1241f1e2261SMartin Matuska static inline unsigned int popcnt(uint64_t x) {
1251f1e2261SMartin Matuska 	unsigned int count = 0;
1261f1e2261SMartin Matuska 
1271f1e2261SMartin Matuska 	while (x != 0) {
1281f1e2261SMartin Matuska 		count += 1;
1291f1e2261SMartin Matuska 		x &= x - 1;
1301f1e2261SMartin Matuska 	}
1311f1e2261SMartin Matuska 
1321f1e2261SMartin Matuska 	return (count);
1331f1e2261SMartin Matuska }
1341f1e2261SMartin Matuska 
1351f1e2261SMartin Matuska /*
1361f1e2261SMartin Matuska  * Largest power of two less than or equal to x.
1371f1e2261SMartin Matuska  * As a special case, returns 1 when x is 0.
1381f1e2261SMartin Matuska  */
round_down_to_power_of_2(uint64_t x)1391f1e2261SMartin Matuska static inline uint64_t round_down_to_power_of_2(uint64_t x) {
1401f1e2261SMartin Matuska 	return (1ULL << highest_one(x | 1));
1411f1e2261SMartin Matuska }
1421f1e2261SMartin Matuska 
counter_low(uint64_t counter)1431f1e2261SMartin Matuska static inline uint32_t counter_low(uint64_t counter) {
1441f1e2261SMartin Matuska 	return ((uint32_t)counter);
1451f1e2261SMartin Matuska }
1461f1e2261SMartin Matuska 
counter_high(uint64_t counter)1471f1e2261SMartin Matuska static inline uint32_t counter_high(uint64_t counter) {
1481f1e2261SMartin Matuska 	return ((uint32_t)(counter >> 32));
1491f1e2261SMartin Matuska }
1501f1e2261SMartin Matuska 
load32(const void * src)1511f1e2261SMartin Matuska static inline uint32_t load32(const void *src) {
1521f1e2261SMartin Matuska 	const uint8_t *p = (const uint8_t *)src;
1531f1e2261SMartin Matuska 	return ((uint32_t)(p[0]) << 0) | ((uint32_t)(p[1]) << 8) |
1541f1e2261SMartin Matuska 	    ((uint32_t)(p[2]) << 16) | ((uint32_t)(p[3]) << 24);
1551f1e2261SMartin Matuska }
1561f1e2261SMartin Matuska 
load_key_words(const uint8_t key[BLAKE3_KEY_LEN],uint32_t key_words[8])1571f1e2261SMartin Matuska static inline void load_key_words(const uint8_t key[BLAKE3_KEY_LEN],
1581f1e2261SMartin Matuska     uint32_t key_words[8]) {
1591f1e2261SMartin Matuska 	key_words[0] = load32(&key[0 * 4]);
1601f1e2261SMartin Matuska 	key_words[1] = load32(&key[1 * 4]);
1611f1e2261SMartin Matuska 	key_words[2] = load32(&key[2 * 4]);
1621f1e2261SMartin Matuska 	key_words[3] = load32(&key[3 * 4]);
1631f1e2261SMartin Matuska 	key_words[4] = load32(&key[4 * 4]);
1641f1e2261SMartin Matuska 	key_words[5] = load32(&key[5 * 4]);
1651f1e2261SMartin Matuska 	key_words[6] = load32(&key[6 * 4]);
1661f1e2261SMartin Matuska 	key_words[7] = load32(&key[7 * 4]);
1671f1e2261SMartin Matuska }
1681f1e2261SMartin Matuska 
store32(void * dst,uint32_t w)1691f1e2261SMartin Matuska static inline void store32(void *dst, uint32_t w) {
1701f1e2261SMartin Matuska 	uint8_t *p = (uint8_t *)dst;
1711f1e2261SMartin Matuska 	p[0] = (uint8_t)(w >> 0);
1721f1e2261SMartin Matuska 	p[1] = (uint8_t)(w >> 8);
1731f1e2261SMartin Matuska 	p[2] = (uint8_t)(w >> 16);
1741f1e2261SMartin Matuska 	p[3] = (uint8_t)(w >> 24);
1751f1e2261SMartin Matuska }
1761f1e2261SMartin Matuska 
store_cv_words(uint8_t bytes_out[32],uint32_t cv_words[8])1771f1e2261SMartin Matuska static inline void store_cv_words(uint8_t bytes_out[32], uint32_t cv_words[8]) {
1781f1e2261SMartin Matuska 	store32(&bytes_out[0 * 4], cv_words[0]);
1791f1e2261SMartin Matuska 	store32(&bytes_out[1 * 4], cv_words[1]);
1801f1e2261SMartin Matuska 	store32(&bytes_out[2 * 4], cv_words[2]);
1811f1e2261SMartin Matuska 	store32(&bytes_out[3 * 4], cv_words[3]);
1821f1e2261SMartin Matuska 	store32(&bytes_out[4 * 4], cv_words[4]);
1831f1e2261SMartin Matuska 	store32(&bytes_out[5 * 4], cv_words[5]);
1841f1e2261SMartin Matuska 	store32(&bytes_out[6 * 4], cv_words[6]);
1851f1e2261SMartin Matuska 	store32(&bytes_out[7 * 4], cv_words[7]);
1861f1e2261SMartin Matuska }
1871f1e2261SMartin Matuska 
1881f1e2261SMartin Matuska #ifdef	__cplusplus
1891f1e2261SMartin Matuska }
1901f1e2261SMartin Matuska #endif
1911f1e2261SMartin Matuska 
1921f1e2261SMartin Matuska #endif				/* BLAKE3_IMPL_H */
193