Annotation of OpenXM_contrib/gmp/tests/mpz/t-aorsmul.c, Revision 1.1
1.1 ! ohara 1: /* Test mpz_addmul, mpz_addmul_ui, mpz_submul, mpz_submul_ui.
! 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:
! 24: #include <stdio.h>
! 25: #include <stdlib.h>
! 26: #include <string.h>
! 27:
! 28: #include "gmp.h"
! 29: #include "gmp-impl.h"
! 30: #include "tests.h"
! 31:
! 32:
! 33: void
! 34: check_one (mpz_srcptr w, mpz_srcptr x, mpz_srcptr y)
! 35: {
! 36: mpz_t want, got;
! 37:
! 38: mpz_init (want);
! 39: mpz_init (got);
! 40:
! 41: mpz_mul (want, x, y);
! 42: mpz_add (want, w, want);
! 43: mpz_set (got, w);
! 44: mpz_addmul (got, x, y);
! 45: MPZ_CHECK_FORMAT (got);
! 46: if (mpz_cmp (want, got) != 0)
! 47: {
! 48: printf ("mpz_addmul fail\n");
! 49: fail:
! 50: mpz_trace ("w", w);
! 51: mpz_trace ("x", x);
! 52: mpz_trace ("y", y);
! 53: mpz_trace ("want", want);
! 54: mpz_trace ("got ", got);
! 55: abort ();
! 56: }
! 57:
! 58: mpz_mul (want, x, y);
! 59: mpz_sub (want, w, want);
! 60: mpz_set (got, w);
! 61: mpz_submul (got, x, y);
! 62: MPZ_CHECK_FORMAT (got);
! 63: if (mpz_cmp (want, got) != 0)
! 64: {
! 65: printf ("mpz_submul fail\n");
! 66: goto fail;
! 67: }
! 68:
! 69: mpz_clear (want);
! 70: mpz_clear (got);
! 71: }
! 72:
! 73:
! 74: void
! 75: check_one_ui (mpz_ptr w, mpz_ptr x, unsigned long y)
! 76: {
! 77: mpz_t want, got;
! 78:
! 79: mpz_init (want);
! 80: mpz_init (got);
! 81:
! 82: mpz_mul_ui (want, x, (unsigned long) y);
! 83: mpz_add (want, w, want);
! 84: mpz_set (got, w);
! 85: mpz_addmul_ui (got, x, (unsigned long) y);
! 86: MPZ_CHECK_FORMAT (got);
! 87: if (mpz_cmp (want, got) != 0)
! 88: {
! 89: printf ("mpz_addmul_ui fail\n");
! 90: fail:
! 91: mpz_trace ("w", w);
! 92: mpz_trace ("x", x);
! 93: printf ("y=0x%lX %lu\n", y, y);
! 94: mpz_trace ("want", want);
! 95: mpz_trace ("got ", got);
! 96: abort ();
! 97: }
! 98:
! 99: mpz_mul_ui (want, x, y);
! 100: mpz_sub (want, w, want);
! 101: mpz_set (got, w);
! 102: mpz_submul_ui (got, x, y);
! 103: MPZ_CHECK_FORMAT (got);
! 104: if (mpz_cmp (want, got) != 0)
! 105: {
! 106: printf ("mpz_submul_ui fail\n");
! 107: goto fail;
! 108: }
! 109:
! 110: mpz_clear (want);
! 111: mpz_clear (got);
! 112: }
! 113:
! 114:
! 115: void
! 116: check_all (mpz_ptr w, mpz_ptr x, mpz_ptr y)
! 117: {
! 118: int swap, wneg, xneg, yneg;
! 119:
! 120: MPZ_CHECK_FORMAT (w);
! 121: MPZ_CHECK_FORMAT (x);
! 122: MPZ_CHECK_FORMAT (y);
! 123:
! 124: for (swap = 0; swap < 2; swap++)
! 125: {
! 126: for (wneg = 0; wneg < 2; wneg++)
! 127: {
! 128: for (xneg = 0; xneg < 2; xneg++)
! 129: {
! 130: for (yneg = 0; yneg < 2; yneg++)
! 131: {
! 132: check_one (w, x, y);
! 133:
! 134: if (mpz_fits_ulong_p (y))
! 135: check_one_ui (w, x, mpz_get_ui (y));
! 136:
! 137: mpz_neg (y, y);
! 138: }
! 139: mpz_neg (x, x);
! 140: }
! 141: mpz_neg (w, w);
! 142: }
! 143: mpz_swap (x, y);
! 144: }
! 145: }
! 146:
! 147: void
! 148: check_data (void)
! 149: {
! 150: static const struct {
! 151: mp_limb_t w[6];
! 152: mp_limb_t x[6];
! 153: mp_limb_t y[6];
! 154:
! 155: } data[] = {
! 156:
! 157: #define M GMP_NUMB_MAX
! 158:
! 159: /* reducing to zero */
! 160: { { 1 }, { 1 }, { 1 } },
! 161: { { 2 }, { 1 }, { 2 } },
! 162: { { 0,1 }, { 0,1 }, { 1 } },
! 163:
! 164: /* reducing to 1 */
! 165: { { 0,1 }, { M }, { 1 } },
! 166: { { 0,0,1 }, { M,M }, { 1 } },
! 167: { { 0,0,0,1 }, { M,M,M }, { 1 } },
! 168: { { 0,0,0,0,1 }, { M,M,M,M }, { 1 } },
! 169:
! 170: /* reducing to -1 */
! 171: { { M }, { 0,1 }, { 1 } },
! 172: { { M,M }, { 0,0,1 }, { 1 } },
! 173: { { M,M,M }, { 0,0,0,1 }, { 1 } },
! 174: { { M,M,M,M }, { 0,0,0,0,1 }, { 1 } },
! 175:
! 176: /* carry out of addmul */
! 177: { { M }, { 1 }, { 1 } },
! 178: { { M,M }, { 1 }, { 1 } },
! 179: { { M,M,M }, { 1 }, { 1 } },
! 180:
! 181: /* borrow from submul */
! 182: { { 0,1 }, { 1 }, { 1 } },
! 183: { { 0,0,1 }, { 1 }, { 1 } },
! 184: { { 0,0,0,1 }, { 1 }, { 1 } },
! 185:
! 186: /* borrow from submul */
! 187: { { 0,0,1 }, { 0,1 }, { 1 } },
! 188: { { 0,0,0,1 }, { 0,1 }, { 1 } },
! 189: { { 0,0,0,0,1 }, { 0,1 }, { 1 } },
! 190:
! 191: /* more borrow from submul */
! 192: { { M }, { 0,1 }, { 1 } },
! 193: { { M }, { 0,0,1 }, { 1 } },
! 194: { { M }, { 0,0,0,1 }, { 1 } },
! 195: { { M }, { 0,0,0,0,1 }, { 1 } },
! 196:
! 197: /* big borrow from submul */
! 198: { { 0,0,1 }, { M,M }, { M } },
! 199: { { 0,0,0,1 }, { M,M }, { M } },
! 200: { { 0,0,0,0,1 }, { M,M }, { M } },
! 201:
! 202: /* small w */
! 203: { { 0,1 }, { M,M }, { M } },
! 204: { { 0,1 }, { M,M,M }, { M } },
! 205: { { 0,1 }, { M,M,M,M }, { M } },
! 206: { { 0,1 }, { M,M,M,M,M }, { M } },
! 207: };
! 208:
! 209: mpz_t w, x, y;
! 210: int i;
! 211:
! 212: mpz_init (w);
! 213: mpz_init (x);
! 214: mpz_init (y);
! 215:
! 216: for (i = 0; i < numberof (data); i++)
! 217: {
! 218: mpz_set_n (w, data[i].w, numberof(data[i].w));
! 219: mpz_set_n (x, data[i].x, numberof(data[i].x));
! 220: mpz_set_n (y, data[i].y, numberof(data[i].y));
! 221: check_all (w, x, y);
! 222: }
! 223:
! 224: mpz_clear (w);
! 225: mpz_clear (x);
! 226: mpz_clear (y);
! 227: }
! 228:
! 229:
! 230: void
! 231: check_random (int argc, char *argv[])
! 232: {
! 233: gmp_randstate_ptr rands = RANDS;
! 234: mpz_t w, x, y;
! 235: int i, reps = 2000;
! 236:
! 237: mpz_init (w);
! 238: mpz_init (x);
! 239: mpz_init (y);
! 240:
! 241: if (argc == 2)
! 242: reps = atoi (argv[1]);
! 243:
! 244: for (i = 0; i < reps; i++)
! 245: {
! 246: mpz_errandomb (w, rands, 5*BITS_PER_MP_LIMB);
! 247: mpz_errandomb (x, rands, 5*BITS_PER_MP_LIMB);
! 248: mpz_errandomb (y, rands, 5*BITS_PER_MP_LIMB);
! 249: check_all (w, x, y);
! 250:
! 251: mpz_errandomb (w, rands, 5*BITS_PER_MP_LIMB);
! 252: mpz_errandomb (x, rands, 5*BITS_PER_MP_LIMB);
! 253: mpz_errandomb (y, rands, BITS_PER_ULONG);
! 254: check_all (w, x, y);
! 255: }
! 256:
! 257: mpz_clear (w);
! 258: mpz_clear (x);
! 259: mpz_clear (y);
! 260: }
! 261:
! 262:
! 263: int
! 264: main (int argc, char *argv[])
! 265: {
! 266: tests_start ();
! 267: mp_trace_base = -16;
! 268:
! 269: check_data ();
! 270: check_random (argc, argv);
! 271:
! 272: tests_end ();
! 273: exit (0);
! 274: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to t-aorsmul.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