Mercurial > repos > public > sbplib
diff +time/+rk/rungekutta_4RV.m @ 846:c6fcee3fcf1b feature/burgers1d
Add generalized RungeKutta and RungeKuttaRV class which extracts its coefficients from a butcher tableau, specified on the scheme.
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 20 Sep 2018 17:51:19 +0200 |
| parents | +time/+rk4/rungekutta_4RV.m@f63b99f0729d |
| children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+time/+rk/rungekutta_4RV.m Thu Sep 20 17:51:19 2018 +0200 @@ -0,0 +1,18 @@ +% Takes one time step of size k using the rungekutta method +% starting from v_0 and where the function F(v,t) gives the +% time derivatives. +function v = rungekutta_4RV(v, t , k, F, RV) + k1 = F(v, t, RV.getViscosity()); + + RV.update(v+0.5*k*k1, v, 0.5*k); + k2 = F(v+0.5*k*k1, t+0.5*k, RV.getViscosity()); + + RV.update(v+0.5*k*k2, v, 0.5*k); + k3 = F(v+0.5*k*k2, t+0.5*k, RV.getViscosity()); + + RV.update(v+k*k3, v, k); + k4 = F(v+k*k3,t+k, RV.getViscosity()); + + RV.update(v + (1/6)*(k1+2*(k2+k3)+k4)*k, v, k); + v = v + (1/6)*(k1+2*(k2+k3)+k4)*k; +end \ No newline at end of file