@@ -94,18 +94,16 @@ class cvodes_integrator_adjoint_vari : public vari_base {
9494 std::vector<Eigen::Matrix<T_Return, Eigen::Dynamic, 1 >> y_return_;
9595 std::tuple<T_Args...> local_args_tuple_;
9696 const std::tuple<
97- promote_scalar_t <partials_type_t <scalar_type_t <T_Args>>, T_Args>...>
98- value_of_args_tuple_;
99-
97+ promote_scalar_t <partials_type_t <scalar_type_t <T_Args>>, T_Args>...>
98+ value_of_args_tuple_;
10099
101100 template <typename FF, typename StateFwd, typename StateBwd, typename Quad,
102101 typename AbsTolFwd, typename AbsTolBwd>
103102 cvodes_solver (const char * function_name, FF&& f, size_t N,
104103 size_t num_args_vars, size_t ts_size, int solver_forward,
105104 StateFwd& state_forward, StateBwd& state_backward, Quad& quad,
106105 AbsTolFwd& absolute_tolerance_forward,
107- AbsTolBwd& absolute_tolerance_backward,
108- const T_Args&... args)
106+ AbsTolBwd& absolute_tolerance_backward, const T_Args&... args)
109107 : chainable_alloc(),
110108 f_ (std::forward<FF>(f)),
111109 function_name_str_(function_name),
@@ -192,7 +190,7 @@ class cvodes_integrator_adjoint_vari : public vari_base {
192190 * @return a vector of states, each state being a vector of the
193191 * same size as the state variable, corresponding to a time in ts.
194192 */
195- template <typename FF, require_eigen_col_vector_t <T_y0>* = nullptr >
193+ template <typename FF, require_eigen_col_vector_t <T_y0>* = nullptr >
196194 cvodes_integrator_adjoint_vari (
197195 const char * function_name, FF&& f, const T_y0& y0, const T_t0& t0,
198196 const std::vector<T_ts>& ts, double relative_tolerance_forward,
@@ -206,8 +204,8 @@ template <typename FF, require_eigen_col_vector_t<T_y0>* = nullptr>
206204 int interpolation_polynomial, int solver_forward, int solver_backward,
207205 std::ostream* msgs, const T_Args&... args)
208206 : vari_base(),
209- // local_args_tuple_(deep_copy_vars(args)...),
210- // value_of_args_tuple_(value_of(args)...),
207+ // local_args_tuple_(deep_copy_vars(args)...),
208+ // value_of_args_tuple_(value_of(args)...),
211209 y_(ts.size()),
212210 ts_(ts.begin(), ts.end()),
213211 y0_(y0),
@@ -597,7 +595,8 @@ template <typename FF, require_eigen_col_vector_t<T_y0>* = nullptr>
597595 */
598596 inline int rhs (double t, const double * y, double *& dy_dt) const {
599597 const Eigen::VectorXd y_vec = Eigen::Map<const Eigen::VectorXd>(y, N_);
600- const Eigen::VectorXd dy_dt_vec = rhs (t, y_vec, solver_->value_of_args_tuple_ );
598+ const Eigen::VectorXd dy_dt_vec
599+ = rhs (t, y_vec, solver_->value_of_args_tuple_ );
601600 check_size_match (solver_->function_name_str_ .c_str (), " dy_dt" ,
602601 dy_dt_vec.size (), " states" , N_);
603602 Eigen::Map<Eigen::VectorXd>(dy_dt, N_) = dy_dt_vec;
0 commit comments