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

support plane comparison is exact now

This commit is contained in:
Sven Oesau 2024-04-24 16:59:38 +02:00
parent e9ed069a97
commit c8e45dd0ed
1 changed files with 3 additions and 32 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

35
Kinetic_space_partition/include/CGAL/KSP_3/Support_plane.h
View File

@ -769,38 +769,9 @@ bool operator==(const Support_plane<GeomTraits, IntersectionKernel>& a, const Su
return false; return false;
} }
using FT = typename GeomTraits::FT; if (a.exact_plane() == b.exact_plane())
const auto& planea = a.plane(); return true;
const auto& planeb = b.plane(); else return false;
const auto va = planea.orthogonal_vector();
const auto vb = planeb.orthogonal_vector();
// Are the planes parallel?
FT aval = approximate_angle(va, vb);
CGAL_assertion(aval >= FT(0) && aval <= FT(180));
if (aval >= FT(90))
aval = FT(180) - aval;
if (aval >= a.angle_tolerance()) {
return false;
}
const auto pa1 = a.to_3d(a.centroid());
const auto pb1 = planeb.projection(pa1);
const auto pb2 = b.to_3d(b.centroid());
const auto pa2 = planea.projection(pb2);
const FT bval1 = KSP::internal::distance(pa1, pb1);
const FT bval2 = KSP::internal::distance(pa2, pb2);
const FT bval = (CGAL::max)(bval1, bval2);
CGAL_assertion(bval >= FT(0));
if (bval >= a.distance_tolerance())
return false;
return true;
} }
#endif //DOXYGEN_RUNNING #endif //DOXYGEN_RUNNING