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/>.
29 #define ROTL(x,b) (u64)( ((x) << (b)) | ( (x) >> (64 - (b))) )
31 #define U32TO8_LE(p, v) \
32 (p)[0] = (u8)((v) ); (p)[1] = (u8)((v) >> 8); \
33 (p)[2] = (u8)((v) >> 16); (p)[3] = (u8)((v) >> 24);
35 #define U64TO8_LE(p, v) \
36 U32TO8_LE((p), (u32)((v) )); \
37 U32TO8_LE((p) + 4, (u32)((v) >> 32));
39 #define U8TO64_LE(p) \
41 ((u64)((p)[1]) << 8) | \
42 ((u64)((p)[2]) << 16) | \
43 ((u64)((p)[3]) << 24) | \
44 ((u64)((p)[4]) << 32) | \
45 ((u64)((p)[5]) << 40) | \
46 ((u64)((p)[6]) << 48) | \
47 ((u64)((p)[7]) << 56))
51 v0 += v1; v1=ROTL(v1,13); v1 ^= v0; v0=ROTL(v0,32); \
52 v2 += v3; v3=ROTL(v3,16); v3 ^= v2; \
53 v0 += v3; v3=ROTL(v3,21); v3 ^= v0; \
54 v2 += v1; v1=ROTL(v1,17); v1 ^= v2; v2=ROTL(v2,32); \
58 void g_siphash24(guint8 out[8], const void *_in, gsize inlen, const guint8 k[16])
60 /* "somepseudorandomlygeneratedbytes" */
61 u64 v0 = 0x736f6d6570736575ULL;
62 u64 v1 = 0x646f72616e646f6dULL;
63 u64 v2 = 0x6c7967656e657261ULL;
64 u64 v3 = 0x7465646279746573ULL;
66 u64 k0 = U8TO64_LE( k );
67 u64 k1 = U8TO64_LE( k + 8 );
70 const u8 *end = in + inlen - ( inlen % sizeof( u64 ) );
71 const int left = inlen & 7;
72 b = ( ( u64 )inlen ) << 56;
78 for ( ; in != end; in += 8 )
82 printf( "(%3d) v0 %08x %08x\n", ( int )inlen, ( u32 )( v0 >> 32 ), ( u32 )v0 );
83 printf( "(%3d) v1 %08x %08x\n", ( int )inlen, ( u32 )( v1 >> 32 ), ( u32 )v1 );
84 printf( "(%3d) v2 %08x %08x\n", ( int )inlen, ( u32 )( v2 >> 32 ), ( u32 )v2 );
85 printf( "(%3d) v3 %08x %08x\n", ( int )inlen, ( u32 )( v3 >> 32 ), ( u32 )v3 );
86 printf( "(%3d) compress %08x %08x\n", ( int )inlen, ( u32 )( m >> 32 ), ( u32 )m );
96 case 7: b |= ( ( u64 )in[ 6] ) << 48;
98 case 6: b |= ( ( u64 )in[ 5] ) << 40;
100 case 5: b |= ( ( u64 )in[ 4] ) << 32;
102 case 4: b |= ( ( u64 )in[ 3] ) << 24;
104 case 3: b |= ( ( u64 )in[ 2] ) << 16;
106 case 2: b |= ( ( u64 )in[ 1] ) << 8;
108 case 1: b |= ( ( u64 )in[ 0] ); break;
114 printf( "(%3d) v0 %08x %08x\n", ( int )inlen, ( u32 )( v0 >> 32 ), ( u32 )v0 );
115 printf( "(%3d) v1 %08x %08x\n", ( int )inlen, ( u32 )( v1 >> 32 ), ( u32 )v1 );
116 printf( "(%3d) v2 %08x %08x\n", ( int )inlen, ( u32 )( v2 >> 32 ), ( u32 )v2 );
117 printf( "(%3d) v3 %08x %08x\n", ( int )inlen, ( u32 )( v3 >> 32 ), ( u32 )v3 );
118 printf( "(%3d) padding %08x %08x\n", ( int )inlen, ( u32 )( b >> 32 ), ( u32 )b );
125 printf( "(%3d) v0 %08x %08x\n", ( int )inlen, ( u32 )( v0 >> 32 ), ( u32 )v0 );
126 printf( "(%3d) v1 %08x %08x\n", ( int )inlen, ( u32 )( v1 >> 32 ), ( u32 )v1 );
127 printf( "(%3d) v2 %08x %08x\n", ( int )inlen, ( u32 )( v2 >> 32 ), ( u32 )v2 );
128 printf( "(%3d) v3 %08x %08x\n", ( int )inlen, ( u32 )( v3 >> 32 ), ( u32 )v3 );
135 b = v0 ^ v1 ^ v2 ^ v3;