# HG changeset patch
# User Tuomo Valkonen <tuomov@iki.fi>
# Date 1733503092 18000
# Node ID 3d5c8ea1522ccdc7143a19f9684ff45ff591e175
# Parent  63318d1b4f00447852bea345aea4d20c5def02a5
assert tolerance

diff -r 63318d1b4f00 -r 3d5c8ea1522c src/cylinder.rs
--- 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))