[MathLib/GeoLib] Made operator== a template function.

GeoLib/Point.cpp

@@ -17,8 +17,6 @@
 
 #include "Point.h"
 
-#include "MathTools.h"
-
 namespace GeoLib {
 
 bool operator<= (const GeoLib::Point& p0, const GeoLib::Point& p1)
@@ -77,11 +75,4 @@ bool lessEq(const GeoLib::Point& p0, const GeoLib::Point& p1, double tol)
 	}
 }
 
-bool operator==(GeoLib::Point const& a, GeoLib::Point const& b)
-{
-	double sqr_dist(MathLib::sqrDist(a,b));
-	return (sqr_dist < sqrt(std::numeric_limits<double>::min()));
-}
-
-
 } // end namespace GeoLib

GeoLib/Point.h

@@ -70,9 +70,6 @@ bool operator<= (GeoLib::Point const & p0, GeoLib::Point const & p1);
 bool lessEq(const GeoLib::Point& p0,
             const GeoLib::Point& p1,
             double tol = std::numeric_limits<double>::epsilon());
-
-bool operator==(GeoLib::Point const& a, GeoLib::Point const& b);
-
 }
 
 #endif /* POINT_H_ */

MathLib/MathTools.h

@@ -121,6 +121,13 @@ typename POINT_T::FP_T sqrDist(POINT_T const& p0, POINT_T const& p1)
 	return MathLib::scalarProduct<typename POINT_T::FP_T,3>(v,v);
 }
 
+template <typename POINT_T>
+bool operator==(POINT_T const& a, POINT_T const& b)
+{
+	typename POINT_T::FP_T const sqr_dist(sqrDist(a,b));
+	return (sqr_dist < pow(std::numeric_limits<typename POINT_T::FP_T>::epsilon(),2));
+}
+
 /** squared dist between double arrays p0 and p1 (size of arrays is 3) */
 double sqrDist(const double* p0, const double* p1);