Mercurial > repos > public > sbplib
comparison +time/Rungekutta4RV.m @ 831:d0934d1143b7 feature/burgers1d
Fix bug in initialization of differential operators
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Tue, 11 Sep 2018 13:24:08 +0200 |
parents | fae41958af4f |
children | 008496ca38f3 |
comparison
equal
deleted
inserted
replaced
815:fae41958af4f | 831:d0934d1143b7 |
---|---|
36 | 36 |
37 function obj = step(obj) | 37 function obj = step(obj) |
38 v_prev = obj.v; | 38 v_prev = obj.v; |
39 F = @(v,t) obj.F(v, t, obj.RV.getViscosity()); | 39 F = @(v,t) obj.F(v, t, obj.RV.getViscosity()); |
40 obj.v = time.rk4.rungekutta_4(obj.v, obj.t, obj.k, F); | 40 obj.v = time.rk4.rungekutta_4(obj.v, obj.t, obj.k, F); |
41 obj.RV.update(obj.v, v_prev, obj.k); | |
41 obj.t = obj.t + obj.k; | 42 obj.t = obj.t + obj.k; |
42 obj.n = obj.n + 1; | 43 obj.n = obj.n + 1; |
43 obj.RV.update(obj.v, v_prev, obj.k); | |
44 end | 44 end |
45 end | 45 end |
46 | 46 |
47 | 47 |
48 methods (Static) | 48 methods (Static) |