49 ynpk1[i] = y[i]+0.5*k1[i];
77 for (
int i=0;i<len;i++)
107 for (i=0;i<y->
len;i++) {
108 k1[i] = dt*ydot->
vec[i];
109 ynpk1[i] = y->
vec[i]+0.5*k1[i];
112 f(ynpk1,ydot,f_data);
114 for (i=0;i<y->
len;i++) {
115 k2[i] = dt*ydot->
vec[i];
136 for (
int i=0;i<y->
len;i++)
137 y->
vec[i] += ydot->
vec[i]*dt;
void ODEint_FE(GlobalData_t *y, GlobalData_t *ydot, void *f, void *f_data, int len, GlobalData_t dt)
void ODEint_RK(GlobalData_t *y, GlobalData_t *ydot, void f(GlobalData_t *, GlobalData_t *, void *), void *f_data, int len, GlobalData_t dt)
void d_ODEint_RK(d_OdeVec *y, d_OdeVec *ydot, void f(GlobalData_t *, d_OdeVec *, void *), void *f_data, GlobalData_t dt)
void d_ODEint_FE(d_OdeVec *y, d_OdeVec *ydot, void *f, void *f_data, GlobalData_t dt)