ENH: add filament group strategy

1.When capacity is greater than the num of filaments, always choose the
map result that can be accommodated
2.In BestFit strategy,always try to fill up the existing capacity
3.In BestCost strategy, just try the group with fewest flush

jira:NEW

Signed-off-by: xun.zhang <xun.zhang@bambulab.com>
Change-Id: Ifd6d64b77774039e57ffff26cf2243a4d3f89054
This commit is contained in:
xun.zhang 2024-07-12 18:04:01 +08:00 committed by lane.wei
parent 8232a716e6
commit cddf8cae27
3 changed files with 82 additions and 48 deletions

View File

@ -3,7 +3,7 @@
namespace Slic3r namespace Slic3r
{ {
int FilamentGroup::calc_filament_group(const std::vector<std::vector<unsigned int>>& layer_filaments) int FilamentGroup::calc_filament_group(const std::vector<std::vector<unsigned int>>& layer_filaments,const FGStrategy& g_strategy)
{ {
std::set<unsigned int>used_filaments; std::set<unsigned int>used_filaments;
for (const auto& lf : layer_filaments) for (const auto& lf : layer_filaments)
@ -17,15 +17,15 @@ namespace Slic3r
if (m_filament_num <= 1) if (m_filament_num <= 1)
return 0; return 0;
if (m_filament_num < 10) if (m_filament_num < 10)
return calc_filament_group_by_enum(layer_filaments); return calc_filament_group_by_enum(layer_filaments,g_strategy);
else else
return calc_filament_group_by_pam(layer_filaments,300); return calc_filament_group_by_pam(layer_filaments,g_strategy,300);
} }
int FilamentGroup::calc_filament_group_by_enum(const std::vector<std::vector<unsigned int>>& layer_filaments) int FilamentGroup::calc_filament_group_by_enum(const std::vector<std::vector<unsigned int>>& layer_filaments,const FGStrategy& g_strategy)
{ {
auto bit_count_one = [](int n) auto bit_count_one = [](uint64_t n)
{ {
int count = 0; int count = 0;
while (n != 0) while (n != 0)
@ -36,14 +36,26 @@ namespace Slic3r
return count; return count;
}; };
bool have_enough_size = (m_filament_num <= (m_max_group_size[0] + m_max_group_size[1]));
uint64_t max_group_num = static_cast<uint64_t>(1 << m_filament_num); uint64_t max_group_num = static_cast<uint64_t>(1 << m_filament_num);
int best_cost = std::numeric_limits<int>::max(); int best_cost = std::numeric_limits<int>::max();
std::vector<int>best_label; std::vector<int>best_label;
for (uint64_t i = 0; i < max_group_num; ++i) { for (uint64_t i = 0; i < max_group_num; ++i) {
int num_to_group_1 = bit_count_one(i); int num_to_group_1 = bit_count_one(i);
if (num_to_group_1 > m_max_group_size[1] || (m_filament_num - num_to_group_1) > m_max_group_size[0]) int num_to_group_0 = m_filament_num - num_to_group_1;
bool should_accept = false;
if (have_enough_size)
should_accept = (num_to_group_0 <= m_max_group_size[0] && num_to_group_1 <= m_max_group_size[1]);
else if (g_strategy == FGStrategy::BestCost)
should_accept = true;
else if (g_strategy == FGStrategy::BestFit)
should_accept = (num_to_group_0 >= m_max_group_size[0] && num_to_group_1 >= m_max_group_size[1]);
if (!should_accept)
continue; continue;
std::set<int>group_0, group_1; std::set<int>group_0, group_1;
for (int j = 0; j < m_filament_num; ++j) { for (int j = 0; j < m_filament_num; ++j) {
if (i & static_cast<uint64_t>(1 << j)) if (i & static_cast<uint64_t>(1 << j))
@ -52,29 +64,26 @@ namespace Slic3r
group_0.insert(m_used_filaments[j]); group_0.insert(m_used_filaments[j]);
} }
if (group_0.size() < m_max_group_size[0] && group_1.size() < m_max_group_size[1]){ std::vector<int>filament_maps(m_filament_num);
for (int i = 0; i < m_filament_num; ++i) {
if (group_0.find(m_used_filaments[i]) != group_0.end())
filament_maps[i] = 0;
if (group_1.find(m_used_filaments[i]) != group_1.end())
filament_maps[i] = 1;
}
std::vector<int>filament_maps(m_filament_num); int total_cost = reorder_filaments_for_minimum_flush_volume(
for (int i = 0; i < m_filament_num; ++i) { m_used_filaments,
if (group_0.find(m_used_filaments[i]) != group_0.end()) filament_maps,
filament_maps[i] = 0; layer_filaments,
if (group_1.find(m_used_filaments[i]) != group_1.end()) m_flush_matrix,
filament_maps[i] = 1; get_custom_seq,
} nullptr
);
int total_cost = reorder_filaments_for_minimum_flush_volume( if (total_cost < best_cost) {
m_used_filaments, best_cost = total_cost;
filament_maps, best_label = filament_maps;
layer_filaments,
m_flush_matrix,
get_custom_seq,
nullptr
);
if (total_cost < best_cost) {
best_cost = total_cost;
best_label = filament_maps;
}
} }
} }
@ -83,7 +92,7 @@ namespace Slic3r
return best_cost; return best_cost;
} }
int FilamentGroup::calc_filament_group_by_pam(const std::vector<std::vector<unsigned int>>& layer_filaments, int timeout_ms) int FilamentGroup::calc_filament_group_by_pam(const std::vector<std::vector<unsigned int>>& layer_filaments,const FGStrategy& g_strategy, int timeout_ms)
{ {
//calc pair counts //calc pair counts
std::vector<std::vector<int>>count_matrix(m_filament_num,std::vector<int>(m_filament_num)); std::vector<std::vector<int>>count_matrix(m_filament_num,std::vector<int>(m_filament_num));
@ -116,7 +125,7 @@ namespace Slic3r
} }
KMediods PAM(distance_matrix, m_filament_num,m_max_group_size); KMediods PAM(distance_matrix, m_filament_num,m_max_group_size);
PAM.fit(timeout_ms); PAM.fit(g_strategy,timeout_ms);
this->m_filament_labels = PAM.get_filament_labels(); this->m_filament_labels = PAM.get_filament_labels();
int cost = reorder_filaments_for_minimum_flush_volume( int cost = reorder_filaments_for_minimum_flush_volume(
@ -132,7 +141,7 @@ namespace Slic3r
} }
void KMediods::fit( int timeout_ms) void KMediods::fit(const FGStrategy&g_strategy , int timeout_ms)
{ {
std::vector<int>best_medoids; std::vector<int>best_medoids;
std::vector<int>best_labels; std::vector<int>best_labels;
@ -151,7 +160,7 @@ namespace Slic3r
else else
medoids = initialize(INIT_TYPE::Random); medoids = initialize(INIT_TYPE::Random);
labels = assign_label(medoids); labels = assign_label(medoids,g_strategy);
int cost = calc_cost(labels, medoids); int cost = calc_cost(labels, medoids);
for (int i = 0; i < m_filament_num; ++i) { for (int i = 0; i < m_filament_num; ++i) {
@ -161,7 +170,7 @@ namespace Slic3r
for (int j = 0; j < 2; ++j) { for (int j = 0; j < 2; ++j) {
std::vector<int> new_medoids = medoids; std::vector<int> new_medoids = medoids;
new_medoids[j] = i; new_medoids[j] = i;
std::vector<int> new_labels = assign_label(new_medoids); std::vector<int> new_labels = assign_label(new_medoids,g_strategy);
int new_cost = calc_cost(new_labels, new_medoids); int new_cost = calc_cost(new_labels, new_medoids);
if (new_cost < cost) if (new_cost < cost)
@ -182,14 +191,14 @@ namespace Slic3r
} }
count += 1; count += 1;
if (T.time_machine_end() > timeout_ms) if (T.time_machine_end() > timeout_ms || m_medoids_set.size() == (m_filament_num * (m_filament_num - 1) / 2))
break; break;
} }
this->m_filament_labels = best_labels; this->m_filament_labels = best_labels;
} }
std::vector<int> KMediods::assign_label(const std::vector<int>& medoids) const std::vector<int> KMediods::assign_label(const std::vector<int>& medoids,const FGStrategy&g_strategy) const
{ {
std::vector<int>labels(m_filament_num); std::vector<int>labels(m_filament_num);
struct Comp { struct Comp {
@ -205,15 +214,34 @@ namespace Slic3r
min_heap.push({ i,distancec_to_0 - distancec_to_1 }); min_heap.push({ i,distancec_to_0 - distancec_to_1 });
} }
std::set<int> group_0, group_1; std::set<int> group_0, group_1;
while (!min_heap.empty()) { bool have_enough_size = (m_filament_num <= (m_max_group_size[0] + m_max_group_size[1]));
auto top = min_heap.top(); if (have_enough_size || g_strategy == FGStrategy::BestFit) {
min_heap.pop(); while (!min_heap.empty()) {
if (group_0.size() < m_max_group_size[0] && (top.second <= 0 || group_1.size() >= m_max_group_size[1])) auto top = min_heap.top();
group_0.insert(top.first); min_heap.pop();
else if (group_0.size() < m_max_group_size[0] && (top.second <= 0 || group_1.size() >= m_max_group_size[1]))
group_1.insert(top.first); group_0.insert(top.first);
else if (group_1.size() < m_max_group_size[1] && (top.second > 0 || group_0.size() >= m_max_group_size[0]))
group_1.insert(top.first);
else {
if (top.second <= 0)
group_0.insert(top.first);
else
group_1.insert(top.first);
}
}
} }
else if (g_strategy == FGStrategy::BestCost) {
while (!min_heap.empty()) {
auto top = min_heap.top();
min_heap.pop();
if (top.second <= 0)
group_0.insert(top.first);
else
group_1.insert(top.first);
}
}
for (auto& item : group_0) for (auto& item : group_0)
labels[item] = 0; labels[item] = 0;
for (auto& item : group_1) for (auto& item : group_1)

View File

@ -26,9 +26,15 @@ namespace Slic3r
} }
}; };
enum FGStrategy {
BestCost,
BestFit
};
class FilamentGroup class FilamentGroup
{ {
public:
public: public:
FilamentGroup(const std::vector<FlushMatrix>& flush_matrix, const int filament_num, const std::vector<int>& max_group_size) : FilamentGroup(const std::vector<FlushMatrix>& flush_matrix, const int filament_num, const std::vector<int>& max_group_size) :
m_flush_matrix{ flush_matrix }, m_flush_matrix{ flush_matrix },
@ -36,9 +42,9 @@ namespace Slic3r
m_max_group_size{ max_group_size } m_max_group_size{ max_group_size }
{} {}
int calc_filament_group(const std::vector<std::vector<unsigned int>>& layer_filaments); int calc_filament_group(const std::vector<std::vector<unsigned int>>& layer_filaments, const FGStrategy& g_strategy = FGStrategy::BestFit);
int calc_filament_group_by_enum(const std::vector<std::vector<unsigned int>>& layer_filaments); int calc_filament_group_by_enum(const std::vector<std::vector<unsigned int>>& layer_filaments,const FGStrategy& g_strategy);
int calc_filament_group_by_pam(const std::vector<std::vector<unsigned int>>& layer_filaments, int timeout_ms = 300); int calc_filament_group_by_pam(const std::vector<std::vector<unsigned int>>& layer_filaments,const FGStrategy& g_strategy,int timeout_ms = 300);
std::vector<int> get_filament_map() const {return m_filament_labels;} std::vector<int> get_filament_map() const {return m_filament_labels;}
@ -68,14 +74,14 @@ namespace Slic3r
m_filament_num{ filament_num }, m_filament_num{ filament_num },
m_max_group_size{ max_group_size } {} m_max_group_size{ max_group_size } {}
void fit(int timeout_ms = 300); void fit(const FGStrategy& g_strategy,int timeout_ms = 300);
std::vector<int>get_filament_labels()const { std::vector<int>get_filament_labels()const {
return m_filament_labels; return m_filament_labels;
} }
private: private:
std::vector<int>initialize(INIT_TYPE type)const; std::vector<int>initialize(INIT_TYPE type)const;
std::vector<int>assign_label(const std::vector<int>& medoids)const; std::vector<int>assign_label(const std::vector<int>& medoids,const FGStrategy&g_strategy)const;
int calc_cost(const std::vector<int>& labels, const std::vector<int>& medoids)const; int calc_cost(const std::vector<int>& labels, const std::vector<int>& medoids)const;
private: private:
std::vector<std::vector<float>>m_distance_matrix; std::vector<std::vector<float>>m_distance_matrix;

View File

@ -1036,7 +1036,7 @@ std::vector<int> ToolOrdering::get_recommended_filament_maps(const std::vector<s
{ 16,16 } { 16,16 }
); );
fg.get_custom_seq = get_custom_seq; fg.get_custom_seq = get_custom_seq;
fg.calc_filament_group(layer_filaments); fg.calc_filament_group(layer_filaments,FGStrategy::BestFit);
auto filament_map = fg.get_filament_map(); auto filament_map = fg.get_filament_map();
for (size_t idx = 0; idx < filament_map.size(); ++idx) { for (size_t idx = 0; idx < filament_map.size(); ++idx) {