openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-08-01 20:40:28 +02:00
Code Issues Releases Wiki Activity

optimization

Browse Source
This commit is contained in:
Justin Lin
2022-05-05 12:05:17 +08:00
parent e714fb0bc3
commit 65a19f5f42
1 changed files with 15 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

28
src/polyhedra/_impl/_geom_platonic_polyhedra.scad
View File

@@ -1,22 +1,24 @@
function _tri_subdivide(points, detail) = function _tri_subdivide(points, detail) =
let( let(
rows = detail + 1, rows = detail + 1,
dr = (points[2] - points[0]) / rows, p0 = points[0],
dc = (points[1] - points[0]) / rows, basis = [points[2] - p0, points[1] - p0] / rows,
pts = [ pts = [
for(ri = [0:rows], ci = [0:rows - ri]) for(ri = [0:rows], ci = [0:rows - ri])
points[0] + ci * dc + ri * dr p0 + [ri, ci] * basis
], ],
ri_base = [for(ri = 0, acc = 0; ri <= rows; ri = ri + 1, acc = acc + rows - ri + 2) acc], ri_base = [for(ri = 0, row_2 = rows + 2; ri <= rows; ri = ri + 1) ri * row_2 - ri * (ri + 1) * 0.5],
idx = function(ci, ri) ci + ri_base[ri],
faces = [ faces = [
for(ri = [0:rows - 1]) for(ri = [0:rows - 1], base_i = ri_base[ri], base_i1 = ri_base[ri + 1])
let(cols = rows - ri - 1) let(cols = rows - ri - 1)
for(ci = [0:cols]) for(ci = [0:cols])
each [ if(ci != cols)
[idx(ci + 1, ri), idx(ci, ri + 1), idx(ci, ri)], each [
if(ci != cols) [idx(ci + 1, ri + 1), idx(ci, ri + 1), idx(ci + 1, ri)] [ci + 1 + base_i, ci + base_i1, ci + base_i],
] [ci + 1 + base_i1, ci + base_i1, ci + 1 + base_i]
]
else
[ci + 1 + base_i, ci + base_i1, ci + base_i]
] ]
) )
[pts, faces]; [pts, faces];
@@ -26,18 +28,18 @@ function _subdivide_project(points, faces, radius, detail) =
let( let(
subdivided_all = [ subdivided_all = [
for(face = faces) for(face = faces)
_tri_subdivide([for(i = face) points[i]], detail) _tri_subdivide([for(i = face) points[i]], detail)
], ],
flatten_points = [ flatten_points = [
for(pts_faces = subdivided_all, p = pts_faces[0]) for(pts_faces = subdivided_all, p = pts_faces[0])
p / norm(p) * radius p / norm(p) * radius
], ],
pts_number_per_tri = len(subdivided_all[0][0]), pts_number_per_tri = len(subdivided_all[0][0]),
flatten_faces = [ flatten_faces = [
for(i = [0:len(subdivided_all) - 1]) for(i = [0:len(subdivided_all) - 1])
let(faces = subdivided_all[i][1], off = [i, i, i] * pts_number_per_tri) let(faces = subdivided_all[i][1], off = [i, i, i] * pts_number_per_tri)
for(face = faces) for(face = faces)
face + off face + off
] ]
) )
[flatten_points, flatten_faces]; [flatten_points, flatten_faces];