diff options
author | markm <markm@FreeBSD.org> | 1995-11-12 14:40:41 +0000 |
---|---|---|
committer | markm <markm@FreeBSD.org> | 1995-11-12 14:40:41 +0000 |
commit | d938bb78837327b46dc34873a6eb60c002605808 (patch) | |
tree | d720c213b0c2f268c17d14276e809e876bd37dc8 /gnu/lib/libgmp/tests/tst-mdm_ui.c | |
parent | 7e662b2efdf7b11bed05bdd944746d7b7d0e56af (diff) | |
download | FreeBSD-src-d938bb78837327b46dc34873a6eb60c002605808.zip FreeBSD-src-d938bb78837327b46dc34873a6eb60c002605808.tar.gz |
GNU MP (Multiprecision) library. This is needed by secure RPC (being
done by Bill Paul) and various other BSD programs.
Obtained from:FSF
Diffstat (limited to 'gnu/lib/libgmp/tests/tst-mdm_ui.c')
-rw-r--r-- | gnu/lib/libgmp/tests/tst-mdm_ui.c | 116 |
1 files changed, 116 insertions, 0 deletions
diff --git a/gnu/lib/libgmp/tests/tst-mdm_ui.c b/gnu/lib/libgmp/tests/tst-mdm_ui.c new file mode 100644 index 0000000..26befbe --- /dev/null +++ b/gnu/lib/libgmp/tests/tst-mdm_ui.c @@ -0,0 +1,116 @@ +/* Test mpz_abs, mpz_add, mpz_cmp. mpz_cmp_si, mpz_cmp_ui, mpz_divmod, + mpz_mdiv_ui, mpz_mdivmod_ui, mpz_mmod_ui, mpz_mul, mpz_mul_ui. + +Copyright (C) 1993 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 General Public License as published by +the Free Software Foundation; either version 2, 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 General Public License for more details. + +You should have received a copy of the GNU General Public License +along with the GNU MP Library; see the file COPYING. If not, write to +the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. */ + +#include <stdio.h> +#include "gmp.h" +#include "gmp-impl.h" +#include "urandom.h" + +void debug_mp (); + +#ifndef SIZE +#define SIZE 8 +#endif + +main (argc, argv) + int argc; + char **argv; +{ + MP_INT dividend; + MP_INT quotient, remainder; + MP_INT quotient2, remainder2; + MP_INT temp; + mp_size dividend_size; + mp_limb divisor; + int i; + int reps = 100000; + + if (argc == 2) + reps = atoi (argv[1]); + + mpz_init (÷nd); + mpz_init ("ient); + mpz_init (&remainder); + mpz_init ("ient2); + mpz_init (&remainder2); + mpz_init (&temp); + + for (i = 0; i < reps; i++) + { + dividend_size = urandom () % SIZE - SIZE/2; + mpz_random2 (÷nd, dividend_size); + + divisor = urandom (); + if (divisor == 0) + continue; + + mpz_mdivmod_ui ("ient, &remainder, ÷nd, divisor); + mpz_mdiv_ui ("ient2, ÷nd, divisor); + mpz_mmod_ui (&remainder2, ÷nd, divisor); + + /* First determine that the quotients and remainders computed + with different functions are equal. */ + if (mpz_cmp ("ient, "ient2) != 0) + dump_abort (÷nd, divisor); + if (mpz_cmp (&remainder, &remainder2) != 0) + dump_abort (÷nd, divisor); + + /* Check if the sign of the quotient is correct. */ + if (mpz_cmp_ui ("ient, 0) != 0) + if ((mpz_cmp_ui ("ient, 0) < 0) + != (mpz_cmp_ui (÷nd, 0) < 0)) + dump_abort (÷nd, divisor); + + /* Check if the remainder has the same sign as the divisor + (quotient rounded towards minus infinity). */ + if (mpz_cmp_ui (&remainder, 0) != 0) + if (mpz_cmp_ui (&remainder, 0) < 0) + dump_abort (÷nd, divisor); + + mpz_mul_ui (&temp, "ient, divisor); + mpz_add (&temp, &temp, &remainder); + if (mpz_cmp (&temp, ÷nd) != 0) + dump_abort (÷nd, divisor); + + mpz_abs (&remainder, &remainder); + if (mpz_cmp_ui (&remainder, divisor) >= 0) + dump_abort (÷nd, divisor); + } + + exit (0); +} + +dump_abort (dividend, divisor) + MP_INT *dividend; + mp_limb divisor; +{ + fprintf (stderr, "ERROR\n"); + fprintf (stderr, "dividend = "); debug_mp (dividend, -16); + fprintf (stderr, "divisor = %lX\n", divisor); + abort(); +} + +void +debug_mp (x, base) + MP_INT *x; +{ + mpz_out_str (stderr, base, x); fputc ('\n', stderr); +} |