dump/visualising-siphash/generate.c at master · majek/dump · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>
#include <math.h>


#ifdef __x86_64__
#  define RDRAND_LONG ".byte 0x48,0x0f,0xc7,0xf0"
#else
#  define RDRAND_INT  ".byte 0x0f,0xc7,0xf0"
#  define RDRAND_LONG   RDRAND_INT
#endif

#define RDRAND_RETRY_LOOPS  10

static inline long rdrand_long(unsigned long *v) {
    int ok;
    asm volatile("1: " RDRAND_LONG "\n\t"
             "jc 2f\n\t"
             "decl %0\n\t"
             "jnz 1b\n\t"
             "2:"
             : "=r" (ok), "=a" (*v)
             : "0" (RDRAND_RETRY_LOOPS));
    return ok;
}

static void fill_random(void *buf, unsigned buf_sz) {
	unsigned long *b = buf;
	unsigned b_sz = buf_sz / sizeof(long);

	unsigned i;
	for (i = 0; i < b_sz; i++) {
		rdrand_long(&b[i]);
	}
}


#define ROTATE(x,b) (uint64_t)( ((x) << (b)) | ( (x) >> (64 - (b))) )

#define HALF_ROUND(a,b,c,d,s,t)			\
	a += b; c += d;				\
	b = ROTATE(b, s) ^ a;	      		\
	d = ROTATE(d, t) ^ c;			\
	a = ROTATE(a, 32);

#define DOUBLE_ROUND(v0,v1,v2,v3)		\
	HALF_ROUND(v0,v1,v2,v3,13,16);		\
	HALF_ROUND(v2,v1,v0,v3,17,21);		\
	HALF_ROUND(v0,v1,v2,v3,13,16);		\
	HALF_ROUND(v2,v1,v0,v3,17,21);

#define SINGLE_ROUND(v0,v1,v2,v3)		\
	HALF_ROUND(v0,v1,v2,v3,13,16);		\
	HALF_ROUND(v2,v1,v0,v3,17,21);


/* static void sipround(uint64_t v[4]) { */
/* 	DOUBLE_ROUND(v[0],v[1],v[2],v[3]); */
/* 	v[0] ^= b; v[2] ^= 0xff; */
/* 	DOUBLE_ROUND(v[0],v[1],v[2],v[3]); */
/* 	DOUBLE_ROUND(v[0],v[1],v[2],v[3]); */

/* //	SINGLE_ROUND(v[0],v[1],v[2],v[3]); */
/* //	DOUBLE_ROUND(v[0],v[1],v[2],v[3]); */
/* } */

void siphash24(uint64_t *v) {
	uint64_t k0 = v[0];
	uint64_t k1 = v[1];
	uint64_t s = v[2];


	uint64_t len = 7;
	//s =  0;
	// len = 0
	uint64_t b = len << 56;
	b |= s & 0x00ffffffffffffffULL;


	uint64_t v0 = k0 ^ 0x736f6d6570736575ULL;
	uint64_t v1 = k1 ^ 0x646f72616e646f6dULL;
	uint64_t v2 = k0 ^ 0x6c7967656e657261ULL;
	uint64_t v3 = k1 ^ 0x7465646279746573ULL;

	v3 ^= b;
	DOUBLE_ROUND(v0,v1,v2,v3);
	v0 ^= b; v2 ^= 0xff;
	DOUBLE_ROUND(v0,v1,v2,v3);
//	DOUBLE_ROUND(v0,v1,v2,v3);
	SINGLE_ROUND(v0,v1,v2,v3);

	v[0] = v0;
	v[1] = v1;
	v[2] = v2;
	v[3] = v3;
}


int main() {
	const int rounds = 1000000000;

	int bit;
	for (bit = 0; bit < 256; bit++) {
		uint64_t c[256] = {0};

		if (bit >= 192-8) {
			unsigned i;
			for (i=0; i < 256; i++)
				c[i] = rounds / 2;
		} else {
			uint64_t bitmask = 1ULL << (bit & 63);

			int r;
			for (r = 0; r < rounds; r++) {
				uint64_t a[4], b[4], diff[4];
				fill_random(a, sizeof(a));
				memcpy(b, a, sizeof(b));

				a[bit / 64] |= bitmask;
				b[bit / 64] &= ~bitmask;

				siphash24(a);
				siphash24(b);

				unsigned i, j;
				for (i=0; i < 4; i++)
					diff[i] = a[i]^b[i];
				for (i=0; i < 64; i++)
					for (j=0; j < 4; j++)
						c[j*64+i] += !!(diff[j] & (1ULL << i));
			}
		}

		unsigned i;
		for(i=0; i < 256; i++) {
			printf("%i\t%i\t%.20f\n",
			       bit, i,
			       c[i] / (double)rounds);
		}
		printf("\n");
	}
}