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

Assertion if two points are equal

This commit is contained in:
Andreas Fabri 2024-12-09 07:43:24 +00:00
parent 9a837f4a5d
commit b9e4140351
1 changed files with 12 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
Frechet_distance/examples/Frechet_distance/Frechet_distance_3.cpp
View File

@ -1,17 +1,20 @@
#include <CGAL/Frechet_distance.h> #include <CGAL/Frechet_distance.h>
#include <CGAL/Frechet_distance_traits_3.h> #include <CGAL/Frechet_distance_traits_3.h>
#include <CGAL/Simple_cartesian.h> #include <CGAL/Simple_cartesian.h>
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/IO/WKT.h> #include <CGAL/IO/WKT.h>
#include <ostream> #include <ostream>
#include <fstream> #include <fstream>
using Kernel = CGAL::Simple_cartesian<double>; using Kernel = CGAL::Exact_predicates_inexact_constructions_kernel;
//using Kernel = CGAL::Simple_cartesian<double>;
using Traits = CGAL::Frechet_distance_traits_3<Kernel>; using Traits = CGAL::Frechet_distance_traits_3<Kernel>;
using Point = Traits::Point_d; using Point = Traits::Point_d;
int main(int argc, char* argv[]) int main(int argc, char* argv[])
{ {
#if 0
std::vector<Point> A, B; std::vector<Point> A, B;
{ {
std::ifstream in(CGAL::data_file_path("wkt/moebius.wkt")); std::ifstream in(CGAL::data_file_path("wkt/moebius.wkt"));
@ -21,6 +24,14 @@ int main(int argc, char* argv[])
std::ifstream in(CGAL::data_file_path("wkt/moebius2.wkt")); std::ifstream in(CGAL::data_file_path("wkt/moebius2.wkt"));
CGAL::IO::read_linestring_WKT(in, B); CGAL::IO::read_linestring_WKT(in, B);
} }
#else
// two identical points fails
std::array<Point,4> A = { Point(0,0,0), Point(0,0,0), Point(1,0,1), Point(1,1,0)};
std::array<Point,4> B = { Point(0,0,0), Point(0,0,0), Point(1,0,0), Point(1,1,0)};
#endif
std::pair<double, double> res = CGAL::approximate_Frechet_distance<Traits>(A, B, 0.000001); std::pair<double, double> res = CGAL::approximate_Frechet_distance<Traits>(A, B, 0.000001);
std::cout << "The Frechet distance between the polylines is between " << res.first << " and " << res.second << std::endl; std::cout << "The Frechet distance between the polylines is between " << res.first << " and " << res.second << std::endl;
return 0; return 0;