From 11da511c784eca003deb90c23570f0873954e0de Mon Sep 17 00:00:00 2001 From: Duncan Wilkie Date: Sat, 18 Nov 2023 06:11:09 -0600 Subject: Initial commit. --- gmp-6.3.0/mini-gmp/tests/t-div.c | 258 +++++++++++++++++++++++++++++++++++++++ 1 file changed, 258 insertions(+) create mode 100644 gmp-6.3.0/mini-gmp/tests/t-div.c (limited to 'gmp-6.3.0/mini-gmp/tests/t-div.c') diff --git a/gmp-6.3.0/mini-gmp/tests/t-div.c b/gmp-6.3.0/mini-gmp/tests/t-div.c new file mode 100644 index 0000000..be2f6b1 --- /dev/null +++ b/gmp-6.3.0/mini-gmp/tests/t-div.c @@ -0,0 +1,258 @@ +/* + +Copyright 2012, 2013, 2018 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 +#include +#include + +#include "testutils.h" + +#define MAXBITS 400 +#define COUNT 10000 + +typedef void div_qr_func (mpz_t, mpz_t, const mpz_t, const mpz_t); +typedef unsigned long div_qr_ui_func (mpz_t, mpz_t, const mpz_t, unsigned long); +typedef void div_func (mpz_t, const mpz_t, const mpz_t); +typedef unsigned long div_x_ui_func (mpz_t, const mpz_t, unsigned long); +typedef unsigned long div_ui_func (const mpz_t, unsigned long); + +void +testmain (int argc, char **argv) +{ + unsigned i; + mpz_t a, b, q, r, rq, rr; + int div_p; + + mpz_init (a); + mpz_init (b); + mpz_init (r); + mpz_init (q); + mpz_init (rr); + mpz_init (rq); + + for (i = 0; i < COUNT; i++) + { + unsigned j; + for (j = 0; j < 3; j++) + { + static const enum hex_random_op ops[3] = { OP_CDIV, OP_FDIV, OP_TDIV }; + static const char name[3] = { 'c', 'f', 't'}; + static div_qr_func * const div_qr [3] = + { + mpz_cdiv_qr, mpz_fdiv_qr, mpz_tdiv_qr + }; + static div_qr_ui_func * const div_qr_ui[3] = + { + mpz_cdiv_qr_ui, mpz_fdiv_qr_ui, mpz_tdiv_qr_ui + }; + static div_func * const div_q [3] = + { + mpz_cdiv_q, mpz_fdiv_q, mpz_tdiv_q + }; + static div_x_ui_func * const div_q_ui[3] = + { + mpz_cdiv_q_ui, mpz_fdiv_q_ui, mpz_tdiv_q_ui + }; + static div_func * const div_r [3] = + { + mpz_cdiv_r, mpz_fdiv_r, mpz_tdiv_r + }; + static div_x_ui_func * const div_r_ui[3] = + { + mpz_cdiv_r_ui, mpz_fdiv_r_ui, mpz_tdiv_r_ui + }; + static div_ui_func * const div_ui[3] = + { + mpz_cdiv_ui, mpz_fdiv_ui, mpz_tdiv_ui + }; + + mini_random_op4 (ops[j], MAXBITS, a, b, rq, rr); + div_qr[j] (q, r, a, b); + if (mpz_cmp (r, rr) || mpz_cmp (q, rq)) + { + fprintf (stderr, "mpz_%cdiv_qr failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + dump ("rref", rr); + dump ("q ", q); + dump ("qref", rq); + abort (); + } + mpz_set_si (q, -5); + div_q[j] (q, a, b); + if (mpz_cmp (q, rq)) + { + fprintf (stderr, "mpz_%cdiv_q failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + dump ("q ", q); + dump ("qref", rq); + abort (); + } + mpz_set_ui (r, ~5); + div_r[j] (r, a, b); + if (mpz_cmp (r, rr)) + { + fprintf (stderr, "mpz_%cdiv_r failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + dump ("rref", rr); + abort (); + } + + if (j == 0) /* do this once, not for all roundings */ + { + div_p = mpz_divisible_p (a, b); + if ((mpz_sgn (r) == 0) ^ (div_p != 0)) + { + fprintf (stderr, "mpz_divisible_p failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + abort (); + } + } + + mpz_set_si (r, -6); + if (j == 0 && mpz_sgn (b) < 0) /* ceil, negative divisor */ + { + mpz_mod (r, a, b); + if (mpz_cmp (r, rr)) + { + fprintf (stderr, "mpz_mod failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + dump ("rref", rr); + abort (); + } + } + + if (j == 1 && mpz_sgn (b) > 0) /* floor, positive divisor */ + { + mpz_mod (r, a, b); + if (mpz_cmp (r, rr)) + { + fprintf (stderr, "mpz_mod failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + dump ("rref", rr); + abort (); + } + } + + if (mpz_fits_ulong_p (b)) + { + unsigned long rl; + + mpz_set_si (r, -7); + mpz_set_ui (q, ~7); + rl = div_qr_ui[j] (q, r, a, mpz_get_ui (b)); + if (rl != mpz_get_ui (rr) + || mpz_cmp (r, rr) || mpz_cmp (q, rq)) + { + fprintf (stderr, "mpz_%cdiv_qr_ui failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + fprintf(stderr, "rl = %lx\n", rl); + dump ("r ", r); + dump ("rref", rr); + dump ("q ", q); + dump ("qref", rq); + abort (); + } + + mpz_set_si (q, 3); + rl = div_q_ui[j] (q, a, mpz_get_ui (b)); + if (rl != mpz_get_ui (rr) || mpz_cmp (q, rq)) + { + fprintf (stderr, "mpz_%cdiv_q_ui failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + fprintf(stderr, "rl = %lx\n", rl); + dump ("rref", rr); + dump ("q ", q); + dump ("qref", rq); + abort (); + } + + mpz_set_ui (r, 7); + rl = div_r_ui[j] (r, a, mpz_get_ui (b)); + if (rl != mpz_get_ui (rr) || mpz_cmp (r, rr)) + { + fprintf (stderr, "mpz_%cdiv_r_ui failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + fprintf(stderr, "rl = %lx\n", rl); + dump ("r ", r); + dump ("rref", rr); + abort (); + } + + rl = div_ui[j] (a, mpz_get_ui (b)); + if (rl != mpz_get_ui (rr)) + { + fprintf (stderr, "mpz_%cdiv_ui failed:\n", name[j]); + dump ("a", a); + dump ("b", b); + fprintf(stderr, "rl = %lx\n", rl); + dump ("rref", rr); + abort (); + } + + if (j == 0) /* do this once, not for all roundings */ + { + div_p = mpz_divisible_ui_p (a, mpz_get_ui (b)); + if ((mpz_sgn (r) == 0) ^ (div_p != 0)) + { + fprintf (stderr, "mpz_divisible_ui_p failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + abort (); + } + } + + if (j == 1) /* floor */ + { + mpz_set_si (r, -2); + mpz_mod_ui (r, a, mpz_get_ui (b)); + if (mpz_cmp (r, rr)) + { + fprintf (stderr, "mpz_mod failed:\n"); + dump ("a", a); + dump ("b", b); + dump ("r ", r); + dump ("rref", rr); + abort (); + } + } + } + } + } + mpz_clear (a); + mpz_clear (b); + mpz_clear (r); + mpz_clear (q); + mpz_clear (rr); + mpz_clear (rq); +} -- cgit v1.2.3