diff --git a/MathLib/ODE/CVodeSolver.cpp b/MathLib/ODE/CVodeSolver.cpp
index 6b9cb14f2971ec01525980755f4ff825790dea64..c6f242809394fff8dcbc6c6c87e1edd737592488 100644
--- a/MathLib/ODE/CVodeSolver.cpp
+++ b/MathLib/ODE/CVodeSolver.cpp
@@ -256,8 +256,7 @@ void CVodeSolverImpl::preSolve()
 			bool successful =
 			    static_cast<FunctionHandles*>(function_handles)
 			        ->callJacobian(t, NV_DATA_S(y), NV_DATA_S(ydot),
-			                       DENSE_COL(jac, 0),
-			                       MathLib::StorageOrder::ColumnMajor);
+			                       DENSE_COL(jac, 0));
 			return successful ? 0 : 1;
 		};
 
diff --git a/MathLib/ODE/OdeSolverFactory.h b/MathLib/ODE/OdeSolverFactory.h
index 5b9d7c7beec3f179125914c93c3e2979e3f18559..4330db613c75a061012b4c1b6421bba672f8632c 100644
--- a/MathLib/ODE/OdeSolverFactory.h
+++ b/MathLib/ODE/OdeSolverFactory.h
@@ -46,9 +46,8 @@ struct Handles<N, FunctionArgument> : public MathLib::FunctionHandles
 		return true;
 	}
 
-	bool callJacobian(const double t, const double* const y,
-	                  const double* const ydot, double* const jac,
-	                  MathLib::StorageOrder order) override
+	bool callJacobian(const double t, const double* const y, double* const ydot,
+	                  double* const jac) override
 	{
 		if (df)
 			return df(t,
@@ -95,9 +94,8 @@ struct Handles<N> : public MathLib::FunctionHandles
 		return true;
 	}
 
-	bool callJacobian(const double t, const double* const y,
-	                  const double* const ydot, double* const jac,
-	                  MathLib::StorageOrder order) override
+	bool callJacobian(const double t, const double* const y, double* const ydot,
+	                  double* const jac) override
 	{
 		if (df)
 			return df(t,
diff --git a/MathLib/ODE/declarations.h b/MathLib/ODE/declarations.h
index 000f4246a4a6e10350ff82aaa259c47090252d2d..753e56ef393dff320319b18c75e1b567e5c9933c 100644
--- a/MathLib/ODE/declarations.h
+++ b/MathLib/ODE/declarations.h
@@ -14,12 +14,6 @@
 
 namespace MathLib
 {
-enum class StorageOrder
-{
-	ColumnMajor,
-	RowMajor
-};
-
 template <unsigned N, typename... FunctionArguments>
 using Function = bool (*)(const double t,
                           Eigen::Map<const Eigen::Matrix<double, N, 1>> const y,
@@ -43,8 +37,7 @@ public:
 	virtual bool callJacobian(const double t,
 	                          double const* const y,
 	                          double* const ydot,
-	                          double* const jac,
-	                          StorageOrder order) = 0;
+	                          double* const jac) = 0;
 
 	virtual bool hasJacobian() const = 0;