Hi, > Can this be faster than the current code? I don't know. I have to say that in MPFR I tried to implement Newton's method for 1/sqrt(x), which is theoretically faster than mpn_sqrtrem, by implementing Algorithm ApproximateRecSquareRoot from Modern Computer Arithmetic, but I failed miserably. Good luck! Paul