Function sqrt_binary_search

src/math/mysqrt.rs:1–20

(x: f64)

pub fn sqrt_binary_search(x: f64) -> f64 {
    assert!(x >= 0.0);
    let mut low = 0.0;
    let mut up = x;
    let mut last = 0.0;
    loop {
        let mid = (low + up) / 2.0;
        if (mid - last).abs() <= f64::EPSILON {
            return mid;
        }

        if mid * mid > x {
            up = mid;
        } else {
            low = mid;
        }

        last = mid;
    }
}

pub fn sqrt_newton(x: f64) -> f64 {
    assert!(x >= 0.0);