2 SipHash reference C implementation
5 Jean-Philippe Aumasson <jeanphilippe.aumasson@gmail.com>
6 Daniel J. Bernstein <djb@cr.yp.to>
8 To the extent possible under law, the author(s) have dedicated all copyright
9 and related and neighboring rights to this software to the public domain
10 worldwide. This software is distributed without any warranty.
12 You should have received a copy of the CC0 Public Domain Dedication along with
13 this software. If not, see <http://creativecommons.org/publicdomain/zero/1.0/>.
15 (Minimal changes made by Lennart Poettering, to make clean for inclusion in systemd)
21 #include "siphash24.h"
27 #define ROTL(x,b) (u64)( ((x) << (b)) | ( (x) >> (64 - (b))) )
29 #define U32TO8_LE(p, v) \
30 (p)[0] = (u8)((v) ); (p)[1] = (u8)((v) >> 8); \
31 (p)[2] = (u8)((v) >> 16); (p)[3] = (u8)((v) >> 24);
33 #define U64TO8_LE(p, v) \
34 U32TO8_LE((p), (u32)((v) )); \
35 U32TO8_LE((p) + 4, (u32)((v) >> 32));
37 #define U8TO64_LE(p) \
39 ((u64)((p)[1]) << 8) | \
40 ((u64)((p)[2]) << 16) | \
41 ((u64)((p)[3]) << 24) | \
42 ((u64)((p)[4]) << 32) | \
43 ((u64)((p)[5]) << 40) | \
44 ((u64)((p)[6]) << 48) | \
45 ((u64)((p)[7]) << 56))
49 v0 += v1; v1=ROTL(v1,13); v1 ^= v0; v0=ROTL(v0,32); \
50 v2 += v3; v3=ROTL(v3,16); v3 ^= v2; \
51 v0 += v3; v3=ROTL(v3,21); v3 ^= v0; \
52 v2 += v1; v1=ROTL(v1,17); v1 ^= v2; v2=ROTL(v2,32); \
56 void siphash24(uint8_t out[8], const void *_in, size_t inlen, const uint8_t k[16])
58 /* "somepseudorandomlygeneratedbytes" */
59 u64 v0 = 0x736f6d6570736575ULL;
60 u64 v1 = 0x646f72616e646f6dULL;
61 u64 v2 = 0x6c7967656e657261ULL;
62 u64 v3 = 0x7465646279746573ULL;
64 u64 k0 = U8TO64_LE( k );
65 u64 k1 = U8TO64_LE( k + 8 );
68 const u8 *end = in + inlen - ( inlen % sizeof( u64 ) );
69 const int left = inlen & 7;
70 b = ( ( u64 )inlen ) << 56;
76 for ( ; in != end; in += 8 )
80 printf( "(%3d) v0 %08x %08x\n", ( int )inlen, ( u32 )( v0 >> 32 ), ( u32 )v0 );
81 printf( "(%3d) v1 %08x %08x\n", ( int )inlen, ( u32 )( v1 >> 32 ), ( u32 )v1 );
82 printf( "(%3d) v2 %08x %08x\n", ( int )inlen, ( u32 )( v2 >> 32 ), ( u32 )v2 );
83 printf( "(%3d) v3 %08x %08x\n", ( int )inlen, ( u32 )( v3 >> 32 ), ( u32 )v3 );
84 printf( "(%3d) compress %08x %08x\n", ( int )inlen, ( u32 )( m >> 32 ), ( u32 )m );
94 case 7: b |= ( ( u64 )in[ 6] ) << 48;
96 case 6: b |= ( ( u64 )in[ 5] ) << 40;
98 case 5: b |= ( ( u64 )in[ 4] ) << 32;
100 case 4: b |= ( ( u64 )in[ 3] ) << 24;
102 case 3: b |= ( ( u64 )in[ 2] ) << 16;
104 case 2: b |= ( ( u64 )in[ 1] ) << 8;
106 case 1: b |= ( ( u64 )in[ 0] ); break;
112 printf( "(%3d) v0 %08x %08x\n", ( int )inlen, ( u32 )( v0 >> 32 ), ( u32 )v0 );
113 printf( "(%3d) v1 %08x %08x\n", ( int )inlen, ( u32 )( v1 >> 32 ), ( u32 )v1 );
114 printf( "(%3d) v2 %08x %08x\n", ( int )inlen, ( u32 )( v2 >> 32 ), ( u32 )v2 );
115 printf( "(%3d) v3 %08x %08x\n", ( int )inlen, ( u32 )( v3 >> 32 ), ( u32 )v3 );
116 printf( "(%3d) padding %08x %08x\n", ( int )inlen, ( u32 )( b >> 32 ), ( u32 )b );
123 printf( "(%3d) v0 %08x %08x\n", ( int )inlen, ( u32 )( v0 >> 32 ), ( u32 )v0 );
124 printf( "(%3d) v1 %08x %08x\n", ( int )inlen, ( u32 )( v1 >> 32 ), ( u32 )v1 );
125 printf( "(%3d) v2 %08x %08x\n", ( int )inlen, ( u32 )( v2 >> 32 ), ( u32 )v2 );
126 printf( "(%3d) v3 %08x %08x\n", ( int )inlen, ( u32 )( v3 >> 32 ), ( u32 )v3 );
133 b = v0 ^ v1 ^ v2 ^ v3;