diff --git a/MathLib/Vector3.cpp b/MathLib/Vector3.cpp
index 040b473a9ab7eb2b2e7fc8e779e71cd4f15866d3..f9e80edb2cda9bcf4175c24b8774ee31555fead5 100644
--- a/MathLib/Vector3.cpp
+++ b/MathLib/Vector3.cpp
@@ -12,18 +12,6 @@
 namespace MathLib
 {
 
-double scalarTriple(MathLib::Vector3 const& u, MathLib::Vector3 const& v,
-                    MathLib::Vector3 const& w)
-{
-    auto const pu =
-        Eigen::Map<Eigen::Vector3d>(const_cast<double*>(u.getCoords()));
-    auto const pv =
-        Eigen::Map<Eigen::Vector3d>(const_cast<double*>(v.getCoords()));
-    auto const pw =
-        Eigen::Map<Eigen::Vector3d>(const_cast<double*>(w.getCoords()));
-    return pu.cross(pv).dot(pw);
-}
-
 double scalarTriple(Eigen::Vector3d const& u, Eigen::Vector3d const& v,
                     Eigen::Vector3d const& w)
 {
diff --git a/MathLib/Vector3.h b/MathLib/Vector3.h
index d7b35db9a4a5d125205902db3a1493de186e3bf8..51ccb328613d7e071cb776dfa63a5c197722a1ab 100644
--- a/MathLib/Vector3.h
+++ b/MathLib/Vector3.h
@@ -19,8 +19,8 @@
 #include <Eigen/Eigen>
 #include <cmath>
 
-#include "TemplatePoint.h"
 #include "MathTools.h"
+#include "TemplatePoint.h"
 
 namespace MathLib
 {
@@ -198,10 +198,6 @@ template <typename T1> TemplateVector3<T1> operator*(
 using Vector3 = TemplateVector3<double>;
 
 /// Calculates the scalar triple (u x v) . w
-double scalarTriple(MathLib::Vector3 const& u, MathLib::Vector3 const& v,
-                    MathLib::Vector3 const& w);
-
-/// Calculates the scalar triple (u x v) . w using Eigen functionality
 double scalarTriple(Eigen::Vector3d const& u, Eigen::Vector3d const& v,
                     Eigen::Vector3d const& w);
 }  // namespace MathLib