openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-08-14 02:34:12 +02:00
Code Issues Releases Wiki Activity

remove density param

Browse Source
This commit is contained in:
Justin Lin
2022-06-03 10:37:45 +08:00
parent e0b32b96b7
commit 8371943c25
4 changed files with 33 additions and 46 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
examples/packing_circles.scad
View File

@@ -2,7 +2,6 @@ use <util/rand.scad>;
use <experimental/circle_packing.scad>;
size = [200, 100];
density = 2;
min_radius = 1;
point_numbers = 100;
@@ -11,7 +10,7 @@ points = [
[rand(0, size.x), rand(0, size.y)]
];
circles = circle_packing(points, density, min_radius);
circles = circle_packing(points, min_radius);
mr = max([for(c = circles) c[1]]);
translate([0, 0, mr])
for(c = circles) {

12
src/experimental/circle_packing.scad
View File

@@ -2,19 +2,23 @@ use <experimental/tri_circle_packing.scad>;
use <triangle/tri_delaunay.scad>;
function circle_packing(points, density = 1, min_r = 1) =
function circle_packing(points, min_r = 1) =
[
for(t = tri_delaunay(points, ret = "TRI_SHAPES"))
each tri_circle_packing(t, density, min_r)
each tri_circle_packing(t, min_r)
];
$fn = 24;
density = 3;
min_r = 1;
points = [for(i = [0:100]) rands(-100, 100, 2)];
for(c = circle_packing(points, density, min_r)) {
for(c = circle_packing(points, min_r)) {
translate(c[0])
sphere(c[1]);
}
for(t = tri_delaunay(points, ret = "TRI_SHAPES")) {
offset(-.1)
polygon(t);
}

2
src/experimental/circle_packing3.scad
View File

@@ -13,7 +13,7 @@ function circle_packing3(points, density = 1, min_r = 1) =
function circle_packing_triangle3(t, density, min_r) =
[
for(st = tri_subdivide(t, density))
each tri_circle_packing(st, density, min_r)
each tri_circle_packing(st, min_r)
];
$fn = 24;

62
src/experimental/tri_circle_packing.scad
View File

@@ -1,7 +1,7 @@
use <triangle/tri_incenter.scad>;
use <ptf/ptf_rotate.scad>;
function r(sin, leng_cv, pre_R) = sin * (leng_cv - pre_R) / (1 + sin);
function r2(sinv, leng_cv, pre_R) = sinv * (leng_cv - pre_R) / (1 + sinv);
// the 3rd circle
function c3_r(r1, r2) = (r1 * r2) / (r1 + r2 + 2 * sqrt(r1 * r2));
@@ -13,7 +13,7 @@ function c3_a(r1, r2, r3) =
)
acos((b ^ 2 + c ^ 2 - a ^ 2) / (2 * b * c));
function tri_circle_packing(t, density, min_r) =
function tri_circle_packing(t, min_r) =
let(
center = tri_incenter(t),
s1 = t[1] - t[0],
@@ -35,45 +35,29 @@ function tri_circle_packing(t, density, min_r) =
leng_cc = norm(cc),
unit_cc = cc / leng_cc,
sinc = R / leng_cc,
_small_circles = function(density, pre_leng_a = R, pre_leng_b = R, pre_leng_c = R, pre_ra = R, pre_rb = R, pre_rc = R)
density <= 0 ? [] :
pack_a = function(sinv, leng_cv, unit_cv, pre_leng = R, pre_r = R)
let(
ra = r(sina, leng_ca, pre_leng_a),
Ra = pre_leng_a + ra,
ct_a = center - unit_ca * Ra,
r2 = r2(sinv, leng_cv, pre_leng),
leng = pre_leng + r2,
r2_ct = center - unit_cv * leng,
rb = r(sinb, leng_cb, pre_leng_b),
Rb = pre_leng_b + rb,
ct_b = center - unit_cb * Rb,
rc = r(sinc, leng_cc, pre_leng_c),
Rc = pre_leng_c + rc,
ct_c = center - unit_cc * Rc,
r3a = c3_r(pre_ra, ra),
alpha3a = c3_a(pre_ra, ra, r3a),
vta = unit_ca * (r3a + ra),
ct1a = ct_a + ptf_rotate(vta, alpha3a),
ct2a = ct_a + ptf_rotate(vta, -alpha3a),
r3b = c3_r(pre_rb, rb),
alpha3b = c3_a(pre_rb, rb, r3b),
vtb = unit_cb * (r3b + rb),
ct1b = ct_b + ptf_rotate(vtb, alpha3b),
ct2b = ct_b + ptf_rotate(vtb, -alpha3b),
r3c = c3_r(pre_rc, rc),
alpha3c = c3_a(pre_rc, rc, r3c),
vtc = unit_cc * (r3c + rc),
ct1c = ct_c + ptf_rotate(vtc, alpha3c),
ct2c = ct_c + ptf_rotate(vtc, -alpha3c)
r3 = c3_r(pre_r, r2),
a = c3_a(pre_r, r2, r3),
vta = unit_cv * (r3 + r2),
r3_ct1 = r2_ct + ptf_rotate(vta, a),
r3_ct2 = r2_ct + ptf_rotate(vta, -a)
)
r2 > min_r ?
concat(
[[r2_ct, r2], if(r3 > min_r) each [[r3_ct1, r3], [r3_ct2, r3]]],
pack_a(sinv, leng_cv, unit_cv, r2 + leng, r2)
)
: []
)
[
if(ra > min_r) each [[ct_a, ra], if(r3a > min_r) each [[ct1a, r3a], [ct2a, r3a]]],
if(rb > min_r) each [[ct_b, rb], if(r3b > min_r) each [[ct1b, r3b], [ct2b, r3b]]],
if(rc > min_r) each [[ct_c, rc], if(r3c > min_r) each [[ct1c, r3c], [ct2c, r3c]]],
each _small_circles(density - 1, ra + Ra, rb + Rb, rc + Rc, ra, rb, rc)
]
)
[[center, R], each _small_circles(density - 1)];
[center, R],
each pack_a(sina, leng_ca, unit_ca),
each pack_a(sinb, leng_cb, unit_cb),
each pack_a(sinc, leng_cc, unit_cc),
];