Mercurial > repos > public > sbplib

comparison +time/RungekuttaRV.m @ 1011:e0560bc4fb7d feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add todo:s for time stepping with RV
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 15 Nov 2018 13:49:11 -0800
parents cda996e64925
children
comparison
equal deleted inserted replaced
1010:f753bada1a46 1011:e0560bc4fb7d
35 35
36 function obj = step(obj) 36 function obj = step(obj)
37 obj.v = time.rk.rungekuttaRV(obj.v, obj.t, obj.k, obj.F, obj.RV, obj.coeffs); 37 obj.v = time.rk.rungekuttaRV(obj.v, obj.t, obj.k, obj.F, obj.RV, obj.coeffs);
38 obj.t = obj.t + obj.k; 38 obj.t = obj.t + obj.k;
39 obj.n = obj.n + 1; 39 obj.n = obj.n + 1;
40 % TBD: Add option for updating the residual inside or outside? Decide on best way to do it?
41 % v_prev = obj.v;
42 % F = @(v,t)obj.F(v,t,obj.RV.getViscosity());
43 % obj.v = time.rk.rungekutta(obj.v, obj.t, obj.k, F, obj.coeffs);
44 % obj.RV.update(obj.v,v_prev,obj.k);
45 % obj.t = obj.t + obj.k;
46 % obj.n = obj.n + 1;
40 end 47 end
41 end 48 end
42 49
43 50
44 methods (Static) 51 methods (Static)