diff options
Diffstat (limited to 'vendor/gmp-6.3.0/mini-gmp/tests/t-logops.c')
-rw-r--r-- | vendor/gmp-6.3.0/mini-gmp/tests/t-logops.c | 112 |
1 files changed, 112 insertions, 0 deletions
diff --git a/vendor/gmp-6.3.0/mini-gmp/tests/t-logops.c b/vendor/gmp-6.3.0/mini-gmp/tests/t-logops.c new file mode 100644 index 0000000..b4e0239 --- /dev/null +++ b/vendor/gmp-6.3.0/mini-gmp/tests/t-logops.c @@ -0,0 +1,112 @@ +/* + +Copyright 2012, 2013 Free Software Foundation, Inc. + +This file is part of the GNU MP Library test suite. + +The GNU MP Library test suite is free software; you can redistribute it +and/or modify it under the terms of the GNU 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 test suite 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 General +Public License for more details. + +You should have received a copy of the GNU General Public License along with +the GNU MP Library test suite. If not, see https://www.gnu.org/licenses/. */ + +#include <limits.h> +#include <stdlib.h> +#include <stdio.h> + +#include "testutils.h" + +#define MAXBITS 400 +#define COUNT 10000 + +void +testlogops (int count) +{ + int i; + mpz_t a, b, res, ref; + mp_bitcnt_t c; + + mpz_init (a); + mpz_init (b); + mpz_init (res); + mpz_init (ref); + + for (i = 0; i < count; i++) + { + mini_random_op3 (OP_AND, MAXBITS, a, b, ref); + mpz_and (res, a, b); + if (mpz_cmp (res, ref)) + { + fprintf (stderr, "mpz_and failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r", res); + dump ("ref", ref); + abort (); + } + + mini_random_op3 (OP_IOR, MAXBITS, a, b, ref); + mpz_ior (res, a, b); + if (mpz_cmp (res, ref)) + { + fprintf (stderr, "mpz_ior failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r", res); + dump ("ref", ref); + abort (); + } + + mini_random_op3 (OP_XOR, MAXBITS, a, b, ref); + mpz_xor (res, a, b); + if (mpz_cmp (res, ref)) + { + fprintf (stderr, "mpz_xor failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r", res); + dump ("ref", ref); + abort (); + } + + if (i % 8) { + c = 0; + mpz_mul_2exp (res, res, i % 8); + } else if (mpz_sgn (res) >= 0) { + c = mpz_odd_p (res) != 0; + mpz_tdiv_q_2exp (res, res, 1); + } else { + c = (~ (mp_bitcnt_t) 0) - 3; + mpz_set_ui (res, 11 << ((i >> 3)%4)); /* set 3 bits */ + } + + if (mpz_popcount (res) + c != mpz_hamdist (a, b)) + { + fprintf (stderr, "mpz_popcount(r) + %lu and mpz_hamdist(a,b) differ:\n", c); + dump ("a", a); + dump ("b", b); + dump ("r", res); + fprintf (stderr, "mpz_popcount(r) = %lu:\n", mpz_popcount (res)); + fprintf (stderr, "mpz_hamdist(a,b) = %lu:\n", mpz_hamdist (a, b)); + abort (); + } + } + mpz_clear (a); + mpz_clear (b); + mpz_clear (res); + mpz_clear (ref); +} + +void +testmain (int argc, char **argv) +{ + testhalves (COUNT*2/3, testlogops); + testlogops (COUNT/3); +} |