openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-01-16 13:50:05 +01:00
Code Issues Releases Wiki Activity

optimization

Browse Source
This commit is contained in:
Justin Lin 2022-08-22 14:34:43 +08:00
parent 87ec19fb7a
commit 824168a006
1 changed files with 18 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/voronoi/_impl/_convex_intersection.scad
View File

@ -5,6 +5,23 @@ use <../../util/set/hashset_elems.scad>
include <../../__comm__/_pt2_hash.scad>
function __line_intersection2(line_pts1, line_pts2, epsilon = 0.0001) =
let(
a1 = line_pts1[0],
a2 = line_pts1[1],
b1 = line_pts2[0],
b2 = line_pts2[1],
a = a2 - a1,
b = b2 - b1,
c = cross(a, b)
)
abs(c) < epsilon ? [] : // they are parallel or conincident edges
let(
t = cross(b1 - a1, b) / c,
u = -cross(a1 - b1, a) / c
)
t >= 0 && t <= 1 && u >= 0 && u <= 1 ? a1 + a * t : [];
function _in_convex_r(convex_pts, pt, leng, i = 0) =
let(j = i + 1)
j == leng ||
@ -25,7 +42,7 @@ function _intersection_ps(closed_shape, line_pts, epsilon) =
let(
npts = [
for(i = [0:len(closed_shape) - 2])
let(p = lines_intersection(line_pts, [closed_shape[i], closed_shape[i + 1]], epsilon = epsilon))
let(p = __line_intersection2(line_pts, [closed_shape[i], closed_shape[i + 1]], epsilon = epsilon))
if(p != []) p
],
leng = len(npts)