Annotation of OpenXM_contrib/gmp/tests/mpz/t-cong.c, Revision 1.1.1.1
1.1 ohara 1: /* test mpz_congruent_p and mpz_congruent_ui_p */
2:
3: /*
4: Copyright 2001 Free Software Foundation, Inc.
5:
6: This file is part of the GNU MP Library.
7:
8: The GNU MP Library is free software; you can redistribute it and/or modify
9: it under the terms of the GNU Lesser General Public License as published by
10: the Free Software Foundation; either version 2.1 of the License, or (at your
11: option) any later version.
12:
13: The GNU MP Library is distributed in the hope that it will be useful, but
14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
16: License for more details.
17:
18: You should have received a copy of the GNU Lesser General Public License
19: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
21: MA 02111-1307, USA.
22: */
23:
24: #include <stdio.h>
25: #include <stdlib.h>
26: #include "gmp.h"
27: #include "gmp-impl.h"
28: #include "tests.h"
29:
30:
31: void
32: check_one (mpz_srcptr a, mpz_srcptr c, mpz_srcptr d, int want)
33: {
34: int got;
35: int swap;
36:
37: for (swap = 0; swap <= 1; swap++)
38: {
39: got = (mpz_congruent_p (a, c, d) != 0);
40: if (want != got)
41: {
42: printf ("mpz_congruent_p wrong\n");
43: printf (" expected %d got %d\n", want, got);
44: mpz_trace (" a", a);
45: mpz_trace (" c", c);
46: mpz_trace (" d", d);
47: mp_trace_base = -16;
48: mpz_trace (" a", a);
49: mpz_trace (" c", c);
50: mpz_trace (" d", d);
51: abort ();
52: }
53:
54: if (mpz_fits_ulong_p (c) && mpz_fits_ulong_p (d))
55: {
56: unsigned long uc = mpz_get_ui (c);
57: unsigned long ud = mpz_get_ui (d);
58: got = (mpz_congruent_ui_p (a, uc, ud) != 0);
59: if (want != got)
60: {
61: printf ("mpz_congruent_ui_p wrong\n");
62: printf (" expected %d got %d\n", want, got);
63: mpz_trace (" a", a);
64: printf (" c=%lu\n", uc);
65: printf (" d=%lu\n", ud);
66: mp_trace_base = -16;
67: mpz_trace (" a", a);
68: printf (" c=0x%lX\n", uc);
69: printf (" d=0x%lX\n", ud);
70: abort ();
71: }
72: }
73:
74: MPZ_SRCPTR_SWAP (a, c);
75: }
76: }
77:
78:
79: void
80: check_data (void)
81: {
82: static const struct {
83: const char *a;
84: const char *c;
85: const char *d;
86: int want;
87:
88: } data[] = {
89:
90: /* anything congruent mod 1 */
91: { "0", "0", "1", 1 },
92: { "1", "0", "1", 1 },
93: { "0", "1", "1", 1 },
94: { "123", "456", "1", 1 },
95: { "0x123456789123456789", "0x987654321987654321", "1", 1 },
96:
97: /* csize==1, dsize==2 changing to 1 after stripping 2s */
98: { "0x3333333333333333", "0x33333333",
99: "0x180000000", 1 },
100: { "0x33333333333333333333333333333333", "0x3333333333333333",
101: "0x18000000000000000", 1 },
102:
103: /* another dsize==2 becoming 1, with opposite signs this time */
104: { "0x444444441",
105: "-0x22222221F",
106: "0x333333330", 1 },
107: { "0x44444444444444441",
108: "-0x2222222222222221F",
109: "0x33333333333333330", 1 },
110: };
111:
112: mpz_t a, c, d;
113: int i;
114:
115: mpz_init (a);
116: mpz_init (c);
117: mpz_init (d);
118:
119: for (i = 0; i < numberof (data); i++)
120: {
121: mpz_set_str_or_abort (a, data[i].a, 0);
122: mpz_set_str_or_abort (c, data[i].c, 0);
123: mpz_set_str_or_abort (d, data[i].d, 0);
124: check_one (a, c, d, data[i].want);
125: }
126:
127: mpz_clear (a);
128: mpz_clear (c);
129: mpz_clear (d);
130: }
131:
132:
133: void
134: check_random (int argc, char *argv[])
135: {
136: gmp_randstate_ptr rands = RANDS;
137: mpz_t a, c, d, ra, rc;
138: int i;
139: int want;
140: int reps = 2000;
141:
142: if (argc >= 2)
143: reps = atoi (argv[1]);
144:
145: mpz_init (a);
146: mpz_init (c);
147: mpz_init (d);
148: mpz_init (ra);
149: mpz_init (rc);
150:
151: for (i = 0; i < reps; i++)
152: {
153: mpz_errandomb (a, rands, 8*BITS_PER_MP_LIMB);
154: MPZ_CHECK_FORMAT (a);
155: mpz_errandomb (c, rands, 8*BITS_PER_MP_LIMB);
156: MPZ_CHECK_FORMAT (c);
157: mpz_errandomb_nonzero (d, rands, 8*BITS_PER_MP_LIMB);
158:
159: mpz_negrandom (a, rands);
160: MPZ_CHECK_FORMAT (a);
161: mpz_negrandom (c, rands);
162: MPZ_CHECK_FORMAT (c);
163: mpz_negrandom (d, rands);
164:
165: mpz_fdiv_r (ra, a, d);
166: mpz_fdiv_r (rc, c, d);
167:
168: want = (mpz_cmp (ra, rc) == 0);
169: check_one (a, c, d, want);
170:
171: mpz_sub (ra, ra, rc);
172: mpz_sub (a, a, ra);
173: MPZ_CHECK_FORMAT (a);
174: check_one (a, c, d, 1);
175:
176: if (! mpz_pow2abs_p (d))
177: {
178: mpz_flipbit (a, urandom() % (8*BITS_PER_MP_LIMB));
179: check_one (a, c, d, 0);
180: }
181: }
182:
183: mpz_clear (a);
184: mpz_clear (c);
185: mpz_clear (d);
186: mpz_clear (ra);
187: mpz_clear (rc);
188: }
189:
190:
191: int
192: main (int argc, char *argv[])
193: {
194: tests_start ();
195:
196: check_data ();
197: check_random (argc, argv);
198:
199: tests_end ();
200: exit (0);
201: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to t-cong.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / tests / mpz