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
|
/* This file generated by gen-psqr.c - DO NOT EDIT. */
#if GMP_LIMB_BITS != 32 || GMP_NAIL_BITS != 0
Error, error, this data is for 32 bit limb and 0 bit nail
#endif
/* Non-zero bit indicates a quadratic residue mod 0x100.
This test identifies 82.81% as non-squares (212/256). */
static const mp_limb_t
sq_res_0x100[8] = {
CNST_LIMB(0x2030213),
CNST_LIMB(0x2020212),
CNST_LIMB(0x2020213),
CNST_LIMB(0x2020212),
CNST_LIMB(0x2030212),
CNST_LIMB(0x2020212),
CNST_LIMB(0x2020212),
CNST_LIMB(0x2020212),
};
/* 2^24-1 = 3^2 * 5 * 7 * 13 * 17 ... */
#define PERFSQR_MOD_BITS 25
/* This test identifies 95.66% as non-squares. */
#define PERFSQR_MOD_TEST(up, usize) \
do { \
mp_limb_t r; \
PERFSQR_MOD_34 (r, up, usize); \
\
/* 73.33% */ \
PERFSQR_MOD_2 (r, CNST_LIMB(45), CNST_LIMB(0xfa4fa5), \
CNST_LIMB(0x920), CNST_LIMB(0x1a442481)); \
\
/* 47.06% */ \
PERFSQR_MOD_1 (r, CNST_LIMB(17), CNST_LIMB(0xf0f0f1), \
CNST_LIMB(0x1a317)); \
\
/* 46.15% */ \
PERFSQR_MOD_1 (r, CNST_LIMB(13), CNST_LIMB(0xec4ec5), \
CNST_LIMB(0x9e5)); \
\
/* 42.86% */ \
PERFSQR_MOD_1 (r, CNST_LIMB( 7), CNST_LIMB(0xdb6db7), \
CNST_LIMB(0x69)); \
} while (0)
/* Grand total sq_res_0x100 and PERFSQR_MOD_TEST, 99.25% non-squares. */
/* helper for tests/mpz/t-perfsqr.c */
#define PERFSQR_DIVISORS { 256, 45, 17, 13, 7, }
|
