hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-05 09:21:38 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

improve speed (hopefully) of division with GMP

Browse source
This commit is contained in:
Fredrik Johansson 2014-12-09 13:58:45 +01:00
parent d72a4b1cd0
commit 210337d95b
1 changed files with 15 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

15
arf/div.c
View file

@ -26,6 +26,12 @@
#include "arf.h"
#include "mpn_extras.h"
#ifdef mpn_tdiv
#define USE_GMP_DIV_Q 1
#else
#define USE_GMP_DIV_Q 0
#endif
void
arf_div_special(arf_t z, const arf_t x, const arf_t y)
{
@ -89,6 +95,10 @@ arf_div(arf_ptr z, arf_srcptr x, arf_srcptr y, long prec, arf_rnd_t rnd)
tn = xn + sn;
zn = tn - yn + 1;
alloc = zn + (tn + 1);
/* need tn + 1 extra temporary limbs, which we store at the end of tptr */
#if USE_GMP_DIV_Q
alloc += tn + 1;
#endif
ARF_MUL_TMP_ALLOC(tmp, alloc)
@ -97,7 +107,12 @@ arf_div(arf_ptr z, arf_srcptr x, arf_srcptr y, long prec, arf_rnd_t rnd)
flint_mpn_zero(tptr, sn);
flint_mpn_copyi(tptr + sn, xptr, xn);
/* uses tn + 1 extra temporary limbs, tn limbs after tptr */
#if USE_GMP_DIV_Q
mpn_div_q(zptr, tptr, tn, yptr, yn, tptr + tn);
#else
mpn_tdiv_q(zptr, tptr, tn, yptr, yn);
#endif
if (zptr[zn - 1] == 0)
{