|
Packit |
5c3484 |
/* mpz_divexact_gcd -- exact division optimized for GCDs.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
THE FUNCTIONS IN THIS FILE ARE FOR INTERNAL USE AND ARE ALMOST CERTAIN TO
|
|
Packit |
5c3484 |
BE SUBJECT TO INCOMPATIBLE CHANGES IN FUTURE GNU MP RELEASES.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Copyright 2000, 2005, 2011, 2012 Free Software Foundation, Inc.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
This file is part of the GNU MP Library.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The GNU MP Library is free software; you can redistribute it and/or modify
|
|
Packit |
5c3484 |
it under the terms of either:
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
* the GNU Lesser General Public License as published by the Free
|
|
Packit |
5c3484 |
Software Foundation; either version 3 of the License, or (at your
|
|
Packit |
5c3484 |
option) any later version.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
or
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
* the GNU General Public License as published by the Free Software
|
|
Packit |
5c3484 |
Foundation; either version 2 of the License, or (at your option) any
|
|
Packit |
5c3484 |
later version.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
or both in parallel, as here.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The GNU MP Library is distributed in the hope that it will be useful, but
|
|
Packit |
5c3484 |
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
|
|
Packit |
5c3484 |
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
Packit |
5c3484 |
for more details.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
You should have received copies of the GNU General Public License and the
|
|
Packit |
5c3484 |
GNU Lesser General Public License along with the GNU MP Library. If not,
|
|
Packit |
5c3484 |
see https://www.gnu.org/licenses/. */
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#include "gmp.h"
|
|
Packit |
5c3484 |
#include "gmp-impl.h"
|
|
Packit |
5c3484 |
#include "longlong.h"
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* Set q to a/d, expecting d to be from a GCD and therefore usually small.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The distribution of GCDs of random numbers can be found in Knuth volume 2
|
|
Packit |
5c3484 |
section 4.5.2 theorem D.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
GCD chance
|
|
Packit |
5c3484 |
1 60.8%
|
|
Packit |
5c3484 |
2^k 20.2% (1<=k<32)
|
|
Packit |
5c3484 |
3*2^k 9.0% (1<=k<32)
|
|
Packit |
5c3484 |
other 10.1%
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Only the low limb is examined for optimizations, since GCDs bigger than
|
|
Packit |
5c3484 |
2^32 (or 2^64) will occur very infrequently.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Future: This could change to an mpn_divexact_gcd, possibly partly
|
|
Packit |
5c3484 |
inlined, if/when the relevant mpq functions change to an mpn based
|
|
Packit |
5c3484 |
implementation. */
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if GMP_NUMB_BITS % 2 == 0
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
mpz_divexact_by3 (mpz_ptr q, mpz_srcptr a)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_size_t size = SIZ(a);
|
|
Packit |
5c3484 |
mp_size_t abs_size = ABS(size);
|
|
Packit |
5c3484 |
mp_ptr qp;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
qp = MPZ_REALLOC (q, abs_size);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mpn_bdiv_dbm1 (qp, PTR(a), abs_size, GMP_NUMB_MASK / 3);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
abs_size -= (qp[abs_size-1] == 0);
|
|
Packit |
5c3484 |
SIZ(q) = (size>0 ? abs_size : -abs_size);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if GMP_NUMB_BITS % 4 == 0
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
mpz_divexact_by5 (mpz_ptr q, mpz_srcptr a)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_size_t size = SIZ(a);
|
|
Packit |
5c3484 |
mp_size_t abs_size = ABS(size);
|
|
Packit |
5c3484 |
mp_ptr qp;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
qp = MPZ_REALLOC (q, abs_size);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mpn_bdiv_dbm1 (qp, PTR(a), abs_size, GMP_NUMB_MASK / 5);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
abs_size -= (qp[abs_size-1] == 0);
|
|
Packit |
5c3484 |
SIZ(q) = (size>0 ? abs_size : -abs_size);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
mpz_divexact_limb (mpz_ptr q, mpz_srcptr a, mp_limb_t d)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_size_t size = SIZ(a);
|
|
Packit |
5c3484 |
mp_size_t abs_size = ABS(size);
|
|
Packit |
5c3484 |
mp_ptr qp;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
qp = MPZ_REALLOC (q, abs_size);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mpn_divexact_1 (qp, PTR(a), abs_size, d);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
abs_size -= (qp[abs_size-1] == 0);
|
|
Packit |
5c3484 |
SIZ(q) = (size>0 ? abs_size : -abs_size);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
void
|
|
Packit |
5c3484 |
mpz_divexact_gcd (mpz_ptr q, mpz_srcptr a, mpz_srcptr d)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
ASSERT (mpz_sgn (d) > 0);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (SIZ(a) == 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
SIZ(q) = 0;
|
|
Packit |
5c3484 |
return;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (SIZ(d) == 1)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_limb_t dl = PTR(d)[0];
|
|
Packit |
5c3484 |
int twos;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if ((dl & 1) == 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
count_trailing_zeros (twos, dl);
|
|
Packit |
5c3484 |
dl >>= twos;
|
|
Packit |
5c3484 |
mpz_tdiv_q_2exp (q, a, twos);
|
|
Packit |
5c3484 |
a = q;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (dl == 1)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (q != a)
|
|
Packit |
5c3484 |
mpz_set (q, a);
|
|
Packit |
5c3484 |
return;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
#if GMP_NUMB_BITS % 2 == 0
|
|
Packit |
5c3484 |
if (dl == 3)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mpz_divexact_by3 (q, a);
|
|
Packit |
5c3484 |
return;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if GMP_NUMB_BITS % 4 == 0
|
|
Packit |
5c3484 |
if (dl == 5)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mpz_divexact_by5 (q, a);
|
|
Packit |
5c3484 |
return;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mpz_divexact_limb (q, a, dl);
|
|
Packit |
5c3484 |
return;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mpz_divexact (q, a, d);
|
|
Packit |
5c3484 |
}
|