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

Disambiguation of saddle points

Browse Source
This commit is contained in:
Justin Lin
2020-03-21 08:32:49 +08:00
parent 7d3443cfc6
commit 5cd8fb04fe
1 changed files with 24 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
src/experimental/_impl/_marching_squares_impl.scad
View File

@@ -39,10 +39,18 @@ function case4_contours(cell_pts, sigma) = [
[interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)] [interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)]
]; ];
function case5_contours(cell_pts, sigma) = [ function case5_contours(cell_pts, sigma) =
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[1], cell_pts[2], sigma)], let(mdpz = ((cell_pts[0] + cell_pts[1]) / 2)[2])
[interpolated_pt(cell_pts[0], cell_pts[3], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)] mdpz >= sigma ?
]; [
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[1], cell_pts[2], sigma)],
[interpolated_pt(cell_pts[0], cell_pts[3], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)]
]
:
[
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[0], cell_pts[3], sigma)],
[interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)]
];
function case6_contours(cell_pts, sigma) = [ function case6_contours(cell_pts, sigma) = [
[interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[0], cell_pts[3], sigma)] [interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[0], cell_pts[3], sigma)]
@@ -56,10 +64,18 @@ function case8_contours(cell_pts, sigma) = case7_contours(cell_pts, sigma);
function case9_contours(cell_pts, sigma) = case6_contours(cell_pts, sigma); function case9_contours(cell_pts, sigma) = case6_contours(cell_pts, sigma);
function case10_contours(cell_pts, sigma) = [ function case10_contours(cell_pts, sigma) =
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[0], cell_pts[3], sigma)], let(mdpz = ((cell_pts[0] + cell_pts[1]) / 2)[2])
[interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)] mdpz >= sigma ?
]; [
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[0], cell_pts[3], sigma)],
[interpolated_pt(cell_pts[1], cell_pts[2], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)]
]
:
[
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[1], cell_pts[2], sigma)],
[interpolated_pt(cell_pts[0], cell_pts[3], sigma), interpolated_pt(cell_pts[2], cell_pts[3], sigma)]
];
function case11_contours(cell_pts, sigma) = case4_contours(cell_pts, sigma); function case11_contours(cell_pts, sigma) = case4_contours(cell_pts, sigma);