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

fix d error

Browse Source
This commit is contained in:
Justin Lin
2021-06-01 17:18:03 +08:00
parent 009230c8f8
commit 32c407ff0b
1 changed files with 9 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
src/_impl/_bspline_curve_impl.scad
View File

@@ -21,29 +21,29 @@ function _bspline_curve_s(s, end, t, knots) =
function _bspline_curve_alpha(i, l, t, degree, knots) = function _bspline_curve_alpha(i, l, t, degree, knots) =
(t - knots[i]) / (knots[i + degree + 1 - l] - knots[i]); (t - knots[i]) / (knots[i + degree + 1 - l] - knots[i]);
function _bspline_curve_nvi(v, i, l, t, degree, knots) = function _bspline_curve_nvi(v, i, l, t, degree, knots, d) =
let( let(
alpha = _bspline_curve_alpha(i, l, t, degree, knots) alpha = _bspline_curve_alpha(i, l, t, degree, knots)
) )
[[for(j = 0; j< degree + 1; j = j + 1) ((1 - alpha) * v[i - 1][j] + alpha * v[i][j])]]; [[for(j = 0; j< d + 1; j = j + 1) ((1 - alpha) * v[i - 1][j] + alpha * v[i][j])]];
function _bspline_curve_nvl(v, l, s, t, degree, knots, i) = function _bspline_curve_nvl(v, l, s, t, degree, knots, d, i) =
i == (s - degree - 1 + l) ? v : i == (s - degree - 1 + l) ? v :
let( let(
leng_v = len(v), leng_v = len(v),
nvi = _bspline_curve_nvi(v, i, l, t, degree, knots), nvi = _bspline_curve_nvi(v, i, l, t, degree, knots, d),
nv = concat( nv = concat(
[for(j = 0; j < i; j = j + 1) v[j]], [for(j = 0; j < i; j = j + 1) v[j]],
nvi, nvi,
[for(j = i + 1; j < leng_v; j = j + 1) v[j]] [for(j = i + 1; j < leng_v; j = j + 1) v[j]]
) )
) )
_bspline_curve_nvl(nv, l, s, t, degree, knots, i - 1); _bspline_curve_nvl(nv, l, s, t, degree, knots, d, i - 1);
function _bspline_curve_v(v, s, t, degree, knots, l = 1) = function _bspline_curve_v(v, s, t, degree, knots, d, l = 1) =
l > degree + 1 ? v : l > degree + 1 ? v :
let(nv = _bspline_curve_nvl(v, l, s, t, degree, knots, s)) let(nv = _bspline_curve_nvl(v, l, s, t, degree, knots, d, s))
_bspline_curve_v(nv, s, t, degree, knots, l + 1); _bspline_curve_v(nv, s, t, degree, knots, d, l + 1);
function _bspline_curve_interpolate(t, degree, points, knots, weights) = function _bspline_curve_interpolate(t, degree, points, knots, weights) =
let( let(
@@ -58,7 +58,7 @@ function _bspline_curve_interpolate(t, degree, points, knots, weights) =
], ],
ts = _bspline_curve_ts(t, degree, kts), ts = _bspline_curve_ts(t, degree, kts),
s = ts[1], s = ts[1],
nv = _bspline_curve_v(v, s, ts[0], degree, kts) nv = _bspline_curve_v(v, s, ts[0], degree, kts, d)
) )
[for(i = 0; i < d; i = i + 1) nv[s][i] / nv[s][d]]; [for(i = 0; i < d; i = i + 1) nv[s][i] / nv[s][d]];