1 1.1 mrg /* mpz_cdiv_qr -- Division rounding the quotient towards +infinity. The 2 1.1 mrg remainder gets the opposite sign as the denominator. 3 1.1 mrg 4 1.1.1.3 mrg Copyright 1994-1996, 2000, 2001, 2005, 2012 Free Software Foundation, Inc. 5 1.1 mrg 6 1.1 mrg This file is part of the GNU MP Library. 7 1.1 mrg 8 1.1 mrg The GNU MP Library is free software; you can redistribute it and/or modify 9 1.1.1.3 mrg it under the terms of either: 10 1.1.1.3 mrg 11 1.1.1.3 mrg * the GNU Lesser General Public License as published by the Free 12 1.1.1.3 mrg Software Foundation; either version 3 of the License, or (at your 13 1.1.1.3 mrg option) any later version. 14 1.1.1.3 mrg 15 1.1.1.3 mrg or 16 1.1.1.3 mrg 17 1.1.1.3 mrg * the GNU General Public License as published by the Free Software 18 1.1.1.3 mrg Foundation; either version 2 of the License, or (at your option) any 19 1.1.1.3 mrg later version. 20 1.1.1.3 mrg 21 1.1.1.3 mrg or both in parallel, as here. 22 1.1 mrg 23 1.1 mrg The GNU MP Library is distributed in the hope that it will be useful, but 24 1.1 mrg WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 1.1.1.3 mrg or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 26 1.1.1.3 mrg for more details. 27 1.1 mrg 28 1.1.1.3 mrg You should have received copies of the GNU General Public License and the 29 1.1.1.3 mrg GNU Lesser General Public License along with the GNU MP Library. If not, 30 1.1.1.3 mrg see https://www.gnu.org/licenses/. */ 31 1.1 mrg 32 1.1 mrg #include "gmp-impl.h" 33 1.1 mrg 34 1.1 mrg void 35 1.1 mrg mpz_cdiv_qr (mpz_ptr quot, mpz_ptr rem, mpz_srcptr dividend, mpz_srcptr divisor) 36 1.1 mrg { 37 1.1.1.2 mrg mp_size_t divisor_size = SIZ (divisor); 38 1.1 mrg mp_size_t xsize; 39 1.1 mrg mpz_t temp_divisor; /* N.B.: lives until function returns! */ 40 1.1 mrg TMP_DECL; 41 1.1 mrg 42 1.1 mrg TMP_MARK; 43 1.1 mrg 44 1.1 mrg /* We need the original value of the divisor after the quotient and 45 1.1 mrg remainder have been preliminary calculated. We have to copy it to 46 1.1 mrg temporary space if it's the same variable as either QUOT or REM. */ 47 1.1 mrg if (quot == divisor || rem == divisor) 48 1.1 mrg { 49 1.1 mrg MPZ_TMP_INIT (temp_divisor, ABS (divisor_size)); 50 1.1 mrg mpz_set (temp_divisor, divisor); 51 1.1 mrg divisor = temp_divisor; 52 1.1 mrg } 53 1.1 mrg 54 1.1.1.2 mrg xsize = SIZ (dividend) ^ divisor_size;; 55 1.1 mrg mpz_tdiv_qr (quot, rem, dividend, divisor); 56 1.1 mrg 57 1.1.1.2 mrg if (xsize >= 0 && SIZ (rem) != 0) 58 1.1 mrg { 59 1.1 mrg mpz_add_ui (quot, quot, 1L); 60 1.1 mrg mpz_sub (rem, rem, divisor); 61 1.1 mrg } 62 1.1 mrg 63 1.1 mrg TMP_FREE; 64 1.1 mrg } 65
Indexes created Sat Feb 28 05:31:39 UTC 2026