divappr tests

Niels Möller nisse at lysator.liu.se
Mon Dec 14 22:32:07 CET 2009

bodrato at mail.dm.unipi.it writes:

> I tried your patch, but "make tune" aborted with:
> #define DC_DIVAPPR_Q_THRESHOLD        sbpi1_divappr_q.c:181: GNU MP
> assertion failed: np[1] == n1
> *** Signal 6

It's the patch of the final case which doesn't satisfy the assumptions
of that assert (by storing a result directly in np[1] instead of going
via the n1 variable). So I'll delete that assert and check in.

I also have some concerns over the divappr interface (basically, to
produce qn quotient words it should need only qn+1 limbs of N and at
most qn+1 limbs of d, I'll have to get the details down in a separate


Niels Möller. PGP-encrypted email is preferred. Keyid C0B98E26.
Internet email is subject to wholesale government surveillance.

More information about the gmp-devel mailing list