div_qr_1 interface

Torbjorn Granlund tg at gmplib.org
Thu Oct 24 23:49:41 CEST 2013


I pushed initial C versions of these functions:

    mpn_div_qr_1n_pi2
    mpn_div_qr_1u_pi2

I have had these for a long time, judging from the file time stamps.

These accept n-limb dividends in a single consecutive operand and
generate n-limb quotients also in a consecutive operand.  I now realise
that you use a n+1 limb dividend passed as separate operands for the low
n limbs and the high limb.  I had misunderstood your intended
interface...

The new pi2 functions follow your suggestion to stuff the dividend
together with the various precomputations in a single operand.  Here it
is a structure, perhaps we could use a simple array instead (like cps).

-- 
Torbjörn


More information about the gmp-devel mailing list