songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

prevent linearity in median splitters due to duplicated points

This commit is contained in:
Sven Oesau 2025-01-30 16:55:51 +01:00
parent 5b28489b16
commit 48569ff99c
1 changed files with 18 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
Spatial_searching/include/CGAL/Point_container.h
View File

@ -419,6 +419,24 @@ public:
typename Traits::Cartesian_const_iterator_d mpit = construct_it((*(*mid)));
FT val1 = *(mpit+split_coord);
// Avoid using the low coord value as it results in an empty split
if (val1 == tbox.min_coord(split_coord)) {
iterator it = std::min_element(mid, end(), [=](const Point_d* a, const Point_d* b) -> bool {
FT a_c = *(construct_it(*a) + split_coord);
FT b_c = *(construct_it(*b) + split_coord);
if (a_c == val1)
return false;
if (b_c == val1)
return true;
return a_c < b_c;
});
return *(construct_it(**it) + split_coord);
}
mid++;
mpit = construct_it((*(*mid)));
FT val2 = *(mpit+split_coord);