Annotation of OpenXM_contrib/gmp/tests/mpn/t-aors_1.c, Revision 1.1.1.1
1.1 ohara 1: /* Test mpn_add_1 and mpn_sub_1.
2:
3: Copyright 2001, 2002 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: #define M GMP_NUMB_MAX
32: #define ASIZE 10
33: #define MAGIC 0x1234
34:
35: #define SETUP() \
36: do { \
37: refmpn_random (got, data[i].size); \
38: got[data[i].size] = MAGIC; \
39: } while (0)
40:
41: #define SETUP_INPLACE() \
42: do { \
43: refmpn_copyi (got, data[i].src, data[i].size); \
44: got[data[i].size] = MAGIC; \
45: } while (0)
46:
47: #define VERIFY(name) \
48: do { \
49: verify (name, i, data[i].src, data[i].n, \
50: got_c, data[i].want_c, \
51: got, data[i].want, data[i].size); \
52: } while (0)
53:
54: typedef mp_limb_t (*mpn_aors_1_t)
55: _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
56: mpn_aors_1_t fudge _PROTO ((mpn_aors_1_t));
57:
58:
59: void
60: verify (const char *name, int i,
61: mp_srcptr src, mp_limb_t n,
62: mp_limb_t got_c, mp_limb_t want_c,
63: mp_srcptr got, mp_srcptr want, mp_size_t size)
64: {
65: if (got[size] != MAGIC)
66: {
67: printf ("Overwrite at %s i=%d\n", name, i);
68: abort ();
69: }
70:
71: if (got_c != want_c || ! refmpn_equal_anynail (got, want, size))
72: {
73: printf ("Wrong at %s i=%d size=%ld\n", name, i, size);
74: mpn_trace (" src", src, size);
75: mpn_trace (" n", &n, (mp_size_t) 1);
76: mpn_trace (" got", got, size);
77: mpn_trace (" want", want, size);
78: mpn_trace (" got c", &got_c, (mp_size_t) 1);
79: mpn_trace ("want c", &want_c, (mp_size_t) 1);
80: abort ();
81: }
82: }
83:
84:
85: void
86: check_add_1 (void)
87: {
88: static const struct {
89: mp_size_t size;
90: mp_limb_t n;
91: const mp_limb_t src[ASIZE];
92: mp_limb_t want_c;
93: const mp_limb_t want[ASIZE];
94: } data[] = {
95: { 1, 0, { 0 }, 0, { 0 } },
96: { 1, 0, { 1 }, 0, { 1 } },
97: { 1, 1, { 0 }, 0, { 1 } },
98: { 1, 0, { M }, 0, { M } },
99: { 1, M, { 0 }, 0, { M } },
100: { 1, 1, { 123 }, 0, { 124 } },
101:
102: { 1, 1, { M }, 1, { 0 } },
103: { 1, M, { 1 }, 1, { 0 } },
104: { 1, M, { M }, 1, { M-1 } },
105:
106: { 2, 0, { 0, 0 }, 0, { 0, 0 } },
107: { 2, 0, { 1, 0 }, 0, { 1, 0 } },
108: { 2, 1, { 0, 0 }, 0, { 1, 0 } },
109: { 2, 0, { M, 0 }, 0, { M, 0 } },
110: { 2, M, { 0, 0 }, 0, { M, 0 } },
111: { 2, 1, { M, 0 }, 0, { 0, 1 } },
112: { 2, M, { 1, 0 }, 0, { 0, 1 } },
113: { 2, M, { M, 0 }, 0, { M-1, 1 } },
114: { 2, M, { M, 0 }, 0, { M-1, 1 } },
115:
116: { 2, 1, { M, M }, 1, { 0, 0 } },
117: { 2, M, { 1, M }, 1, { 0, 0 } },
118: { 2, M, { M, M }, 1, { M-1, 0 } },
119: { 2, M, { M, M }, 1, { M-1, 0 } },
120:
121: { 3, 1, { M, M, M }, 1, { 0, 0, 0 } },
122: { 3, M, { 1, M, M }, 1, { 0, 0, 0 } },
123: { 3, M, { M, M, M }, 1, { M-1, 0, 0 } },
124: { 3, M, { M, M, M }, 1, { M-1, 0, 0 } },
125:
126: { 4, 1, { M, M, M, M }, 1, { 0, 0, 0, 0 } },
127: { 4, M, { 1, M, M, M }, 1, { 0, 0, 0, 0 } },
128: { 4, M, { M, M, M, M }, 1, { M-1, 0, 0, 0 } },
129: { 4, M, { M, M, M, M }, 1, { M-1, 0, 0, 0 } },
130:
131: { 4, M, { M, 0, M, M }, 0, { M-1, 1, M, M } },
132: { 4, M, { M, M-1, M, M }, 0, { M-1, M, M, M } },
133:
134: { 4, M, { M, M, 0, M }, 0, { M-1, 0, 1, M } },
135: { 4, M, { M, M, M-1, M }, 0, { M-1, 0, M, M } },
136: };
137:
138: mp_limb_t got[ASIZE];
139: mp_limb_t got_c;
140: int i;
141:
142: for (i = 0; i < numberof (data); i++)
143: {
144: SETUP ();
145: got_c = mpn_add_1 (got, data[i].src, data[i].size, data[i].n);
146: VERIFY ("check_add_1 (separate)");
147:
148: SETUP_INPLACE ();
149: got_c = mpn_add_1 (got, got, data[i].size, data[i].n);
150: VERIFY ("check_add_1 (in-place)");
151:
152: if (data[i].n == 1)
153: {
154: SETUP ();
155: got_c = mpn_add_1 (got, data[i].src, data[i].size, CNST_LIMB(1));
156: VERIFY ("check_add_1 (separate, const 1)");
157:
158: SETUP_INPLACE ();
159: got_c = mpn_add_1 (got, got, data[i].size, CNST_LIMB(1));
160: VERIFY ("check_add_1 (in-place, const 1)");
161: }
162:
163: /* Same again on functions, not inlines. */
164: SETUP ();
165: got_c = (*fudge(mpn_add_1)) (got, data[i].src, data[i].size, data[i].n);
166: VERIFY ("check_add_1 (function, separate)");
167:
168: SETUP_INPLACE ();
169: got_c = (*fudge(mpn_add_1)) (got, got, data[i].size, data[i].n);
170: VERIFY ("check_add_1 (function, in-place)");
171: }
172: }
173:
174: void
175: check_sub_1 (void)
176: {
177: static const struct {
178: mp_size_t size;
179: mp_limb_t n;
180: const mp_limb_t src[ASIZE];
181: mp_limb_t want_c;
182: const mp_limb_t want[ASIZE];
183: } data[] = {
184: { 1, 0, { 0 }, 0, { 0 } },
185: { 1, 0, { 1 }, 0, { 1 } },
186: { 1, 1, { 1 }, 0, { 0 } },
187: { 1, 0, { M }, 0, { M } },
188: { 1, 1, { M }, 0, { M-1 } },
189: { 1, 1, { 123 }, 0, { 122 } },
190:
191: { 1, 1, { 0 }, 1, { M } },
192: { 1, M, { 0 }, 1, { 1 } },
193:
194: { 2, 0, { 0, 0 }, 0, { 0, 0 } },
195: { 2, 0, { 1, 0 }, 0, { 1, 0 } },
196: { 2, 1, { 1, 0 }, 0, { 0, 0 } },
197: { 2, 0, { M, 0 }, 0, { M, 0 } },
198: { 2, 1, { M, 0 }, 0, { M-1, 0 } },
199: { 2, 1, { 123, 0 }, 0, { 122, 0 } },
200:
201: { 2, 1, { 0, 0 }, 1, { M, M } },
202: { 2, M, { 0, 0 }, 1, { 1, M } },
203:
204: { 3, 0, { 0, 0, 0 }, 0, { 0, 0, 0 } },
205: { 3, 0, { 123, 0, 0 }, 0, { 123, 0, 0 } },
206:
207: { 3, 1, { 0, 0, 0 }, 1, { M, M, M } },
208: { 3, M, { 0, 0, 0 }, 1, { 1, M, M } },
209:
210: { 4, 1, { 0, 0, 0, 0 }, 1, { M, M, M, M } },
211: { 4, M, { 0, 0, 0, 0 }, 1, { 1, M, M, M } },
212:
213: { 4, 1, { 0, 0, 1, 42 }, 0, { M, M, 0, 42 } },
214: { 4, M, { 0, 0, 123, 24 }, 0, { 1, M, 122, 24 } },
215: };
216:
217: mp_limb_t got[ASIZE];
218: mp_limb_t got_c;
219: int i;
220:
221: for (i = 0; i < numberof (data); i++)
222: {
223: SETUP ();
224: got_c = mpn_sub_1 (got, data[i].src, data[i].size, data[i].n);
225: VERIFY ("check_sub_1 (separate)");
226:
227: SETUP_INPLACE ();
228: got_c = mpn_sub_1 (got, got, data[i].size, data[i].n);
229: VERIFY ("check_sub_1 (in-place)");
230:
231: if (data[i].n == 1)
232: {
233: SETUP ();
234: got_c = mpn_sub_1 (got, data[i].src, data[i].size, CNST_LIMB(1));
235: VERIFY ("check_sub_1 (separate, const 1)");
236:
237: SETUP_INPLACE ();
238: got_c = mpn_sub_1 (got, got, data[i].size, CNST_LIMB(1));
239: VERIFY ("check_sub_1 (in-place, const 1)");
240: }
241:
242: /* Same again on functions, not inlines. */
243: SETUP ();
244: got_c = (*fudge(mpn_sub_1)) (got, data[i].src, data[i].size, data[i].n);
245: VERIFY ("check_sub_1 (function, separate)");
246:
247: SETUP_INPLACE ();
248: got_c = (*fudge(mpn_sub_1)) (got, got, data[i].size, data[i].n);
249: VERIFY ("check_sub_1 (function, in-place)");
250: }
251: }
252:
253: /* Try to prevent the optimizer inlining. */
254: mpn_aors_1_t
255: fudge (mpn_aors_1_t f)
256: {
257: return f;
258: }
259:
260: int
261: main (void)
262: {
263: tests_start ();
264: mp_trace_base = -16;
265:
266: check_add_1 ();
267: check_sub_1 ();
268:
269: tests_end ();
270: exit (0);
271: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to t-aors_1.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / tests / mpn