Homec4science
Diffusion SpecMiCP / ReactMiCP (gasphase) Browse

specmicp/src/odeint/odeint_types.hpp
881a9cd5b582gasphase

Last Change
Raw File

odeint_types.hpp
View Options

/*-------------------------------------------------------
- Module : odeint
- File : odeint_types.hpp
- Author : Fabien Georget
Copyright (c) 2014, Fabien Georget, Princeton University
---------------------------------------------------------*/
#ifndef SPECMICP_ODEINT_ODEINTTYPES_HPP
#define SPECMICP_ODEINT_ODEINTTYPES_HPP
//! \file odeint_types.hpp common types for odeint
#include <functional>
#include <Eigen/Dense>
namespace specmicp {
namespace odeint {
//! A scalar
using scalar_t = double;
namespace internal {
//! \brief implementation for the vector type (which may not be a vector for dim=1)
template <int dim>
struct vector_trait
{
using type = Eigen::Matrix<scalar_t, Eigen::Dynamic, 1>;
};
//! Specialisation for dimension 1 (scalar)
template <>
struct vector_trait<1>
{
using type = scalar_t;
};
//! Specialization for dimension 2 (fixed-size vector)
template <>
struct vector_trait<2>
{
using type = Eigen::Matrix<scalar_t, 2, 1>;
};
//! Specialization for dimension 3 (fixed-size vector)
template <>
struct vector_trait<3>
{
using type = Eigen::Matrix<scalar_t, 3, 1>;
};
} // end namespace internal
//! The vector type
template <int dim>
using vector_t = typename internal::vector_trait<dim>::type;
//! Return the maximum coefficient of the vector
template <int dim>
scalar_t max_coeff(const vector_t<dim>& vec) {return vec.maxCoeff();}
template <>
scalar_t max_coeff<1>(const vector_t<1>& vec) {return vec;}
//! Type of the function f in dy/dx=f(x,y)
template <int dim>
using rhs_f = std::function<void (scalar_t, const vector_t<dim>&, vector_t<dim>&)>;
} // end namespace odeint
} // end namespace specmicp
#endif // SPECMICP_ODEINT_ODEINTTYPES_HPP
c4science · Help