[MaL] ODE: Add MappedVector/Matrix template alias.

MathLib/ODE/OdeSolver.h

@@ -28,8 +28,7 @@ class OdeSolver
 {
 public:
 	using Arr = std::array<double, NumEquations>;
-	using ConstArrRef =
-	    Eigen::Map<const Eigen::Matrix<double, NumEquations, 1>>;
+	using ConstArrRef = MappedConstVector<NumEquations>;
 	using Function = MathLib::Function<NumEquations, FunctionArguments...>;
 	using JacobianFunction =
 	    MathLib::JacobianFunction<NumEquations, FunctionArguments...>;

MathLib/ODE/OdeSolverFactory.h

@@ -44,10 +44,7 @@ struct Handles<N, FunctionArgument> : public MathLib::FunctionHandles
 		// consider omission of data pointer and switch to std::function or
 		// alike
 		if (f)
-			return f(t,
-			         Eigen::Map<const Eigen::Matrix<double, N, 1>>{y},
-			         Eigen::Map<Eigen::Matrix<double, N, 1>>{ydot},
-			         _data);
+			return f(t, MappedConstVector<N>{y}, MappedVector<N>{ydot}, _data);
 		return true;
 	}
 
@@ -56,9 +53,9 @@ struct Handles<N, FunctionArgument> : public MathLib::FunctionHandles
 	{
 		if (df)
 			return df(t,
-			          Eigen::Map<const Eigen::Matrix<double, N, 1>>{y},
-			          Eigen::Map<Eigen::Matrix<double, N, 1>>{ydot},
-			          Eigen::Map<Eigen::Matrix<double, N, N>>{jac /*, order*/},
+			          MappedConstVector<N>{y},
+			          MappedVector<N>{ydot},
+			          MappedMatrix<N, N>{jac /*, order*/},
 			          _data);
 		return true;
 	}
@@ -85,9 +82,7 @@ struct Handles<N> : public MathLib::FunctionHandles
 		{
 			// auto ydot_ = Eigen::Map<Eigen::Matrix<double, N, 1>>{y};
 			// auto ydot_ = Eigen::Map<Eigen::Matrix<double, N, 1>>{ydot};
-			return f(t,
-			         Eigen::Map<const Eigen::Matrix<double, N, 1>>{y},
-			         Eigen::Map<Eigen::Matrix<double, N, 1>>{ydot});
+			return f(t, MappedConstVector<N>{y}, MappedVector<N>{ydot});
 		}
 		return true;
 	}
@@ -97,9 +92,9 @@ struct Handles<N> : public MathLib::FunctionHandles
 	{
 		if (df)
 			return df(t,
-			          Eigen::Map<const Eigen::Matrix<double, N, 1>>{y},
-			          Eigen::Map<Eigen::Matrix<double, N, 1>>{ydot},
-			          Eigen::Map<Eigen::Matrix<double, N, N>>{jac /*, order*/});
+			          MappedConstVector<N>{y},
+			          MappedVector<N>{ydot},
+			          MappedMatrix<N, N>{jac /*, order*/});
 		return true;
 	}
 

MathLib/ODE/declarations.h

@@ -14,19 +14,30 @@
 
 namespace MathLib
 {
+template <int M, int N>
+using MappedMatrix = Eigen::Map<Eigen::Matrix<double, M, N>>;
+
+template <int M, int N>
+using MappedConstMatrix = Eigen::Map<const Eigen::Matrix<double, M, N>>;
+
+template <int N>
+using MappedVector = MappedMatrix<N, 1>;
+
+template <int N>
+using MappedConstVector = MappedConstMatrix<N, 1>;
+
 template <unsigned N, typename... FunctionArguments>
 using Function = bool (*)(const double t,
-                          Eigen::Map<const Eigen::Matrix<double, N, 1>> const y,
-                          Eigen::Map<Eigen::Matrix<double, N, 1>> ydot,
+                          MappedConstVector<N> const y,
+                          MappedVector<N> ydot,
                           FunctionArguments&... arg);
 
 template <unsigned N, typename... FunctionArguments>
-using JacobianFunction =
-    bool (*)(const double t,
-             Eigen::Map<const Eigen::Matrix<double, N, 1>> const y,
-             Eigen::Map<Eigen::Matrix<double, N, 1>> ydot,
-             Eigen::Map<Eigen::Matrix<double, N, N>> jac,
-             FunctionArguments&... arg);
+using JacobianFunction = bool (*)(const double t,
+                                  MappedConstVector<N> const y,
+                                  MappedVector<N> ydot,
+                                  MappedMatrix<N, N> jac,
+                                  FunctionArguments&... arg);
 
 // This is an internal detail
 class FunctionHandles

Tests/MathLib/TestCVode.cpp

@@ -16,8 +16,8 @@
 #include <cstdio>
 
 bool f(const double,
-       Eigen::Map<const Eigen::Matrix<double, 1, 1>> const y,
-       Eigen::Map<Eigen::Matrix<double, 1, 1>> ydot)
+       MathLib::MappedConstVector<1> const y,
+       MathLib::MappedVector<1> ydot)
 {
 	if (y[0] <= 0.0) return false;
 
@@ -26,9 +26,9 @@ bool f(const double,
 }
 
 bool df(const double /*t*/,
-        Eigen::Map<const Eigen::Matrix<double, 1, 1>> const y,
-        Eigen::Map<Eigen::Matrix<double, 1, 1>> /*ydot*/,
-        Eigen::Map<Eigen::Matrix<double, 1, 1>> jac)
+        MathLib::MappedConstVector<1> const y,
+        MathLib::MappedVector<1> /*ydot*/,
+        MathLib::MappedMatrix<1, 1> jac)
 {
 	if (y[0] <= 0.0) return false;
 
@@ -42,8 +42,8 @@ struct ExtraData
 };
 
 bool f_extra(const double,
-             Eigen::Map<const Eigen::Matrix<double, 1, 1>> const y,
-             Eigen::Map<Eigen::Matrix<double, 1, 1>> ydot,
+             MathLib::MappedConstVector<1> const y,
+             MathLib::MappedVector<1> ydot,
              ExtraData& data)
 {
 	if (y[0] <= 0.0) return false;