Simplification of "reciprocal word 3by2" (aka invert_pi1)

Paul Zimmermann Paul.Zimmermann at inria.fr
Mon Jun 22 07:01:27 UTC 2020


       Dear Pawel,

> are you able also to generate a case for t0 == d0?

after thought, this is not possible. Indeed, if we write d=d1*beta+d0, we have
(p,t0) = (beta+v)*d % beta^2 at line 14, thus it would mean that
(beta+v)*d = beta^3 + d, thus d would divide beta^3, which implies that d is
a power of two. The only possible case is d1=beta/2 and d0=0. But when d0=0
the second part of the algorithm at lines 10-15 does no correction.

Paul


More information about the gmp-discuss mailing list