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

Fix assertion: global lbound is init. w/ the initial bound, not the error bound

This commit is contained in:
Mael Rouxel-Labbé 2022-03-31 12:50:05 +02:00
parent f46d4d1fae
commit 14a9abcca6
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/distance.h
View File

@ -1733,14 +1733,15 @@ bounded_error_squared_Hausdorff_distance_impl(const TriangleMesh1& tm1,
std::cout << "Top of the queue is lower than the lowest!" << std::endl; std::cout << "Top of the queue is lower than the lowest!" << std::endl;
#endif #endif
// This can happen if the actual Hausdorff distance is smaller than sq_initial_bound,
// which can happen either if the error bound is really large, or if the initial bound
// is non-zero (symmetric Hausdorff distance). In this case, we return 'sq_initial_bound'
CGAL_assertion(global_bounds.lower == sq_initial_bound);
global_bounds.upper = global_bounds.lower; // not really needed since lower is returned but doesn't hurt global_bounds.upper = global_bounds.lower; // not really needed since lower is returned but doesn't hurt
global_bounds.upair.first = global_bounds.lpair.first; global_bounds.upair.first = global_bounds.lpair.first;
global_bounds.upair.second = global_bounds.lpair.second; global_bounds.upair.second = global_bounds.lpair.second;
// Current upmost being equal to the lower is fine, but if it's strictly below, it must
// be because we crossed the error bound, or there is some issue...
CGAL_assertion(global_bounds.lower == sq_error_bound);
break; break;
} }