Fri, 06 Dec 2024 11:38:12 -0500
assert tolerance
src/cylinder.rs | file | annotate | diff | comparison | revisions |
--- a/src/cylinder.rs Thu Dec 05 16:50:08 2024 -0500 +++ b/src/cylinder.rs Fri Dec 06 11:38:12 2024 -0500 @@ -104,7 +104,7 @@ let nt2 = t.norm2_squared(); let r2 = r * r; let d = p.dot(t); - assert!(np2 <= r2); + assert!(np2 <= r2 + f64::EPSILON, "‖{p}‖ = {} > {r}", np2.sqrt()); let s = (-d + (d*d + nt2 * (r2 - np2)).sqrt()) / nt2; if s < 1.0 { (Self::from_cartesian(p + s * t), Some((1.0-s) * t))