openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-01-29 11:58:39 +01:00
Code Issues Releases Wiki Activity

optimization

Browse Source
This commit is contained in:
Justin Lin 2022-05-13 21:45:00 +08:00
parent db376e5290
commit d978128f8f
2 changed files with 19 additions and 32 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/triangle/_impl/_tri_delaunay_voronoi_impl.scad
View File

@ -1,13 +1,12 @@
use <../tri_delaunay.scad>;
use <../../util/map/hashmap_get.scad>;
function indicesOfCell(iTris, triIndices) =
_indicesOfCell(iTris, triIndices, len(iTris), [], iTris[0][0]);
function _indicesOfCell(iTris, triIndices, leng, indices, vi, i = 0) =
i == leng ? indices :
let(
t = iTris[search([vi], iTris)[0]],
nIndices = concat(indices, hashmap_get(triIndices, t))
)
_indicesOfCell(iTris, triIndices, leng, nIndices, t[1], i + 1);
let(leng = len(iTris))
[
for(i = 0, vi = iTris[0][0], t = iTris[search([vi], iTris)[0]];
i < leng;
i = i + 1,
vi = t[1],
t = i < leng ? iTris[search([vi], iTris)[0]] : undef)
hashmap_get(triIndices, t)
];

32
src/triangle/tri_delaunay_voronoi.scad
View File

@ -24,36 +24,24 @@ function tri_delaunay_voronoi(d) =
tris = hashmap_keys(delaunay_triangles(d)),
// circumcircle centers
vertices = [for(t = tris) hashmap_get(circles, t, hash = _indices_hash)[0]],
i_range = [0:len(tris) - 1],
i_rts = [
for(i = i_range)
let(
tris_i = tris[i],
a = tris_i[0],
b = tris_i[1],
c = tris_i[2],
rt1 = [b, c, a],
rt2 = [c, a, b],
rt3 = [a, b, c]
)
each [[a, rt1], [b, rt2], [c, rt3]]
for(tri = tris)
let(a = tri[0], b = tri[1], c = tri[2])
each [[a, [b, c, a]], [b, [c, a, b]], [c, tri]]
],
connectedTris = [
for(i = [0:coords_leng - 1])
[for(i_rt = i_rts) if(i_rt[0] == i) i_rt[1]]
[for(j = search(i, i_rts, num_returns_per_match=0)) i_rts[j][1]]
],
triIndices = hashmap([
for(i = i_range)
for(i = [0:len(tris) - 1])
let(
tris_i = tris[i],
a = tris_i[0],
b = tris_i[1],
c = tris_i[2],
rt1 = [b, c, a],
rt2 = [c, a, b],
rt3 = [a, b, c]
tri = tris[i],
a = tri[0],
b = tri[1],
c = tri[2]
)
each [[rt1, i], [rt2, i], [rt3, i]]
each [[[b, c, a], i], [[c, a, b], i], [tri, i]]
]),
cells = [
for(i = [4:coords_leng - 1])