Re: mpn_sqrtrem{1,2} - patch for pure C implem

Dear Marco, > Finally we can totally skip the first approximation and directly compute > > i' = T1[h(a)] + l(a) * T2[h(a)] > > using a piece-wise linear function. yes, we noticed with Niels that we can replace most of the arithmetic in the first Newton iteration by another table lookup.

Re: mpn_sqrtrem{1,2} - patch for pure C implem

"Marco Bodrato" writes: > The constant (c) I added to the first iteration should be incorporated in > T1, and of course it is not necessary to really use a constant. A > different value can be used for each entry in the table... In Newton iteration, one usually