1 #ifndef STAN_MATH_REV_MAT_FUNCTOR_INTEGRATOR_DAE_HPP 2 #define STAN_MATH_REV_MAT_FUNCTOR_INTEGRATOR_DAE_HPP 36 template <
typename F,
typename Tpar>
38 const F& f,
const std::vector<double>& yy0,
const std::vector<double>& yp0,
39 double t0,
const std::vector<double>& ts,
const std::vector<Tpar>& theta,
40 const std::vector<double>& x_r,
const std::vector<int>& x_i,
41 const double rtol,
const double atol,
43 std::ostream* msgs =
nullptr) {
46 const std::vector<int> dummy_eq_id(yy0.size(), 0);
50 f, dummy_eq_id, yy0, yp0, theta, x_r, x_i, msgs};
54 return solver.integrate(dae, t0, ts);
IDAS DAE system with forward sensitivity calculation.
static constexpr int IDAS_MAX_STEPS
void check_ic_consistency(const double &t0, const double &tol)
std::vector< std::vector< Tpar > > integrate_dae(const F &f, const std::vector< double > &yy0, const std::vector< double > &yp0, double t0, const std::vector< double > &ts, const std::vector< Tpar > &theta, const std::vector< double > &x_r, const std::vector< int > &x_i, const double rtol, const double atol, const int64_t max_num_steps=idas_integrator::IDAS_MAX_STEPS, std::ostream *msgs=nullptr)
Return the solutions for a semi-explicit DAE system with residual specified by functor F...