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
|
/* Test mpn_and, mpn_ior, mpn_xor, mpn_andn, mpn_iorn, mpn_xnor, mpn_nand, and
mpn_nior.
Copyright 2011-2013 Free Software Foundation, Inc.
This file is part of the GNU MP Library test suite.
The GNU MP Library test suite is free software; you can redistribute it
and/or modify it under the terms of the GNU General Public License as
published by the Free Software Foundation; either version 3 of the License,
or (at your option) any later version.
The GNU MP Library test suite is distributed in the hope that it will be
useful, but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
Public License for more details.
You should have received a copy of the GNU General Public License along with
the GNU MP Library test suite. If not, see https://www.gnu.org/licenses/. */
#include <stdlib.h>
#include <stdio.h>
/* Fake native prevalence of the tested operations, so that we actually test
the compiled functions, i.e., the ones which users will reach. The inlined
variants will be tested through tests/mpz/logic.c. */
#define HAVE_NATIVE_mpn_com 1
#define HAVE_NATIVE_mpn_and_n 1
#define HAVE_NATIVE_mpn_andn_n 1
#define HAVE_NATIVE_mpn_nand_n 1
#define HAVE_NATIVE_mpn_ior_n 1
#define HAVE_NATIVE_mpn_iorn_n 1
#define HAVE_NATIVE_mpn_nior_n 1
#define HAVE_NATIVE_mpn_xor_n 1
#define HAVE_NATIVE_mpn_xnor_n 1
#include "gmp-impl.h"
#include "tests.h"
void
check_one (mp_srcptr refp, mp_srcptr rp, mp_srcptr ap, mp_srcptr bp, mp_size_t n, const char *funcname)
{
if (mpn_cmp (refp, rp, n))
{
printf ("ERROR in mpn_%s\n", funcname);
printf ("a: "); mpn_dump (ap, n);
printf ("b: "); mpn_dump (bp, n);
printf ("r: "); mpn_dump (rp, n);
printf ("ref: "); mpn_dump (refp, n);
abort();
}
}
int
main (int argc, char **argv)
{
mpz_t a, b;
mp_ptr ap, bp, rp, refp;
mp_size_t max_n, n, i;
gmp_randstate_ptr rands;
long test, reps = 1000;
TMP_DECL;
TMP_MARK;
tests_start ();
TESTS_REPS (reps, argv, argc);
mpz_inits (a, b, NULL);
rands = RANDS; /* FIXME: not used */
max_n = 100;
rp = TMP_ALLOC_LIMBS (1 + max_n * 8 / GMP_LIMB_BITS);
refp = TMP_ALLOC_LIMBS (1 + max_n * 8 / GMP_LIMB_BITS);
for (test = 0; test < reps; test++)
{
for (i = 1; i <= max_n; i++)
{
mpz_rrandomb (a, rands, i * 8);
mpz_rrandomb (b, rands, i * 8);
mpz_setbit (a, i * 8 - 1);
mpz_setbit (b, i * 8 - 1);
ap = PTR(a);
bp = PTR(b);
n = SIZ(a);
refmpn_and_n (refp, ap, bp, n);
mpn_and_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "and_n");
refmpn_ior_n (refp, ap, bp, n);
mpn_ior_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "ior_n");
refmpn_xor_n (refp, ap, bp, n);
mpn_xor_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "xor_n");
refmpn_andn_n (refp, ap, bp, n);
mpn_andn_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "andn_n");
refmpn_iorn_n (refp, ap, bp, n);
mpn_iorn_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "iorn_n");
refmpn_nand_n (refp, ap, bp, n);
mpn_nand_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "nand_n");
refmpn_nior_n (refp, ap, bp, n);
mpn_nior_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "nior_n");
refmpn_xnor_n (refp, ap, bp, n);
mpn_xnor_n (rp, ap, bp, n);
check_one (refp, rp, ap, bp, n, "xnor_n");
refmpn_com (refp, ap, n);
mpn_com (rp, ap, n);
check_one (refp, rp, ap, bp, n, "com");
}
}
TMP_FREE;
mpz_clears (a, b, NULL);
tests_end ();
return 0;
}
|
