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

Use interval arithmetic where needed

This commit is contained in:
Andreas Fabri 2025-01-28 12:01:05 +00:00
parent 35af188a03
commit 415645e8eb
2 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Frechet_distance/examples/Frechet_distance/Frechet_DS_2.cpp
View File

@ -21,7 +21,7 @@ int main()
const std::filesystem::path data{"./data_2d"};
std::vector<std::string> filenames;
for (auto const& dir_entry : std::filesystem::directory_iterator{data}) {
filenames.push_back(dir_entry.path());
filenames.push_back(dir_entry.path().string());
}
std::sort(filenames.begin(), filenames.end());

10
Frechet_distance/include/CGAL/Frechet_distance/internal/Neighbor_search.h
View File

@ -143,7 +143,8 @@ private:
for (size_t i = 0; i < 2*dim; ++i, ++pb, ++qb) {
// AF: certainly
// AN: yes, here we need certainly!
if (CGAL::abs(*pb - *qb) > distance) {
CGAL::Interval_nt<true> ip = to_interval(*pb), iq = to_interval(*qb);
if (certainly(CGAL::abs(ip - iq) > distance)) {
return false;
}
}
@ -154,10 +155,9 @@ private:
{
typename Point_d::Cartesian_const_iterator_d pb = p.cartesian_begin();
for (int d = 0; d < D::value; ++d, ++pb) {
if (rect.min_coord(d) > *pb + distance &&
rect.max_coord(d) + distance < *pb) {
// AF: certainly
// AN: yes, here is really certainly!
CGAL::Interval_nt<true> irmin = to_interval(rect.min_coord(d)), ip = to_interval(*pb), irmax = to_interval(rect.max_coord(d));
if (certainly( irmin> ip + distance &&
irmax + distance < ip)){
return false;
}
}