2 * Generate random test vectors for 64-bit operations.
4 * As an independent reference, we use the Catacomb multiprecision integer
10 #include <catacomb/mp.h>
11 #include <catacomb/mprand.h>
12 #include <catacomb/fibrand.h>
19 static mp *mp_rol(mp *d, mp *x, unsigned s)
21 mp *l = mp_lsl(MP_NEW, x, s);
22 mp *r = mp_lsr(MP_NEW, x, 64 - s);
29 static mp *mp_ror(mp *d, mp *x, unsigned s)
31 mp *l = mp_lsl(MP_NEW, x, 64 - s);
32 mp *r = mp_lsr(MP_NEW, x, s);
39 static void putmp(mp *x)
44 mp_storeb(x, buf, sizeof(buf));
45 for (i = 0; i < sizeof(buf); i++)
46 printf("%02x", buf[i]);
49 #define GEN_SHIFT(op) \
51 static void gen_##op(void) \
54 fputs("\n" #op "64 {\n", stdout); \
55 for (i = 0; i < NVEC; i++) { \
56 mp *x = mprand_range(MP_NEW, m64, r, 0); \
57 unsigned s = r->ops->range(r, 70), ss = s & 63; \
58 mp *y = mp_##op(MP_NEW, x, ss); \
59 mp_div(0, &y, y, m64); \
62 putmp(x); printf(" %2u", s); putmp(y); \
63 fputs(";\n", stdout); \
64 mp_drop(x); mp_drop(y); \
66 fputs("}\n", stdout); \
69 #define GEN_ARITH(op) \
71 static void gen_##op(void) \
74 fputs("\n" #op "64 {\n", stdout); \
75 for (i = 0; i < NVEC; i++) { \
76 mp *x = mprand_range(MP_NEW, m64, r, 0); \
77 mp *y = mprand_range(MP_NEW, m64, r, 0); \
78 mp *z = mp_##op(MP_NEW, x, y); \
79 mp_div(0, &z, z, m64); \
82 putmp(x); putmp(y); putmp(z); \
83 fputs(";\n", stdout); \
84 mp_drop(x); mp_drop(y); mp_drop(z); \
86 fputs("}\n", stdout); \
98 m64 = mp_lsl(MP_NEW, MP_ONE, 64);
99 r = fibrand_create(0);
100 fputs("# Test vectors for 64-bit operations [generated]\n", stdout);