Salamin–Brent algorithm