Wydajne algorytmy pierwiastkowania liczb binarnych.

Cześć,
poszukuję najbardziej wydajnych algorytmów pierwiastkowania liczb binarnych rozszerzonego zakresu (powiedzmy około 1000 - bitowych liczb). Zależy mi żeby wyselekcjonować z nich 2 najlepsze i zaimplementować je w assemblerze.
Czy ktoś mógłby podpowiedzieć z doświadczenia jakie algorytmy są najlepsze i najwydajniejsze?

Pozdrawiam i z góry dziękuję za każdą propozycję.

http://en.wikipedia.org/wiki/Methods_of_computing_square_roots

Chyba wystarczy zwyczajna metoda Newtona.
Zaczynasz od pierwiastka z precyzją double, czyli 53 bity, albo 64 dla long double, a Newton potem to podwaja w każdym kroku.

1000/53 < 20,
zatem dla uzyskania 1000 bitów precyzji 5 kroków Newtona chyba wystarczy.

53 -> 106 -> 212 -> 424 -> 848 -> 1696 bitów precyzji dla 5-ciu iteracji.