ODESolverTypes.h

Go to the documentation of this file.

 
#pragma once
 
#include <functional>
#include <Eigen/Core>
 
namespace MathLib
{
namespace ODE
{
 
template <int N, int M>
using MappedMatrix = Eigen::Map<Eigen::Matrix<double, N, M, Eigen::ColMajor>>;
 
template <int N, int M>
using MappedConstMatrix =
    Eigen::Map<const Eigen::Matrix<double, N, M, Eigen::ColMajor>>;
 
template <int N>
using MappedVector = MappedMatrix<N, 1>;
 
template <int N>
using MappedConstVector = MappedConstMatrix<N, 1>;
 
template <unsigned N>
using Function = std::function<bool(
    const double t, MappedConstVector<N> const& y, MappedVector<N>& ydot)>;
 
template <unsigned N>
using JacobianFunction = std::function<bool(const double t,
                                            MappedConstVector<N> const& y,
                                            MappedConstVector<N> const& ydot,
                                            MappedMatrix<N, N>& jac)>;
 
 
}  // namespace ODE
}  // namespace MathLib