From 8226061da437c8dc27557e65c28266b49149c34e Mon Sep 17 00:00:00 2001 From: PavelMikus Date: Tue, 15 Feb 2022 13:33:15 +0100 Subject: [PATCH] so far best version --- src/libslic3r/GCode/SeamPlacerNG.cpp | 104 ++++++++++----------------- src/libslic3r/GCode/SeamPlacerNG.hpp | 16 ++--- 2 files changed, 46 insertions(+), 74 deletions(-) diff --git a/src/libslic3r/GCode/SeamPlacerNG.cpp b/src/libslic3r/GCode/SeamPlacerNG.cpp index f87dbae8e..b6cba310d 100644 --- a/src/libslic3r/GCode/SeamPlacerNG.cpp +++ b/src/libslic3r/GCode/SeamPlacerNG.cpp @@ -15,26 +15,10 @@ #include "libslic3r/SVG.hpp" #include "libslic3r/Layer.hpp" -#include - namespace Slic3r { namespace SeamPlacerImpl { -void atomic_fetch_add_float(std::atomic &atomic, float increment) - { - float old_val; - float new_val; - - do - { - old_val = atomic.load(std::memory_order_relaxed); - new_val = old_val + increment; - } while (!atomic.compare_exchange_weak(old_val, new_val, - std::memory_order_release, - std::memory_order_relaxed)); -} - /// Coordinate frame class Frame { public: @@ -171,20 +155,6 @@ std::vector raycast_visibility(const AABBTreeIndirect::Tree<3, float> & BOOST_LOG_TRIVIAL(debug) << "SeamPlacer: raycast visibility for " << ray_count << " rays: end"; - its_write_obj(triangles, "triangles.obj"); - - Slic3r::CNumericLocalesSetter locales_setter; - FILE *fp = boost::nowide::fopen("hits.obj", "w"); - if (fp == nullptr) { - BOOST_LOG_TRIVIAL(error) - << "Couldn't open " << "hits.obj" << " for writing"; - } - - for (size_t i = 0; i < hit_points.size(); ++i) - fprintf(fp, "v %f %f %f \n", hit_points[i].m_position[0], hit_points[i].m_position[1], - hit_points[i].m_position[2]); - fclose(fp); - return hit_points; } @@ -495,24 +465,31 @@ struct DefaultSeamComparator { auto angle_score = [](float ccw_angle) { if (ccw_angle > 0) { - float normalized = (ccw_angle / float(PI)); - return normalized * normalized * normalized * 0.8; + float normalized = (ccw_angle / float(PI)) * 0.9; + return normalized; } else { - float normalized = (-ccw_angle / float(PI)); - return normalized * normalized * normalized * 1.0; + float normalized = (-ccw_angle / float(PI)) * 1.1; + return normalized; } }; + float angle_weight = 2.0; auto vis_score = [](float visibility) { - return 1.0 - visibility / SeamPlacer::expected_hits_per_area; + return (1.0 - visibility / SeamPlacer::expected_hits_per_area); }; + float vis_weight = 1.2; auto align_score = [](float nearby_seams) { - return nearby_seams / (0.25 * (sqrt(2) * SeamPlacer::seam_align_layer_dist)); + return nearby_seams / SeamPlacer::seam_align_layer_dist; }; + float align_weight = 1.0; - float score_a = angle_score(a.m_ccw_angle) + vis_score(a.m_visibility) + align_score(*a.m_nearby_seam_points); - float score_b = angle_score(b.m_ccw_angle) + vis_score(b.m_visibility) + align_score(*b.m_nearby_seam_points); + float score_a = angle_score(a.m_ccw_angle) * angle_weight + + vis_score(a.m_visibility) * vis_weight + + align_score(*a.m_nearby_seam_points) * align_weight; + float score_b = angle_score(b.m_ccw_angle) * angle_weight + + vis_score(b.m_visibility) * vis_weight + + align_score(*b.m_nearby_seam_points) * align_weight; if (score_a > score_b) return true; @@ -612,48 +589,42 @@ void SeamPlacer::distribute_seam_positions_for_alignment(const PrintObject *po) int other_layer_idx_top = std::min(layer_idx + seam_align_layer_dist, m_perimeter_points_per_object[po].size() - 1); + Vec3d last_point_position = seam_position; for (int other_layer_idx = layer_idx + 1; other_layer_idx <= other_layer_idx_top; ++other_layer_idx) { - - std::vector nearby_point_indexes = find_nearby_points( + Vec3d projected_position { last_point_position.x(), last_point_position.y(), po->get_layer(other_layer_idx)->slice_z}; + size_t closest_point_index = find_closest_point( *m_perimeter_points_trees_per_object[po][other_layer_idx], - seam_position, - seam_align_tolerable_dist * (other_layer_idx - layer_idx)); + projected_position); - if (nearby_point_indexes.empty()) { + SeamCandidate &point_ref = + m_perimeter_points_per_object[po][other_layer_idx][closest_point_index]; + if ((point_ref.m_position - projected_position).norm() + < SeamPlacer::seam_align_tolerable_dist) { + point_ref.m_nearby_seam_points->fetch_add(1, std::memory_order_relaxed); + last_point_position = point_ref.m_position; + } else { break; } - - for (size_t nearby_point_index : nearby_point_indexes) { - SeamCandidate &point_ref = - m_perimeter_points_per_object[po][other_layer_idx][nearby_point_index]; - float distance = (seam_position - point_ref.m_position).norm(); - atomic_fetch_add_float(*point_ref.m_nearby_seam_points, 1.0 / distance); - - } } - + last_point_position = seam_position; if (layer_idx > 0) { for (int other_layer_idx = layer_idx - 1; other_layer_idx >= other_layer_idx_bottom; --other_layer_idx) { - - std::vector nearby_point_indexes = find_nearby_points( + Vec3d projected_position { last_point_position.x(), last_point_position.y(), po->get_layer(other_layer_idx)->slice_z}; + size_t closest_point_index = find_closest_point( *m_perimeter_points_trees_per_object[po][other_layer_idx], - seam_position, - seam_align_tolerable_dist * (layer_idx - other_layer_idx)); + projected_position); - if (nearby_point_indexes.empty()) { + SeamCandidate &point_ref = + m_perimeter_points_per_object[po][other_layer_idx][closest_point_index]; + if ((point_ref.m_position - projected_position).norm() + < SeamPlacer::seam_align_tolerable_dist) { + point_ref.m_nearby_seam_points->fetch_add(1, std::memory_order_relaxed); + last_point_position = point_ref.m_position; + } else { break; } - - for (size_t nearby_point_index : nearby_point_indexes) { - SeamCandidate &point_ref = - m_perimeter_points_per_object[po][other_layer_idx][nearby_point_index]; - float distance = (seam_position - point_ref.m_position).norm(); - atomic_fetch_add_float(*point_ref.m_nearby_seam_points, 1.0 / distance); - - } - } } @@ -707,6 +678,7 @@ void SeamPlacer::init(const Print &print) { m_perimeter_points_per_object[po][layer_idx]; size_t current = 0; while (current < layer_perimeter_points.size()) { + //NOTE: pick seam point function also resets the m_nearby_seam_points count on all passed points pick_seam_point(layer_perimeter_points, current, current + layer_perimeter_points[current].m_polygon_index_reverse, DefaultSeamComparator { }); diff --git a/src/libslic3r/GCode/SeamPlacerNG.hpp b/src/libslic3r/GCode/SeamPlacerNG.hpp index 456b5d321..057640776 100644 --- a/src/libslic3r/GCode/SeamPlacerNG.hpp +++ b/src/libslic3r/GCode/SeamPlacerNG.hpp @@ -39,7 +39,7 @@ struct SeamCandidate { m_position(pos), m_visibility(0.0), m_overhang(0.0), m_polygon_index_reverse(polygon_index_reverse), m_seam_index( 0), m_ccw_angle( ccw_angle), m_type(type) { - m_nearby_seam_points = std::make_unique>(0.0); + m_nearby_seam_points = std::make_unique>(0); } Vec3d m_position; float m_visibility; @@ -47,7 +47,7 @@ struct SeamCandidate { size_t m_polygon_index_reverse; size_t m_seam_index; float m_ccw_angle; - std::unique_ptr> m_nearby_seam_points; + std::unique_ptr> m_nearby_seam_points; EnforcedBlockedSeamPoint m_type; }; @@ -81,14 +81,14 @@ class SeamPlacer { public: using SeamCandidatesTree = KDTreeIndirect<3, coordf_t, SeamPlacerImpl::SeamCandidateCoordinateFunctor>; - static constexpr double considered_hits_distance = 2.0; - static constexpr double expected_hits_per_area = 40.0; + static constexpr double considered_hits_distance = 4.0; + static constexpr double expected_hits_per_area = 250.0; static constexpr float cosine_hemisphere_sampling_power = 1.5; - static constexpr float polygon_angles_arm_distance = 1.0; + static constexpr float polygon_angles_arm_distance = 0.6; static constexpr float enforcer_blocker_sqr_distance_tolerance = 0.4; - static constexpr size_t seam_align_iterations = 5; - static constexpr size_t seam_align_layer_dist = 50; - static constexpr float seam_align_tolerable_dist = 1; + static constexpr size_t seam_align_iterations = 4; + static constexpr size_t seam_align_layer_dist = 30; + static constexpr float seam_align_tolerable_dist = 0.3; //perimeter points per object per layer idx, and their corresponding KD trees std::unordered_map>> m_perimeter_points_per_object; std::unordered_map>> m_perimeter_points_trees_per_object;