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
);