diff --git a/GeoLib/AnalyticalGeometry.cpp b/GeoLib/AnalyticalGeometry.cpp
index 1a22d47307a4b1dac26c296987a76435a027d309..dd029dc5bb1e098afed8d6e3b8f000e6edcbf0c3 100644
--- a/GeoLib/AnalyticalGeometry.cpp
+++ b/GeoLib/AnalyticalGeometry.cpp
@@ -44,9 +44,9 @@ double getOrientation2d(MathLib::Point3d const& a,
 
 namespace GeoLib
 {
-Orientation getOrientation(GeoLib::Point const& p0,
-                           GeoLib::Point const& p1,
-                           GeoLib::Point const& p2)
+Orientation getOrientation(MathLib::Point3d const& p0,
+                           MathLib::Point3d const& p1,
+                           MathLib::Point3d const& p2)
 {
     double const orientation = ExactPredicates::getOrientation2d(p0, p1, p2);
     if (orientation > 0)
diff --git a/GeoLib/AnalyticalGeometry.h b/GeoLib/AnalyticalGeometry.h
index 08b5c8d6b1ba9286324c7d36052579e403e5e82d..d9f82309803495489aba166d2e69794613694783 100644
--- a/GeoLib/AnalyticalGeometry.h
+++ b/GeoLib/AnalyticalGeometry.h
@@ -41,10 +41,9 @@ enum Orientation
  * \returns CW (clockwise), CCW (counterclockwise) or COLLINEAR (points are on a
  * line)
  */
-Orientation getOrientation(GeoLib::Point const& p0,
-                           GeoLib::Point const& p1,
-                           GeoLib::Point const& p2);
-
+Orientation getOrientation(MathLib::Point3d const& p0,
+                           MathLib::Point3d const& p1,
+                           MathLib::Point3d const& p2);
 /**
  * compute a supporting plane (represented by plane_normal and the value d) for the polygon
  * Let \f$n\f$ be the plane normal and \f$d\f$ a parameter. Then for all points \f$p \in R^3\f$ of the plane