X-Git-Url: http://git.megacz.com/?a=blobdiff_plain;f=ghc%2Fruntime%2Fgmp%2Fmpz_mdiv.c;fp=ghc%2Fruntime%2Fgmp%2Fmpz_mdiv.c;h=bb8b9de86b08df8a88e3f3e553968f7ca31127a5;hb=e7d21ee4f8ac907665a7e170c71d59e13a01da09;hp=0000000000000000000000000000000000000000;hpb=e48474bff05e6cfb506660420f025f694c870d38;p=ghc-hetmet.git diff --git a/ghc/runtime/gmp/mpz_mdiv.c b/ghc/runtime/gmp/mpz_mdiv.c new file mode 100644 index 0000000..bb8b9de --- /dev/null +++ b/ghc/runtime/gmp/mpz_mdiv.c @@ -0,0 +1,52 @@ +/* mpz_mdiv -- Mathematical DIVision and MODulo, i.e. division that rounds + the quotient towards -infinity. + +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" + +void +#ifdef __STDC__ +mpz_mdiv (MP_INT *quot, + const MP_INT *dividend, const MP_INT *divisor) +#else +mpz_mdiv (quot, dividend, divisor) + MP_INT *quot; + const MP_INT *dividend; + const MP_INT *divisor; +#endif +{ + if ((dividend->size ^ divisor->size) >= 0) + { + /* When the dividend and the divisor has same sign, this function + gives same result as mpz_div. */ + mpz_div (quot, dividend, divisor); + } + else + { + MP_INT rem; + + MPZ_TMP_INIT (&rem, 1 + ABS (dividend->size)); + + mpz_divmod (quot, &rem, dividend, divisor); + if (rem.size != 0) + mpz_sub_ui (quot, quot, 1); + } +}