pub fn is_boundary_zone(norm: f64, rho: f64, boundary_fraction: f64) -> bool
Check if a residual norm is in the boundary zone (> boundary_fraction * ρ)