t-cong.c
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:5k
- /* test mpz_congruent_p and mpz_congruent_ui_p
- Copyright 2001, 2002 Free Software Foundation, Inc.
- This file is part of the GNU MP Library.
- The GNU MP Library is free software; you can redistribute it and/or modify
- it under the terms of the GNU Lesser 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 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 Lesser General Public
- License for more details.
- You should have received a copy of the GNU Lesser General Public License
- along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */
- #include <stdio.h>
- #include <stdlib.h>
- #include "gmp.h"
- #include "gmp-impl.h"
- #include "tests.h"
- void
- check_one (mpz_srcptr a, mpz_srcptr c, mpz_srcptr d, int want)
- {
- int got;
- int swap;
- for (swap = 0; swap <= 1; swap++)
- {
- got = (mpz_congruent_p (a, c, d) != 0);
- if (want != got)
- {
- printf ("mpz_congruent_p wrongn");
- printf (" expected %d got %dn", want, got);
- mpz_trace (" a", a);
- mpz_trace (" c", c);
- mpz_trace (" d", d);
- mp_trace_base = -16;
- mpz_trace (" a", a);
- mpz_trace (" c", c);
- mpz_trace (" d", d);
- abort ();
- }
- if (mpz_fits_ulong_p (c) && mpz_fits_ulong_p (d))
- {
- unsigned long uc = mpz_get_ui (c);
- unsigned long ud = mpz_get_ui (d);
- got = (mpz_congruent_ui_p (a, uc, ud) != 0);
- if (want != got)
- {
- printf ("mpz_congruent_ui_p wrongn");
- printf (" expected %d got %dn", want, got);
- mpz_trace (" a", a);
- printf (" c=%lun", uc);
- printf (" d=%lun", ud);
- mp_trace_base = -16;
- mpz_trace (" a", a);
- printf (" c=0x%lXn", uc);
- printf (" d=0x%lXn", ud);
- abort ();
- }
- }
- MPZ_SRCPTR_SWAP (a, c);
- }
- }
- void
- check_data (void)
- {
- static const struct {
- const char *a;
- const char *c;
- const char *d;
- int want;
- } data[] = {
- /* anything congruent mod 1 */
- { "0", "0", "1", 1 },
- { "1", "0", "1", 1 },
- { "0", "1", "1", 1 },
- { "123", "456", "1", 1 },
- { "0x123456789123456789", "0x987654321987654321", "1", 1 },
- /* csize==1, dsize==2 changing to 1 after stripping 2s */
- { "0x3333333333333333", "0x33333333",
- "0x180000000", 1 },
- { "0x33333333333333333333333333333333", "0x3333333333333333",
- "0x18000000000000000", 1 },
- /* another dsize==2 becoming 1, with opposite signs this time */
- { "0x444444441",
- "-0x22222221F",
- "0x333333330", 1 },
- { "0x44444444444444441",
- "-0x2222222222222221F",
- "0x33333333333333330", 1 },
- };
- mpz_t a, c, d;
- int i;
- mpz_init (a);
- mpz_init (c);
- mpz_init (d);
- for (i = 0; i < numberof (data); i++)
- {
- mpz_set_str_or_abort (a, data[i].a, 0);
- mpz_set_str_or_abort (c, data[i].c, 0);
- mpz_set_str_or_abort (d, data[i].d, 0);
- check_one (a, c, d, data[i].want);
- }
- mpz_clear (a);
- mpz_clear (c);
- mpz_clear (d);
- }
- void
- check_random (int argc, char *argv[])
- {
- gmp_randstate_ptr rands = RANDS;
- mpz_t a, c, d, ra, rc;
- int i;
- int want;
- int reps = 50000;
- if (argc >= 2)
- reps = atoi (argv[1]);
- mpz_init (a);
- mpz_init (c);
- mpz_init (d);
- mpz_init (ra);
- mpz_init (rc);
- for (i = 0; i < reps; i++)
- {
- mpz_errandomb (a, rands, 8*GMP_LIMB_BITS);
- MPZ_CHECK_FORMAT (a);
- mpz_errandomb (c, rands, 8*GMP_LIMB_BITS);
- MPZ_CHECK_FORMAT (c);
- mpz_errandomb_nonzero (d, rands, 8*GMP_LIMB_BITS);
- mpz_negrandom (a, rands);
- MPZ_CHECK_FORMAT (a);
- mpz_negrandom (c, rands);
- MPZ_CHECK_FORMAT (c);
- mpz_negrandom (d, rands);
- mpz_fdiv_r (ra, a, d);
- mpz_fdiv_r (rc, c, d);
- want = (mpz_cmp (ra, rc) == 0);
- check_one (a, c, d, want);
- mpz_sub (ra, ra, rc);
- mpz_sub (a, a, ra);
- MPZ_CHECK_FORMAT (a);
- check_one (a, c, d, 1);
- if (! mpz_pow2abs_p (d))
- {
- refmpz_combit (a, urandom() % (8*GMP_LIMB_BITS));
- check_one (a, c, d, 0);
- }
- }
- mpz_clear (a);
- mpz_clear (c);
- mpz_clear (d);
- mpz_clear (ra);
- mpz_clear (rc);
- }
- int
- main (int argc, char *argv[])
- {
- tests_start ();
- check_data ();
- check_random (argc, argv);
- tests_end ();
- exit (0);
- }