Function sqrt_newton_

src/fl/math/math.cpp.hpp:151–163

// space and rotational math.
template <typename F>
inline F sqrt_newton_(F value) FL_NOEXCEPT {
    if (value <= F(0)) return F(0);
    // Initial estimate: x_0 = value / 2. Pure-arithmetic, no bit-cast needed.
    // Newton converges quadratically so 5-6 iterations is plenty even from a
    // crude start.
    F x = value;
    if (x > F(1)) x = F(1) + (x - F(1)) * F(0.5);  // bias toward 1 for fast convergence
    for (int i = 0; i < 6; ++i) {
        if (x == F(0)) break;
        x = F(0.5) * (x + value / x);
    }
    return x;
}

// --- Polynomial sin/cos with range reduction ------------------------------
//