return pow 2 of radius

2025-08-16 03:34:42 +02:00 · 2021-04-02 16:08:51 +08:00
parent a9e1add8bb
commit ea8a5d87ed
1 changed files with 6 additions and 3 deletions
--- a/src/experimental/tri_circumcircle.scad
+++ b/src/experimental/tri_circumcircle.scad
@@ -9,10 +9,13 @@ function tri_circumcircle(points) =
      d1 = (p2 + p1) / 2 * v1,
      det = -cross(v0 , v1)
   )
-   det == 0? [] : 
+   det == 0 ? undef : 
             let(
                 x = (d1 * v0[1] - d0 * v1[1]) / det,
                 y = (d0 * v1[0] - d1 * v0[0]) / det,
-                 r = norm(p0 - [x,y])
+                 center = [x, y],
+                 v = p0 - center,
+                 r = norm(v),
+                 rr = pow(v[0], 2) + pow(v[1], 2)
             )
-             [[x,y], r];
+             [center, r, rr];