[MaL] Move scalarTriple to MathTools.

MathLib/MathTools.cpp

@@ -51,4 +51,10 @@ double getAngle (const double p0[3], const double p1[3], const double p2[3])
     return std::acos (scalarProduct<double,3> (v0,v1) / (std::sqrt(scalarProduct<double,3>(v0,v0)) * sqrt(scalarProduct<double,3>(v1,v1))));
 }
 
+double scalarTriple(Eigen::Vector3d const& u, Eigen::Vector3d const& v,
+                    Eigen::Vector3d const& w)
+{
+    return u.cross(v).dot(w);
+}
+
 }  // namespace MathLib

MathLib/MathTools.h

@@ -90,4 +90,8 @@ double calcProjPntToLineAndDists(MathLib::Point3d const& pp,
  * @return the angle between the edges
  */
 double getAngle (const double p0[3], const double p1[3], const double p2[3]);
+
+/// Calculates the scalar triple (u x v) . w
+double scalarTriple(Eigen::Vector3d const& u, Eigen::Vector3d const& v,
+                    Eigen::Vector3d const& w);
 }  // namespace MathLib

MathLib/Vector3.cpp

-/**
- * \file
- * \copyright
- * Copyright (c) 2012-2020, OpenGeoSys Community (http://www.opengeosys.org)
- *            Distributed under a Modified BSD License.
- *              See accompanying file LICENSE.txt or
- *              http://www.opengeosys.org/project/license
- */
-
-#include "Vector3.h"
-
-namespace MathLib
-{
-
-double scalarTriple(Eigen::Vector3d const& u, Eigen::Vector3d const& v,
-                    Eigen::Vector3d const& w)
-{
-    return u.cross(v).dot(w);
-}
-
-}  // end namespace MathLib

MathLib/Vector3.h

@@ -196,8 +196,4 @@ template <typename T1> TemplateVector3<T1> operator*(
 }
 
 using Vector3 = TemplateVector3<double>;
-
-/// Calculates the scalar triple (u x v) . w
-double scalarTriple(Eigen::Vector3d const& u, Eigen::Vector3d const& v,
-                    Eigen::Vector3d const& w);
 }  // namespace MathLib