size_t integrate_const(
Stepper stepper , System system , State &start_state ,
Time start_time , Time end_time , Time dt
)
{
return integrate_const( stepper , system , start_state , start_time , end_time , dt , null_observer() );
}
size_t integrate_const(
Stepper &stepper ,
DynamicalSystem &system ,
typename Stepper::container_type &state ,
typename Stepper::time_type start_time ,
typename Stepper::time_type end_time ,
typename Stepper::time_type dt
)
{
return integrate_const(
stepper , system , state, start_time , end_time , dt ,
do_nothing_observer<
typename Stepper::time_type ,
typename Stepper::container_type ,
DynamicalSystem >
);
}
ode_out sim_cpp (const NumericVector Ainit, double t_start, double t_end, double step_size) {
int n_steps = (int) ceil((t_end-t_start)/step_size)+1;
double t_end_new = t_start + ((n_steps-1) * step_size);
vector<state_type> x_vec;
vector<double> tim;
NumericMatrix A_ret (n_steps, n_comp+2);
std::fill(A_ret.begin(), A_ret.end(), 0);
state_type A = {} ;
for(int j = 0; j < n_comp; j++) {
A[j] = Ainit(j);
}
ode_out tmp;
integrate_const (stepr , ode, A, t_start, t_end_new, step_size, push_back_solution (x_vec, tim));
tmp.y = x_vec;
tmp.time = tim;
return(tmp);
}
