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

refactor: use for

Browse Source
This commit is contained in:
Justin Lin
2022-05-07 16:24:00 +08:00
parent 62ef13a8bd
commit ac262336fb
1 changed files with 13 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

29
src/_impl/_archimedean_spiral_impl.scad
View File

@@ -1,26 +1,23 @@
use <../util/radians.scad>;
use <../util/degrees.scad>;
function _radian_step(b, theta, dist) =
let(r_square = (b * theta) ^ 2)
radians(acos(1 - dist ^ 2 / (2 * r_square)));
function _radian_step(b, radian, pow2_dist) =
radians(acos(1 - pow2_dist / (2 * (b * radian) ^ 2)));
function _find_radians(b, point_distance, rads, n, count = 1) =
let(pre_rads = rads[count - 1])
count == n ? rads : (
_find_radians(
b,
point_distance,
[each rads, pre_rads + _radian_step(b, pre_rads, point_distance)],
n,
count + 1
)
);
function _find_radians(b, pow2_dist, init_radian, n) =
[
for(
count = 0, radian = init_radian;
count < n;
count = count + 1, radian = radian + _radian_step(b, radian, pow2_dist)
)
radian
];
function _archimedean_spiral_impl(arm_distance, init_angle, point_distance, num_of_points, rt_dir) =
let(b = arm_distance / (2 * PI), init_radian = radians(init_angle), sgn = rt_dir == "CT_CLK" ? 1 : -1)
[
for(theta = _find_radians(b, point_distance, [init_radian], num_of_points))
let(r = b * theta, a = degrees(sgn * theta))
for(radian = _find_radians(b, point_distance ^ 2, init_radian, num_of_points))
let(r = b * radian, a = degrees(sgn * radian))
[r * [cos(a), sin(a)], a]
];