Re: Anomaly in mpn_sqrtrem and mpn_rottrem

Il Sab, 13 Giugno 2015 12:27 am, bodr...@mail.dm.unipi.it ha scritto: > I wrote a simple patch (it touches very few lines) that allows skippingthe final I forgot the attachment... -- http://bodrato.it/ gmp.diff Description: plain/text ___ gmp-devel ma

Re: Anomaly in mpn_sqrtrem and mpn_rottrem

Ciao, Il Ven, 12 Giugno 2015 9:04 am, Torbjörn Granlund ha scritto: > We might want to look into the plain division-free sqrt(A) = A*sqrt(1/A) > approach before implementing a tricky division sqrt(A). We can try improving the current implementation, before implementing any other algorithm :-) I

Re: Anomaly in mpn_sqrtrem and mpn_rottrem

ni...@lysator.liu.se (Niels Möller) writes: Hmm. Or maybe this is stupid. I could stop insisting on using a full size inverse (so that A / x or E / x can be computed as a *single* mulhi), and instead work with a half-size inverse, so that the quotient is computed in two steps. Then inverse