diff options
Diffstat (limited to 'gnu/lib/libgmp/mpz_mdm_ui.c')
-rw-r--r-- | gnu/lib/libgmp/mpz_mdm_ui.c | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/gnu/lib/libgmp/mpz_mdm_ui.c b/gnu/lib/libgmp/mpz_mdm_ui.c new file mode 100644 index 0000000..cf41912 --- /dev/null +++ b/gnu/lib/libgmp/mpz_mdm_ui.c @@ -0,0 +1,58 @@ +/* mpz_mdivmod -- Mathematical DIVision and MODulo, i.e. division that rounds + the quotient towards -infinity, and with the remainder non-negative. + +Copyright (C) 1991, 1992 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 "gmp.h" +#include "gmp-impl.h" + +unsigned long int +#ifdef __STDC__ +mpz_mdivmod_ui (MP_INT *quot, MP_INT *rem, + const MP_INT *dividend, unsigned long int divisor) +#else +mpz_mdivmod_ui (quot, rem, dividend, divisor) + MP_INT *quot; + MP_INT *rem; + const MP_INT *dividend; + unsigned long int divisor; +#endif +{ + MP_INT temp_rem; /* N.B.: lives until function returns! */ + + /* If the user doesn't want the remainder to be stored in an integer + object, allocate a scratch variable for it. */ + if (rem == NULL) + { + MPZ_TMP_INIT (&temp_rem, 1 + ABS (dividend->size)); + rem = &temp_rem; + } + + mpz_divmod_ui (quot, rem, dividend, divisor); + + if (rem->size < 0) + { + mpz_sub_ui (quot, quot, 1); + mpz_add_ui (rem, rem, divisor); + } + + if (rem->size == 0) + return 0; + return rem->d[0]; +} |