ode_store_sensitivities.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_PRIM_FUNCTOR_ODE_STORE_SENSITIVITIES_HPP
#define STAN_MATH_PRIM_FUNCTOR_ODE_STORE_SENSITIVITIES_HPP
 
#include <stan/math/prim/meta.hpp>
#include <stan/math/prim/fun/Eigen.hpp>
#include <ostream>
#include <vector>
 
namespace stan {
namespace math {
 
template <
    typename F, typename T_y0_t0, typename T_t0, typename T_t, typename... Args,
    typename
    = require_all_arithmetic_t<T_y0_t0, T_t0, T_t, scalar_type_t<Args>...>>
Eigen::VectorXd ode_store_sensitivities(
    const F& f, const std::vector<double>& coupled_state,
    const Eigen::Matrix<T_y0_t0, Eigen::Dynamic, 1>& y0, T_t0 t0, T_t t,
    std::ostream* msgs, const Args&... args) {
  return Eigen::Map<const Eigen::VectorXd>(coupled_state.data(),
                                           coupled_state.size());
}
 
}  // namespace math
}  // namespace stan
#endif