Salamin-Brent algorithm