diff options
Diffstat (limited to 'blake2b-simd.cpp')
| -rw-r--r-- | blake2b-simd.cpp | 1235 |
1 files changed, 1235 insertions, 0 deletions
diff --git a/blake2b-simd.cpp b/blake2b-simd.cpp new file mode 100644 index 00000000..5a7f65d3 --- /dev/null +++ b/blake2b-simd.cpp @@ -0,0 +1,1235 @@ +
+// blake2-simd.cpp - written and placed in the public domain by
+// Jeffrey Walton, Uri Blumenthal and Marcel Raad.
+//
+// This source file uses intrinsics to gain access to ARMv7a/ARMv8a
+// NEON, Power8 and SSE4.1 instructions. A separate source file is
+// needed because additional CXXFLAGS are required to enable the
+// appropriate instructions sets in some build configurations.
+
+#include "pch.h"
+#include "config.h"
+#include "misc.h"
+#include "blake2.h"
+
+// Uncomment for benchmarking C++ against SSE2 or NEON.
+// Do so in both blake2.cpp and blake2-simd.cpp.
+// #undef CRYPTOPP_SSE41_AVAILABLE
+// #undef CRYPTOPP_ARM_NEON_AVAILABLE
+// #undef CRYPTOPP_POWER8_AVAILABLE
+
+// Disable NEON/ASIMD for Cortex-A53 and A57. The shifts are too slow and C/C++ is about
+// 3 cpb faster than NEON/ASIMD. Also see http://github.com/weidai11/cryptopp/issues/367.
+#if (defined(__aarch32__) || defined(__aarch64__)) && defined(CRYPTOPP_SLOW_ARMV8_SHIFT)
+# undef CRYPTOPP_ARM_NEON_AVAILABLE
+#endif
+
+#if (CRYPTOPP_SSE41_AVAILABLE)
+# include <emmintrin.h>
+# include <tmmintrin.h>
+# include <smmintrin.h>
+#endif
+
+#if (CRYPTOPP_ARM_NEON_AVAILABLE)
+# include <arm_neon.h>
+#endif
+
+// Can't use CRYPTOPP_ARM_XXX_AVAILABLE because too many
+// compilers don't follow ACLE conventions for the include.
+#if (CRYPTOPP_ARM_ACLE_AVAILABLE)
+# include <stdint.h>
+# include <arm_acle.h>
+#endif
+
+#if defined(CRYPTOPP_POWER8_AVAILABLE)
+# include "ppc-simd.h"
+#endif
+
+ANONYMOUS_NAMESPACE_BEGIN
+
+using CryptoPP::word32;
+using CryptoPP::word64;
+
+#if (CRYPTOPP_SSE41_AVAILABLE || CRYPTOPP_ARM_NEON_AVAILABLE || CRYPTOPP_POWER8_AVAILABLE)
+
+CRYPTOPP_ALIGN_DATA(16)
+const word64 BLAKE2B_IV[8] = {
+ W64LIT(0x6a09e667f3bcc908), W64LIT(0xbb67ae8584caa73b),
+ W64LIT(0x3c6ef372fe94f82b), W64LIT(0xa54ff53a5f1d36f1),
+ W64LIT(0x510e527fade682d1), W64LIT(0x9b05688c2b3e6c1f),
+ W64LIT(0x1f83d9abfb41bd6b), W64LIT(0x5be0cd19137e2179)
+};
+
+#endif // CRYPTOPP_SSE41_AVAILABLE || CRYPTOPP_ARM_NEON_AVAILABLE
+
+ANONYMOUS_NAMESPACE_END
+
+NAMESPACE_BEGIN(CryptoPP)
+
+#if CRYPTOPP_SSE41_AVAILABLE
+
+#define LOADU(p) _mm_loadu_si128( (const __m128i *)(const void*)(p) )
+#define STOREU(p,r) _mm_storeu_si128((__m128i *)(void*)(p), r)
+#define TOF(reg) _mm_castsi128_ps((reg))
+#define TOI(reg) _mm_castps_si128((reg))
+
+void BLAKE2_Compress64_SSE4(const byte* input, BLAKE2_State<word64, true>& state)
+{
+ #define BLAKE2B_LOAD_MSG_0_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m0, m1); \
+ b1 = _mm_unpacklo_epi64(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_2(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m0, m1); \
+ b1 = _mm_unpackhi_epi64(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_3(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m4, m5); \
+ b1 = _mm_unpacklo_epi64(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_4(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m4, m5); \
+ b1 = _mm_unpackhi_epi64(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m7, m2); \
+ b1 = _mm_unpackhi_epi64(m4, m6); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_2(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m5, m4); \
+ b1 = _mm_alignr_epi8(m3, m7, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_3(b0, b1) \
+ do { \
+ b0 = _mm_shuffle_epi32(m0, _MM_SHUFFLE(1,0,3,2)); \
+ b1 = _mm_unpackhi_epi64(m5, m2); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_4(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m6, m1); \
+ b1 = _mm_unpackhi_epi64(m3, m1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_1(b0, b1) \
+ do { \
+ b0 = _mm_alignr_epi8(m6, m5, 8); \
+ b1 = _mm_unpackhi_epi64(m2, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_2(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m4, m0); \
+ b1 = _mm_blend_epi16(m1, m6, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_3(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m5, m1, 0xF0); \
+ b1 = _mm_unpackhi_epi64(m3, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_4(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m7, m3); \
+ b1 = _mm_alignr_epi8(m2, m0, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_1(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m3, m1); \
+ b1 = _mm_unpackhi_epi64(m6, m5); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_2(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m4, m0); \
+ b1 = _mm_unpacklo_epi64(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_3(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m1, m2, 0xF0); \
+ b1 = _mm_blend_epi16(m2, m7, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_4(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m3, m5); \
+ b1 = _mm_unpacklo_epi64(m0, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_1(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m4, m2); \
+ b1 = _mm_unpacklo_epi64(m1, m5); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_2(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m0, m3, 0xF0); \
+ b1 = _mm_blend_epi16(m2, m7, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_3(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m7, m5, 0xF0); \
+ b1 = _mm_blend_epi16(m3, m1, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_4(b0, b1) \
+ do { \
+ b0 = _mm_alignr_epi8(m6, m0, 8); \
+ b1 = _mm_blend_epi16(m4, m6, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m1, m3); \
+ b1 = _mm_unpacklo_epi64(m0, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_2(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m6, m5); \
+ b1 = _mm_unpackhi_epi64(m5, m1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_3(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m2, m3, 0xF0); \
+ b1 = _mm_unpackhi_epi64(m7, m0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_4(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m6, m2); \
+ b1 = _mm_blend_epi16(m7, m4, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_1(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m6, m0, 0xF0); \
+ b1 = _mm_unpacklo_epi64(m7, m2); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_2(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m2, m7); \
+ b1 = _mm_alignr_epi8(m5, m6, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_3(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m0, m3); \
+ b1 = _mm_shuffle_epi32(m4, _MM_SHUFFLE(1,0,3,2)); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_4(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m3, m1); \
+ b1 = _mm_blend_epi16(m1, m5, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_1(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m6, m3); \
+ b1 = _mm_blend_epi16(m6, m1, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_2(b0, b1) \
+ do { \
+ b0 = _mm_alignr_epi8(m7, m5, 8); \
+ b1 = _mm_unpackhi_epi64(m0, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_3(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m2, m7); \
+ b1 = _mm_unpacklo_epi64(m4, m1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_4(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m0, m2); \
+ b1 = _mm_unpacklo_epi64(m3, m5); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m3, m7); \
+ b1 = _mm_alignr_epi8(m0, m5, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_2(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m7, m4); \
+ b1 = _mm_alignr_epi8(m4, m1, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_3(b0, b1) \
+ do { \
+ b0 = m6; \
+ b1 = _mm_alignr_epi8(m5, m0, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_4(b0, b1) \
+ do { \
+ b0 = _mm_blend_epi16(m1, m3, 0xF0); \
+ b1 = m2; \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m5, m4); \
+ b1 = _mm_unpackhi_epi64(m3, m0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_2(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m1, m2); \
+ b1 = _mm_blend_epi16(m3, m2, 0xF0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_3(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m7, m4); \
+ b1 = _mm_unpackhi_epi64(m1, m6); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_4(b0, b1) \
+ do { \
+ b0 = _mm_alignr_epi8(m7, m5, 8); \
+ b1 = _mm_unpacklo_epi64(m6, m0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m0, m1); \
+ b1 = _mm_unpacklo_epi64(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_2(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m0, m1); \
+ b1 = _mm_unpackhi_epi64(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_3(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m4, m5); \
+ b1 = _mm_unpacklo_epi64(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_4(b0, b1) \
+ do { \
+ b0 = _mm_unpackhi_epi64(m4, m5); \
+ b1 = _mm_unpackhi_epi64(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_1(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m7, m2); \
+ b1 = _mm_unpackhi_epi64(m4, m6); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_2(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m5, m4); \
+ b1 = _mm_alignr_epi8(m3, m7, 8); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_3(b0, b1) \
+ do { \
+ b0 = _mm_shuffle_epi32(m0, _MM_SHUFFLE(1,0,3,2)); \
+ b1 = _mm_unpackhi_epi64(m5, m2); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_4(b0, b1) \
+ do { \
+ b0 = _mm_unpacklo_epi64(m6, m1); \
+ b1 = _mm_unpackhi_epi64(m3, m1); \
+ } while(0)
+
+#ifdef __XOP__
+# define MM_ROTI_EPI64(r, c) \
+ _mm_roti_epi64(r, c)
+#else
+# define MM_ROTI_EPI64(x, c) \
+ (-(c) == 32) ? _mm_shuffle_epi32((x), _MM_SHUFFLE(2,3,0,1)) \
+ : (-(c) == 24) ? _mm_shuffle_epi8((x), r24) \
+ : (-(c) == 16) ? _mm_shuffle_epi8((x), r16) \
+ : (-(c) == 63) ? _mm_xor_si128(_mm_srli_epi64((x), -(c)), _mm_add_epi64((x), (x))) \
+ : _mm_xor_si128(_mm_srli_epi64((x), -(c)), _mm_slli_epi64((x), 64-(-(c))))
+#endif
+
+#define BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1) \
+ row1l = _mm_add_epi64(_mm_add_epi64(row1l, b0), row2l); \
+ row1h = _mm_add_epi64(_mm_add_epi64(row1h, b1), row2h); \
+ \
+ row4l = _mm_xor_si128(row4l, row1l); \
+ row4h = _mm_xor_si128(row4h, row1h); \
+ \
+ row4l = MM_ROTI_EPI64(row4l, -32); \
+ row4h = MM_ROTI_EPI64(row4h, -32); \
+ \
+ row3l = _mm_add_epi64(row3l, row4l); \
+ row3h = _mm_add_epi64(row3h, row4h); \
+ \
+ row2l = _mm_xor_si128(row2l, row3l); \
+ row2h = _mm_xor_si128(row2h, row3h); \
+ \
+ row2l = MM_ROTI_EPI64(row2l, -24); \
+ row2h = MM_ROTI_EPI64(row2h, -24);
+
+#define BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1) \
+ row1l = _mm_add_epi64(_mm_add_epi64(row1l, b0), row2l); \
+ row1h = _mm_add_epi64(_mm_add_epi64(row1h, b1), row2h); \
+ \
+ row4l = _mm_xor_si128(row4l, row1l); \
+ row4h = _mm_xor_si128(row4h, row1h); \
+ \
+ row4l = MM_ROTI_EPI64(row4l, -16); \
+ row4h = MM_ROTI_EPI64(row4h, -16); \
+ \
+ row3l = _mm_add_epi64(row3l, row4l); \
+ row3h = _mm_add_epi64(row3h, row4h); \
+ \
+ row2l = _mm_xor_si128(row2l, row3l); \
+ row2h = _mm_xor_si128(row2h, row3h); \
+ \
+ row2l = MM_ROTI_EPI64(row2l, -63); \
+ row2h = MM_ROTI_EPI64(row2h, -63); \
+
+#define BLAKE2B_DIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h) \
+ t0 = row4l;\
+ t1 = row2l;\
+ row4l = row3l;\
+ row3l = row3h;\
+ row3h = row4l;\
+ row4l = _mm_unpackhi_epi64(row4h, _mm_unpacklo_epi64(t0, t0)); \
+ row4h = _mm_unpackhi_epi64(t0, _mm_unpacklo_epi64(row4h, row4h)); \
+ row2l = _mm_unpackhi_epi64(row2l, _mm_unpacklo_epi64(row2h, row2h)); \
+ row2h = _mm_unpackhi_epi64(row2h, _mm_unpacklo_epi64(t1, t1))
+
+#define BLAKE2B_UNDIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h) \
+ t0 = row3l;\
+ row3l = row3h;\
+ row3h = t0;\
+ t0 = row2l;\
+ t1 = row4l;\
+ row2l = _mm_unpackhi_epi64(row2h, _mm_unpacklo_epi64(row2l, row2l)); \
+ row2h = _mm_unpackhi_epi64(t0, _mm_unpacklo_epi64(row2h, row2h)); \
+ row4l = _mm_unpackhi_epi64(row4l, _mm_unpacklo_epi64(row4h, row4h)); \
+ row4h = _mm_unpackhi_epi64(row4h, _mm_unpacklo_epi64(t1, t1))
+
+#define BLAKE2B_ROUND(r) \
+ BLAKE2B_LOAD_MSG_ ##r ##_1(b0, b1); \
+ BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_LOAD_MSG_ ##r ##_2(b0, b1); \
+ BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_DIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h); \
+ BLAKE2B_LOAD_MSG_ ##r ##_3(b0, b1); \
+ BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_LOAD_MSG_ ##r ##_4(b0, b1); \
+ BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_UNDIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h);
+
+ __m128i row1l, row1h;
+ __m128i row2l, row2h;
+ __m128i row3l, row3h;
+ __m128i row4l, row4h;
+ __m128i b0, b1;
+ __m128i t0, t1;
+
+ const __m128i r16 = _mm_setr_epi8( 2, 3, 4, 5, 6, 7, 0, 1, 10, 11, 12, 13, 14, 15, 8, 9 );
+ const __m128i r24 = _mm_setr_epi8( 3, 4, 5, 6, 7, 0, 1, 2, 11, 12, 13, 14, 15, 8, 9, 10 );
+
+ const __m128i m0 = LOADU( input + 00 );
+ const __m128i m1 = LOADU( input + 16 );
+ const __m128i m2 = LOADU( input + 32 );
+ const __m128i m3 = LOADU( input + 48 );
+ const __m128i m4 = LOADU( input + 64 );
+ const __m128i m5 = LOADU( input + 80 );
+ const __m128i m6 = LOADU( input + 96 );
+ const __m128i m7 = LOADU( input + 112 );
+
+ row1l = LOADU( &state.h[0] );
+ row1h = LOADU( &state.h[2] );
+ row2l = LOADU( &state.h[4] );
+ row2h = LOADU( &state.h[6] );
+ row3l = LOADU( &BLAKE2B_IV[0] );
+ row3h = LOADU( &BLAKE2B_IV[2] );
+ row4l = _mm_xor_si128( LOADU( &BLAKE2B_IV[4] ), LOADU( &state.t[0] ) );
+ row4h = _mm_xor_si128( LOADU( &BLAKE2B_IV[6] ), LOADU( &state.f[0] ) );
+
+ BLAKE2B_ROUND( 0 );
+ BLAKE2B_ROUND( 1 );
+ BLAKE2B_ROUND( 2 );
+ BLAKE2B_ROUND( 3 );
+ BLAKE2B_ROUND( 4 );
+ BLAKE2B_ROUND( 5 );
+ BLAKE2B_ROUND( 6 );
+ BLAKE2B_ROUND( 7 );
+ BLAKE2B_ROUND( 8 );
+ BLAKE2B_ROUND( 9 );
+ BLAKE2B_ROUND( 10 );
+ BLAKE2B_ROUND( 11 );
+
+ row1l = _mm_xor_si128( row3l, row1l );
+ row1h = _mm_xor_si128( row3h, row1h );
+ STOREU( &state.h[0], _mm_xor_si128( LOADU( &state.h[0] ), row1l ) );
+ STOREU( &state.h[2], _mm_xor_si128( LOADU( &state.h[2] ), row1h ) );
+ row2l = _mm_xor_si128( row4l, row2l );
+ row2h = _mm_xor_si128( row4h, row2h );
+ STOREU( &state.h[4], _mm_xor_si128( LOADU( &state.h[4] ), row2l ) );
+ STOREU( &state.h[6], _mm_xor_si128( LOADU( &state.h[6] ), row2h ) );
+}
+#endif // CRYPTOPP_SSE41_AVAILABLE
+
+#if CRYPTOPP_ARM_NEON_AVAILABLE
+void BLAKE2_Compress64_NEON(const byte* input, BLAKE2_State<word64, true>& state)
+{
+ #define BLAKE2B_LOAD_MSG_0_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m0), vget_low_u64(m1)); b1 = vcombine_u64(vget_low_u64(m2), vget_low_u64(m3)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m0), vget_high_u64(m1)); b1 = vcombine_u64(vget_high_u64(m2), vget_high_u64(m3)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m4), vget_low_u64(m5)); b1 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m4), vget_high_u64(m5)); b1 = vcombine_u64(vget_high_u64(m6), vget_high_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m7), vget_low_u64(m2)); b1 = vcombine_u64(vget_high_u64(m4), vget_high_u64(m6)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m5), vget_low_u64(m4)); b1 = vextq_u64(m7, m3, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_3(b0, b1) \
+ do { b0 = vextq_u64(m0, m0, 1); b1 = vcombine_u64(vget_high_u64(m5), vget_high_u64(m2)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m1)); b1 = vcombine_u64(vget_high_u64(m3), vget_high_u64(m1)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_1(b0, b1) \
+ do { b0 = vextq_u64(m5, m6, 1); b1 = vcombine_u64(vget_high_u64(m2), vget_high_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m4), vget_low_u64(m0)); b1 = vcombine_u64(vget_low_u64(m1), vget_high_u64(m6)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m5), vget_high_u64(m1)); b1 = vcombine_u64(vget_high_u64(m3), vget_high_u64(m4)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m7), vget_low_u64(m3)); b1 = vextq_u64(m0, m2, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m3), vget_high_u64(m1)); b1 = vcombine_u64(vget_high_u64(m6), vget_high_u64(m5)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m4), vget_high_u64(m0)); b1 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m1), vget_high_u64(m2)); b1 = vcombine_u64(vget_low_u64(m2), vget_high_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m3), vget_low_u64(m5)); b1 = vcombine_u64(vget_low_u64(m0), vget_low_u64(m4)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m4), vget_high_u64(m2)); b1 = vcombine_u64(vget_low_u64(m1), vget_low_u64(m5)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m0), vget_high_u64(m3)); b1 = vcombine_u64(vget_low_u64(m2), vget_high_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m7), vget_high_u64(m5)); b1 = vcombine_u64(vget_low_u64(m3), vget_high_u64(m1)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_4(b0, b1) \
+ do { b0 = vextq_u64(m0, m6, 1); b1 = vcombine_u64(vget_low_u64(m4), vget_high_u64(m6)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m1), vget_low_u64(m3)); b1 = vcombine_u64(vget_low_u64(m0), vget_low_u64(m4)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m5)); b1 = vcombine_u64(vget_high_u64(m5), vget_high_u64(m1)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m2), vget_high_u64(m3)); b1 = vcombine_u64(vget_high_u64(m7), vget_high_u64(m0)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m6), vget_high_u64(m2)); b1 = vcombine_u64(vget_low_u64(m7), vget_high_u64(m4)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m6), vget_high_u64(m0)); b1 = vcombine_u64(vget_low_u64(m7), vget_low_u64(m2)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m2), vget_high_u64(m7)); b1 = vextq_u64(m6, m5, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m0), vget_low_u64(m3)); b1 = vextq_u64(m4, m4, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m3), vget_high_u64(m1)); b1 = vcombine_u64(vget_low_u64(m1), vget_high_u64(m5)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m6), vget_high_u64(m3)); b1 = vcombine_u64(vget_low_u64(m6), vget_high_u64(m1)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_2(b0, b1) \
+ do { b0 = vextq_u64(m5, m7, 1); b1 = vcombine_u64(vget_high_u64(m0), vget_high_u64(m4)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m2), vget_high_u64(m7)); b1 = vcombine_u64(vget_low_u64(m4), vget_low_u64(m1)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m0), vget_low_u64(m2)); b1 = vcombine_u64(vget_low_u64(m3), vget_low_u64(m5)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m3), vget_low_u64(m7)); b1 = vextq_u64(m5, m0, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m7), vget_high_u64(m4)); b1 = vextq_u64(m1, m4, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_3(b0, b1) \
+ do { b0 = m6; b1 = vextq_u64(m0, m5, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m1), vget_high_u64(m3)); b1 = m2; } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m5), vget_low_u64(m4)); b1 = vcombine_u64(vget_high_u64(m3), vget_high_u64(m0)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m1), vget_low_u64(m2)); b1 = vcombine_u64(vget_low_u64(m3), vget_high_u64(m2)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m7), vget_high_u64(m4)); b1 = vcombine_u64(vget_high_u64(m1), vget_high_u64(m6)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_4(b0, b1) \
+ do { b0 = vextq_u64(m5, m7, 1); b1 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m0)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m0), vget_low_u64(m1)); b1 = vcombine_u64(vget_low_u64(m2), vget_low_u64(m3)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m0), vget_high_u64(m1)); b1 = vcombine_u64(vget_high_u64(m2), vget_high_u64(m3)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_3(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m4), vget_low_u64(m5)); b1 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_high_u64(m4), vget_high_u64(m5)); b1 = vcombine_u64(vget_high_u64(m6), vget_high_u64(m7)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_1(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m7), vget_low_u64(m2)); b1 = vcombine_u64(vget_high_u64(m4), vget_high_u64(m6)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_2(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m5), vget_low_u64(m4)); b1 = vextq_u64(m7, m3, 1); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_3(b0, b1) \
+ do { b0 = vextq_u64(m0, m0, 1); b1 = vcombine_u64(vget_high_u64(m5), vget_high_u64(m2)); } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_4(b0, b1) \
+ do { b0 = vcombine_u64(vget_low_u64(m6), vget_low_u64(m1)); b1 = vcombine_u64(vget_high_u64(m3), vget_high_u64(m1)); } while(0)
+
+ #define vrorq_n_u64_32(x) vreinterpretq_u64_u32(vrev64q_u32(vreinterpretq_u32_u64((x))))
+
+ #define vrorq_n_u64_24(x) vcombine_u64(\
+ vreinterpret_u64_u8(vext_u8(vreinterpret_u8_u64(vget_low_u64(x)), vreinterpret_u8_u64(vget_low_u64(x)), 3)), \
+ vreinterpret_u64_u8(vext_u8(vreinterpret_u8_u64(vget_high_u64(x)), vreinterpret_u8_u64(vget_high_u64(x)), 3)))
+
+ #define vrorq_n_u64_16(x) vcombine_u64(\
+ vreinterpret_u64_u8(vext_u8(vreinterpret_u8_u64(vget_low_u64(x)), vreinterpret_u8_u64(vget_low_u64(x)), 2)), \
+ vreinterpret_u64_u8(vext_u8(vreinterpret_u8_u64(vget_high_u64(x)), vreinterpret_u8_u64(vget_high_u64(x)), 2)))
+
+ #define vrorq_n_u64_63(x) veorq_u64(vaddq_u64(x, x), vshrq_n_u64(x, 63))
+
+ #define BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1) \
+ do { \
+ row1l = vaddq_u64(vaddq_u64(row1l, b0), row2l); \
+ row1h = vaddq_u64(vaddq_u64(row1h, b1), row2h); \
+ row4l = veorq_u64(row4l, row1l); row4h = veorq_u64(row4h, row1h); \
+ row4l = vrorq_n_u64_32(row4l); row4h = vrorq_n_u64_32(row4h); \
+ row3l = vaddq_u64(row3l, row4l); row3h = vaddq_u64(row3h, row4h); \
+ row2l = veorq_u64(row2l, row3l); row2h = veorq_u64(row2h, row3h); \
+ row2l = vrorq_n_u64_24(row2l); row2h = vrorq_n_u64_24(row2h); \
+ } while(0)
+
+ #define BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1) \
+ do { \
+ row1l = vaddq_u64(vaddq_u64(row1l, b0), row2l); \
+ row1h = vaddq_u64(vaddq_u64(row1h, b1), row2h); \
+ row4l = veorq_u64(row4l, row1l); row4h = veorq_u64(row4h, row1h); \
+ row4l = vrorq_n_u64_16(row4l); row4h = vrorq_n_u64_16(row4h); \
+ row3l = vaddq_u64(row3l, row4l); row3h = vaddq_u64(row3h, row4h); \
+ row2l = veorq_u64(row2l, row3l); row2h = veorq_u64(row2h, row3h); \
+ row2l = vrorq_n_u64_63(row2l); row2h = vrorq_n_u64_63(row2h); \
+ } while(0)
+
+ #define BLAKE2B_DIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h) \
+ do { \
+ uint64x2_t t0 = vextq_u64(row2l, row2h, 1); \
+ uint64x2_t t1 = vextq_u64(row2h, row2l, 1); \
+ row2l = t0; row2h = t1; t0 = row3l; row3l = row3h; row3h = t0; \
+ t0 = vextq_u64(row4h, row4l, 1); t1 = vextq_u64(row4l, row4h, 1); \
+ row4l = t0; row4h = t1; \
+ } while(0)
+
+ #define BLAKE2B_UNDIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h) \
+ do { \
+ uint64x2_t t0 = vextq_u64(row2h, row2l, 1); \
+ uint64x2_t t1 = vextq_u64(row2l, row2h, 1); \
+ row2l = t0; row2h = t1; t0 = row3l; row3l = row3h; row3h = t0; \
+ t0 = vextq_u64(row4l, row4h, 1); t1 = vextq_u64(row4h, row4l, 1); \
+ row4l = t0; row4h = t1; \
+ } while(0)
+
+ #define BLAKE2B_ROUND(r) \
+ do { \
+ uint64x2_t b0, b1; \
+ BLAKE2B_LOAD_MSG_ ##r ##_1(b0, b1); \
+ BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_LOAD_MSG_ ##r ##_2(b0, b1); \
+ BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_DIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h); \
+ BLAKE2B_LOAD_MSG_ ##r ##_3(b0, b1); \
+ BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_LOAD_MSG_ ##r ##_4(b0, b1); \
+ BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_UNDIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h); \
+ } while(0)
+
+ const uint64x2_t m0 = vreinterpretq_u64_u8(vld1q_u8(input + 00));
+ const uint64x2_t m1 = vreinterpretq_u64_u8(vld1q_u8(input + 16));
+ const uint64x2_t m2 = vreinterpretq_u64_u8(vld1q_u8(input + 32));
+ const uint64x2_t m3 = vreinterpretq_u64_u8(vld1q_u8(input + 48));
+ const uint64x2_t m4 = vreinterpretq_u64_u8(vld1q_u8(input + 64));
+ const uint64x2_t m5 = vreinterpretq_u64_u8(vld1q_u8(input + 80));
+ const uint64x2_t m6 = vreinterpretq_u64_u8(vld1q_u8(input + 96));
+ const uint64x2_t m7 = vreinterpretq_u64_u8(vld1q_u8(input + 112));
+
+ uint64x2_t row1l, row1h, row2l, row2h;
+ uint64x2_t row3l, row3h, row4l, row4h;
+
+ const uint64x2_t h0 = row1l = vld1q_u64(&state.h[0]);
+ const uint64x2_t h1 = row1h = vld1q_u64(&state.h[2]);
+ const uint64x2_t h2 = row2l = vld1q_u64(&state.h[4]);
+ const uint64x2_t h3 = row2h = vld1q_u64(&state.h[6]);
+
+ row3l = vld1q_u64(&BLAKE2B_IV[0]);
+ row3h = vld1q_u64(&BLAKE2B_IV[2]);
+ row4l = veorq_u64(vld1q_u64(&BLAKE2B_IV[4]), vld1q_u64(&state.t[0]));
+ row4h = veorq_u64(vld1q_u64(&BLAKE2B_IV[6]), vld1q_u64(&state.f[0]));
+
+ BLAKE2B_ROUND(0);
+ BLAKE2B_ROUND(1);
+ BLAKE2B_ROUND(2);
+ BLAKE2B_ROUND(3);
+ BLAKE2B_ROUND(4);
+ BLAKE2B_ROUND(5);
+ BLAKE2B_ROUND(6);
+ BLAKE2B_ROUND(7);
+ BLAKE2B_ROUND(8);
+ BLAKE2B_ROUND(9);
+ BLAKE2B_ROUND(10);
+ BLAKE2B_ROUND(11);
+
+ vst1q_u64(&state.h[0], veorq_u64(h0, veorq_u64(row1l, row3l)));
+ vst1q_u64(&state.h[2], veorq_u64(h1, veorq_u64(row1h, row3h)));
+ vst1q_u64(&state.h[4], veorq_u64(h2, veorq_u64(row2l, row4l)));
+ vst1q_u64(&state.h[6], veorq_u64(h3, veorq_u64(row2h, row4h)));
+}
+#endif // CRYPTOPP_ARM_NEON_AVAILABLE
+
+#if (CRYPTOPP_POWER8_AVAILABLE)
+
+inline uint64x2_p VectorLoad64(const void* p)
+{
+#if defined(__xlc__) || defined(__xlC__) || defined(__clang__)
+ return (uint64x2_p)vec_xl(0, (uint8_t*)p);
+#else
+ return (uint64x2_p)vec_vsx_ld(0, (uint8_t*)p);
+#endif
+}
+
+inline uint64x2_p VectorLoad64LE(const void* p)
+{
+#if __BIG_ENDIAN__
+ const uint8x16_p m = {7,6,5,4, 3,2,1,0, 15,14,13,12, 11,10,9,8};
+ const uint64x2_p v = VectorLoad64(p);
+ return vec_perm(v, v, m);
+#else
+ return VectorLoad64(p);
+#endif
+}
+
+inline void VectorStore64(void* p, const uint64x2_p x)
+{
+#if defined(__xlc__) || defined(__xlC__) || defined(__clang__)
+ vec_xst((uint8x16_p)x,0,(uint8_t*)p);
+#else
+ vec_vsx_st((uint8x16_p)x,0,(uint8_t*)p);
+#endif
+}
+
+inline void VectorStore64LE(void* p, const uint64x2_p x)
+{
+#if __BIG_ENDIAN__
+ const uint8x16_p m = {7,6,5,4, 3,2,1,0, 15,14,13,12, 11,10,9,8};
+ VectorStore64(p, vec_perm(x, x, m));
+#else
+ VectorStore64(p, x);
+#endif
+}
+
+template <unsigned int C>
+inline uint64x2_p VectorShiftLeftOctet(const uint64x2_p a, const uint64x2_p b)
+{
+#if __BIG_ENDIAN__
+ return (uint64x2_p)vec_sld((uint8x16_p)a, (uint8x16_p)b, C);
+#else
+ return (uint64x2_p)vec_sld((uint8x16_p)b, (uint8x16_p)a, 16-C);
+#endif
+}
+
+#define vec_ext(a,b,c) VectorShiftLeftOctet<c*8>(a, b)
+
+// vec_mergeh(a,b) is equivalent to vec_perm(a,b,HH_MASK); and
+// vec_mergel(a,b) is equivalent vec_perm(a,b,LL_MASK). Benchmarks
+// show vec_mergeh and vec_mergel is faster on little-endian
+// machines by 0.4 cpb. Benchmarks show vec_perm is faster on
+// big-endian machines by 1.5 cpb. The code that uses
+// vec_mergeh and vec_mergel is about 880 bytes shorter.
+
+#if defined(__GNUC__) && (__BIG_ENDIAN__)
+# define vec_merge_hi(a,b) vec_perm(a,b, HH_MASK)
+# define vec_merge_lo(a,b) vec_perm(a,b, LL_MASK)
+#else
+# define vec_merge_hi(a,b) vec_mergeh(a,b)
+# define vec_merge_lo(a,b) vec_mergel(a,b)
+#endif
+
+void BLAKE2_Compress64_POWER8(const byte* input, BLAKE2_State<word64, true>& state)
+{
+ // Permute masks. High is element 0 (most significant),
+ // low is element 1 (least significant).
+
+#if defined(__GNUC__) && (__BIG_ENDIAN__)
+ const uint8x16_p HH_MASK = { 0,1,2,3,4,5,6,7, 16,17,18,19,20,21,22,23 };
+ const uint8x16_p LL_MASK = { 8,9,10,11,12,13,14,15, 24,25,26,27,28,29,30,31 };
+#endif
+
+ const uint8x16_p HL_MASK = { 0,1,2,3,4,5,6,7, 24,25,26,27,28,29,30,31 };
+ const uint8x16_p LH_MASK = { 8,9,10,11,12,13,14,15, 16,17,18,19,20,21,22,23 };
+
+ #define BLAKE2B_LOAD_MSG_0_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m0, m1); \
+ b1 = vec_merge_hi(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_2(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m0, m1); \
+ b1 = vec_merge_lo(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_3(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m4, m5); \
+ b1 = vec_merge_hi(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_0_4(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m4, m5); \
+ b1 = vec_merge_lo(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m7, m2); \
+ b1 = vec_merge_lo(m4, m6); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_2(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m5, m4); \
+ b1 = vec_ext(m7, m3, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_3(b0, b1) \
+ do { \
+ b0 = vec_ext(m0, m0, 1); \
+ b1 = vec_merge_lo(m5, m2); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_1_4(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m6, m1); \
+ b1 = vec_merge_lo(m3, m1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_1(b0, b1) \
+ do { \
+ b0 = vec_ext(m5, m6, 1); \
+ b1 = vec_merge_lo(m2, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_2(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m4, m0); \
+ b1 = vec_perm(m1, m6, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_3(b0, b1) \
+ do { \
+ b0 = vec_perm(m5, m1, HL_MASK); \
+ b1 = vec_merge_lo(m3, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_2_4(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m7, m3); \
+ b1 = vec_ext(m0, m2, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_1(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m3, m1); \
+ b1 = vec_merge_lo(m6, m5); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_2(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m4, m0); \
+ b1 = vec_merge_hi(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_3(b0, b1) \
+ do { \
+ b0 = vec_perm(m1, m2, HL_MASK); \
+ b1 = vec_perm(m2, m7, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_3_4(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m3, m5); \
+ b1 = vec_merge_hi(m0, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_1(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m4, m2); \
+ b1 = vec_merge_hi(m1, m5); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_2(b0, b1) \
+ do { \
+ b0 = vec_perm(m0, m3, HL_MASK); \
+ b1 = vec_perm(m2, m7, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_3(b0, b1) \
+ do { \
+ b0 = vec_perm(m7, m5, HL_MASK); \
+ b1 = vec_perm(m3, m1, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_4_4(b0, b1) \
+ do { \
+ b0 = vec_ext(m0, m6, 1); \
+ b1 = vec_perm(m4, m6, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m1, m3); \
+ b1 = vec_merge_hi(m0, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_2(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m6, m5); \
+ b1 = vec_merge_lo(m5, m1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_3(b0, b1) \
+ do { \
+ b0 = vec_perm(m2, m3, HL_MASK); \
+ b1 = vec_merge_lo(m7, m0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_5_4(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m6, m2); \
+ b1 = vec_perm(m7, m4, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_1(b0, b1) \
+ do { \
+ b0 = vec_perm(m6, m0, HL_MASK); \
+ b1 = vec_merge_hi(m7, m2); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_2(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m2, m7); \
+ b1 = vec_ext(m6, m5, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_3(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m0, m3); \
+ b1 = vec_ext(m4, m4, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_6_4(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m3, m1); \
+ b1 = vec_perm(m1, m5, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_1(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m6, m3); \
+ b1 = vec_perm(m6, m1, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_2(b0, b1) \
+ do { \
+ b0 = vec_ext(m5, m7, 1); \
+ b1 = vec_merge_lo(m0, m4); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_3(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m2, m7); \
+ b1 = vec_merge_hi(m4, m1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_7_4(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m0, m2); \
+ b1 = vec_merge_hi(m3, m5); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m3, m7); \
+ b1 = vec_ext(m5, m0, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_2(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m7, m4); \
+ b1 = vec_ext(m1, m4, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_3(b0, b1) \
+ do { \
+ b0 = m6; \
+ b1 = vec_ext(m0, m5, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_8_4(b0, b1) \
+ do { \
+ b0 = vec_perm(m1, m3, HL_MASK); \
+ b1 = m2; \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m5, m4); \
+ b1 = vec_merge_lo(m3, m0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_2(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m1, m2); \
+ b1 = vec_perm(m3, m2, HL_MASK); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_3(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m7, m4); \
+ b1 = vec_merge_lo(m1, m6); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_9_4(b0, b1) \
+ do { \
+ b0 = vec_ext(m5, m7, 1); \
+ b1 = vec_merge_hi(m6, m0); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m0, m1); \
+ b1 = vec_merge_hi(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_2(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m0, m1); \
+ b1 = vec_merge_lo(m2, m3); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_3(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m4, m5); \
+ b1 = vec_merge_hi(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_10_4(b0, b1) \
+ do { \
+ b0 = vec_merge_lo(m4, m5); \
+ b1 = vec_merge_lo(m6, m7); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_1(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m7, m2); \
+ b1 = vec_merge_lo(m4, m6); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_2(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m5, m4); \
+ b1 = vec_ext(m7, m3, 1); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_3(b0, b1) \
+ do { \
+ b0 = vec_ext(m0, m0, 1); \
+ b1 = vec_merge_lo(m5, m2); \
+ } while(0)
+
+ #define BLAKE2B_LOAD_MSG_11_4(b0, b1) \
+ do { \
+ b0 = vec_merge_hi(m6, m1); \
+ b1 = vec_merge_lo(m3, m1); \
+ } while(0)
+
+ // Power8 has packed 64-bit rotate, but in terms of left rotate
+ const uint64x2_p ROR16_MASK = { 64-16, 64-16 };
+ const uint64x2_p ROR24_MASK = { 64-24, 64-24 };
+ const uint64x2_p ROR32_MASK = { 64-32, 64-32 };
+ const uint64x2_p ROR63_MASK = { 64-63, 64-63 };
+
+ #define vec_ror_32(x) vec_rl(x, ROR32_MASK)
+ #define vec_ror_24(x) vec_rl(x, ROR24_MASK)
+ #define vec_ror_16(x) vec_rl(x, ROR16_MASK)
+ #define vec_ror_63(x) vec_rl(x, ROR63_MASK)
+
+ #define BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1) \
+ do { \
+ row1l = vec_add(vec_add(row1l, b0), row2l); \
+ row1h = vec_add(vec_add(row1h, b1), row2h); \
+ row4l = vec_xor(row4l, row1l); row4h = vec_xor(row4h, row1h); \
+ row4l = vec_ror_32(row4l); row4h = vec_ror_32(row4h); \
+ row3l = vec_add(row3l, row4l); row3h = vec_add(row3h, row4h); \
+ row2l = vec_xor(row2l, row3l); row2h = vec_xor(row2h, row3h); \
+ row2l = vec_ror_24(row2l); row2h = vec_ror_24(row2h); \
+ } while(0)
+
+ #define BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1) \
+ do { \
+ row1l = vec_add(vec_add(row1l, b0), row2l); \
+ row1h = vec_add(vec_add(row1h, b1), row2h); \
+ row4l = vec_xor(row4l, row1l); row4h = vec_xor(row4h, row1h); \
+ row4l = vec_ror_16(row4l); row4h = vec_ror_16(row4h); \
+ row3l = vec_add(row3l, row4l); row3h = vec_add(row3h, row4h); \
+ row2l = vec_xor(row2l, row3l); row2h = vec_xor(row2h, row3h); \
+ row2l = vec_ror_63(row2l); row2h = vec_ror_63(row2h); \
+ } while(0)
+
+ #define BLAKE2B_DIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h) \
+ do { \
+ uint64x2_p t0 = vec_ext(row2l, row2h, 1); \
+ uint64x2_p t1 = vec_ext(row2h, row2l, 1); \
+ row2l = t0; row2h = t1; t0 = row3l; row3l = row3h; row3h = t0; \
+ t0 = vec_ext(row4h, row4l, 1); t1 = vec_ext(row4l, row4h, 1); \
+ row4l = t0; row4h = t1; \
+ } while(0)
+
+ #define BLAKE2B_UNDIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h) \
+ do { \
+ uint64x2_p t0 = vec_ext(row2h, row2l, 1); \
+ uint64x2_p t1 = vec_ext(row2l, row2h, 1); \
+ row2l = t0; row2h = t1; t0 = row3l; row3l = row3h; row3h = t0; \
+ t0 = vec_ext(row4l, row4h, 1); t1 = vec_ext(row4h, row4l, 1); \
+ row4l = t0; row4h = t1; \
+ } while(0)
+
+ #define BLAKE2B_ROUND(r) \
+ do { \
+ uint64x2_p b0, b1; \
+ BLAKE2B_LOAD_MSG_ ##r ##_1(b0, b1); \
+ BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_LOAD_MSG_ ##r ##_2(b0, b1); \
+ BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_DIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h); \
+ BLAKE2B_LOAD_MSG_ ##r ##_3(b0, b1); \
+ BLAKE2B_G1(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_LOAD_MSG_ ##r ##_4(b0, b1); \
+ BLAKE2B_G2(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h,b0,b1); \
+ BLAKE2B_UNDIAGONALIZE(row1l,row2l,row3l,row4l,row1h,row2h,row3h,row4h); \
+ } while(0)
+
+ const uint64x2_p m0 = VectorLoad64LE(input + 00);
+ const uint64x2_p m1 = VectorLoad64LE(input + 16);
+ const uint64x2_p m2 = VectorLoad64LE(input + 32);
+ const uint64x2_p m3 = VectorLoad64LE(input + 48);
+ const uint64x2_p m4 = VectorLoad64LE(input + 64);
+ const uint64x2_p m5 = VectorLoad64LE(input + 80);
+ const uint64x2_p m6 = VectorLoad64LE(input + 96);
+ const uint64x2_p m7 = VectorLoad64LE(input + 112);
+
+ uint64x2_p row1l, row1h, row2l, row2h;
+ uint64x2_p row3l, row3h, row4l, row4h;
+
+ const uint64x2_p h0 = row1l = VectorLoad64LE(&state.h[0]);
+ const uint64x2_p h1 = row1h = VectorLoad64LE(&state.h[2]);
+ const uint64x2_p h2 = row2l = VectorLoad64LE(&state.h[4]);
+ const uint64x2_p h3 = row2h = VectorLoad64LE(&state.h[6]);
+
+ row3l = VectorLoad64(&BLAKE2B_IV[0]);
+ row3h = VectorLoad64(&BLAKE2B_IV[2]);
+ row4l = vec_xor(VectorLoad64(&BLAKE2B_IV[4]), VectorLoad64(&state.t[0]));
+ row4h = vec_xor(VectorLoad64(&BLAKE2B_IV[6]), VectorLoad64(&state.f[0]));
+
+ BLAKE2B_ROUND(0);
+ BLAKE2B_ROUND(1);
+ BLAKE2B_ROUND(2);
+ BLAKE2B_ROUND(3);
+ BLAKE2B_ROUND(4);
+ BLAKE2B_ROUND(5);
+ BLAKE2B_ROUND(6);
+ BLAKE2B_ROUND(7);
+ BLAKE2B_ROUND(8);
+ BLAKE2B_ROUND(9);
+ BLAKE2B_ROUND(10);
+ BLAKE2B_ROUND(11);
+
+ VectorStore64LE(&state.h[0], vec_xor(h0, vec_xor(row1l, row3l)));
+ VectorStore64LE(&state.h[2], vec_xor(h1, vec_xor(row1h, row3h)));
+ VectorStore64LE(&state.h[4], vec_xor(h2, vec_xor(row2l, row4l)));
+ VectorStore64LE(&state.h[6], vec_xor(h3, vec_xor(row2h, row4h)));
+}
+#endif // CRYPTOPP_POWER8_AVAILABLE
+
+NAMESPACE_END
|