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

improved performance

Browse Source
This commit is contained in:
Justin Lin
2019-05-04 10:04:50 +08:00
parent d09926285a
commit 9c57a539b2
1 changed files with 51 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

59
src/along_with.scad
View File

@@ -11,6 +11,11 @@
include <__private__/__angy_angz.scad>; include <__private__/__angy_angz.scad>;
include <__private__/__is_vector.scad>; include <__private__/__is_vector.scad>;
include <__private__/__to3d.scad>; include <__private__/__to3d.scad>;
include <__private__/__m_multiply.scad>;
// Becuase of improving the performance, this module requires m_rotation.scad which doesn't require in dotSCAD 1.0.
// For backward compatibility, I directly include m_rotation here.
include <m_rotation.scad>;
module along_with(points, angles, twist = 0, scale = 1.0) { module along_with(points, angles, twist = 0, scale = 1.0) {
leng_points = len(points); leng_points = len(points);
@@ -27,7 +32,44 @@ module along_with(points, angles, twist = 0, scale = 1.0) {
(scale[1] - 1) / leng_points_minus_one, (scale[1] - 1) / leng_points_minus_one,
scale[2] == undef ? 0 : (scale[2] - 1) / leng_points_minus_one scale[2] == undef ? 0 : (scale[2] - 1) / leng_points_minus_one
] : scale_step(); ] : scale_step();
// get rotation matrice for sections
function local_ang_vects(j) =
j == 0 ? [] : local_ang_vects_sub(j);
function local_ang_vects_sub(j) =
let(
vt0 = points[j] - points[j - 1],
vt1 = points[j + 1] - points[j],
a = acos((vt0 * vt1) / (norm(vt0) * norm(vt1))),
v = cross(vt0, vt1)
)
concat([[a, v]], local_ang_vects(j - 1));
function cumulated_rot_matrice(i, rot_matrice) =
let(
leng_rot_matrice = len(rot_matrice),
leng_rot_matrice_minus_one = leng_rot_matrice - 1,
leng_rot_matrice_minus_two = leng_rot_matrice - 2
)
i == leng_rot_matrice - 2 ?
[
rot_matrice[leng_rot_matrice_minus_one],
__m_multiply(rot_matrice[leng_rot_matrice_minus_two], rot_matrice[leng_rot_matrice_minus_one])
]
: cumulated_rot_matrice_sub(i, rot_matrice);
function cumulated_rot_matrice_sub(i, rot_matrice) =
let(
matrice = cumulated_rot_matrice(i + 1, rot_matrice),
curr_matrix = rot_matrice[i],
prev_matrix = matrice[len(matrice) - 1]
)
concat(matrice, [__m_multiply(curr_matrix, prev_matrix)]);
// align modules
module align_with_pts_angles(i) { module align_with_pts_angles(i) {
translate(points[i]) translate(points[i])
rotate(angles[i]) rotate(angles[i])
@@ -45,17 +87,14 @@ module along_with(points, angles, twist = 0, scale = 1.0) {
children(0); children(0);
} }
module align_with_pts_local_rotate(j, init_a, init_s) { module align_with_pts_local_rotate(j, init_a, init_s, cumu_rot_matrice) {
if(j == 0) { // first child if(j == 0) { // first child
align_with_pts_init(init_a, init_s) align_with_pts_init(init_a, init_s)
children(0); children(0);
} }
else { else {
vt0 = points[j] - points[j - 1]; multmatrix(cumu_rot_matrice[j - 1])
vt1 = points[j + 1] - points[j]; align_with_pts_init(init_a, init_s)
a = acos((vt0 * vt1) / (norm(vt0) * norm(vt1)));
rotate(a, cross(vt0, vt1))
align_with_pts_local_rotate(j - 1, init_a, init_s)
children(0); children(0);
} }
} }
@@ -72,8 +111,13 @@ module along_with(points, angles, twist = 0, scale = 1.0) {
} }
} }
else { else {
cumu_rot_matrice = cumulated_rot_matrice(0, [
for(ang_vect = local_ang_vects(leng_points - 2))
m_rotation(ang_vect[0], ang_vect[1])
]);
translate(points[0]) translate(points[0])
align_with_pts_local_rotate(0, 0, [1, 1, 1]) align_with_pts_local_rotate(0, 0, [1, 1, 1], cumu_rot_matrice)
children(0); children(0);
if($children == 1) { if($children == 1) {
@@ -90,5 +134,4 @@ module along_with(points, angles, twist = 0, scale = 1.0) {
} }
} }
} }
} }