Annotation of OpenXM_contrib/gmp/tests/mpz/t-pow.c, Revision 1.1.1.1
1.1 ohara 1: /* Test mpz_pow_ui and mpz_ui_pow_ui.
2:
3: Copyright 1997, 1999, 2000, 2001 Free Software Foundation, Inc.
4:
5: This file is part of the GNU MP Library.
6:
7: The GNU MP Library is free software; you can redistribute it and/or modify
8: it under the terms of the GNU Lesser General Public License as published by
9: the Free Software Foundation; either version 2.1 of the License, or (at your
10: option) any later version.
11:
12: The GNU MP Library is distributed in the hope that it will be useful, but
13: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
14: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
15: License for more details.
16:
17: You should have received a copy of the GNU Lesser General Public License
18: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
19: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
20: MA 02111-1307, USA.
21: */
22:
23: #include <stdio.h>
24: #include <stdlib.h>
25:
26: #include "gmp.h"
27: #include "gmp-impl.h"
28: #include "tests.h"
29:
30:
31: void
32: check_one (mpz_srcptr want, mpz_srcptr base, unsigned long exp)
33: {
34: mpz_t got;
35:
36: mpz_init (got);
37:
38: mpz_pow_ui (got, base, exp);
39: if (mpz_cmp (got, want))
40: {
41: printf ("mpz_pow_ui wrong\n");
42: mpz_trace (" base", base);
43: printf (" exp = %lu (0x%lX)\n", exp, exp);
44: mpz_trace (" got ", got);
45: mpz_trace (" want", want);
46: abort ();
47: }
48:
49: mpz_set (got, base);
50: mpz_pow_ui (got, got, exp);
51: if (mpz_cmp (got, want))
52: {
53: printf ("mpz_pow_ui wrong\n");
54: mpz_trace (" base", base);
55: printf (" exp = %lu (0x%lX)\n", exp, exp);
56: mpz_trace (" got ", got);
57: mpz_trace (" want", want);
58: abort ();
59: }
60:
61: if (mpz_fits_ulong_p (base))
62: {
63: unsigned long base_u = mpz_get_ui (base);
64: mpz_ui_pow_ui (got, base_u, exp);
65: if (mpz_cmp (got, want))
66: {
67: printf ("mpz_ui_pow_ui wrong\n");
68: printf (" base=%lu (0x%lX)\n", base_u, base_u);
69: printf (" exp = %lu (0x%lX)\n", exp, exp);
70: mpz_trace (" got ", got);
71: mpz_trace (" want", want);
72: abort ();
73: }
74: }
75:
76: mpz_clear (got);
77: }
78:
79: void
80: check_base (mpz_srcptr base)
81: {
82: unsigned long exp;
83: mpz_t want;
84:
85: mpz_init (want);
86: mpz_set_ui (want, 1L);
87:
88: for (exp = 0; exp < 20; exp++)
89: {
90: check_one (want, base, exp);
91: mpz_mul (want, want, base);
92: }
93:
94: mpz_clear (want);
95: }
96:
97: void
98: check_various (void)
99: {
100: static const struct {
101: const char *base;
102: } data[] = {
103: { "0" },
104: { "1" },
105: { "2" },
106: { "3" },
107: { "4" },
108: { "5" },
109: { "6" },
110: { "10" },
111: { "15" },
112: { "16" },
113:
114: { "0x1F" },
115: { "0xFF" },
116: { "0x1001" },
117: { "0xFFFF" },
118: { "0x10000001" },
119: { "0x1000000000000001" },
120:
121: /* actual size closest to estimate */
122: { "0xFFFFFFFF" },
123: { "0xFFFFFFFFFFFFFFFF" },
124:
125: /* same after rshift */
126: { "0xFFFFFFFF0" },
127: { "0xFFFFFFFF00" },
128: { "0xFFFFFFFFFFFFFFFF0" },
129: { "0xFFFFFFFFFFFFFFFF00" },
130:
131: /* change from 2 limbs to 1 after rshift */
132: { "0x180000000" },
133: { "0x18000000000000000" },
134:
135: /* change from 3 limbs to 2 after rshift */
136: { "0x18000000100000000" },
137: { "0x180000000000000010000000000000000" },
138:
139: /* handling of absolute value */
140: { "-0x80000000" },
141: { "-0x8000000000000000" },
142:
143: /* low zero limb, and size>2, checking argument overlap detection */
144: { "0x3000000000000000300000000000000030000000000000000" },
145: };
146:
147: mpz_t base;
148: int i;
149:
150: mpz_init (base);
151:
152: for (i = 0; i < numberof (data); i++)
153: {
154: mpz_set_str_or_abort (base, data[i].base, 0);
155: check_base (base);
156: }
157:
158: mpz_clear (base);
159: }
160:
161: void
162: check_random (int reps)
163: {
164: mpz_t base, want;
165: mp_size_t base_size;
166: int i;
167: unsigned long size_range, exp;
168: gmp_randstate_ptr rands = RANDS;
169:
170: mpz_init (base);
171: mpz_init (want);
172:
173: for (i = 0; i < reps; i++)
174: {
175: /* exponentially random 0 to 2^13 bits for base */
176: mpz_urandomb (want, rands, 32);
177: size_range = mpz_get_ui (want) % 12 + 2;
178: mpz_urandomb (want, rands, size_range);
179: base_size = mpz_get_ui (want);
180: mpz_rrandomb (base, rands, base_size);
181:
182: /* randomly signed base */
183: mpz_urandomb (want, rands, 2);
184: if ((mpz_get_ui (want) & 1) != 0)
185: mpz_neg (base, base);
186:
187: /* random 5 bits for exponent */
188: mpz_urandomb (want, rands, 5L);
189: exp = mpz_get_ui (want);
190:
191: refmpz_pow_ui (want, base, exp);
192: check_one (want, base, exp);
193: }
194:
195: mpz_clear (base);
196: mpz_clear (want);
197: }
198:
199: int
200: main (int argc, char **argv)
201: {
202: int reps = 500;
203:
204: /* dummy call to drag in refmpn.o for testing mpz/n_pow_ui.c with
205: refmpn_mul_2 */
206: refmpn_zero_p (NULL, (mp_size_t) 0);
207:
208: tests_start ();
209: mp_trace_base = -16;
210:
211: if (argc == 2)
212: reps = atoi (argv[1]);
213:
214: check_various ();
215: check_random (reps);
216:
217: tests_end ();
218: exit (0);
219: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to t-pow.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