Mercurial > repos > public > sbplib

comparison +rv/+time/RungekuttaInteriorRv.m @ 1030:78c75c95b7dd feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Rename Rungekutta-Rv classes according to naming convention
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 11 Jan 2019 15:52:48 +0100
parents +rv/+time/RungekuttaInteriorRV.m@2d7c1333bd6c
children 2ef20d00b386
comparison
equal deleted inserted replaced
1029:dce08a74e0ad 1030:78c75c95b7dd
1 classdef RungekuttaInteriorRv < time.Timestepper
2 properties
3 F % RHS of the ODE
4 k % Time step
5 t % Time point
6 v % Solution vector
7 n % Time level
8 coeffs % The coefficents used for the RK time integration
9 RV % Residual Viscosity
10 viscosity % Viscosity vector
11 DvDt % Function for computing the time deriative used for the RV evaluation
12
13 residual
14 dvdt
15 Df
16 end
17
18 methods
19 function obj = RungekuttaInteriorRv(F, k, t0, v0, RV, DvDt, order)
20 obj.F = F;
21 obj.k = k;
22 obj.t = t0;
23 obj.v = v0;
24 obj.n = 0;
25 % Extract the coefficients for the specified order
26 % used for the RK updates from the Butcher tableua.
27 [s,a,b,c] = time.rk.butcherTableau(order);
28 obj.coeffs = struct('s',s,'a',a,'b',b,'c',c);
29 obj.RV = RV;
30 obj.DvDt = DvDt;
31 obj.dvdt = obj.DvDt(obj.v);
32 [obj.viscosity, obj.Df] = obj.RV.evaluate(obj.v,obj.dvdt);
33 obj.residual = obj.dvdt + obj.Df;
34 end
35
36 function [v, t] = getV(obj)
37 v = obj.v;
38 t = obj.t;
39 end
40
41 function state = getState(obj)
42 state = struct('v', obj.v, 'residual', obj.residual, 'dvdt', obj.dvdt, 'Df', obj.Df, 'viscosity', obj.viscosity, 't', obj.t);
43 end
44
45 function obj = step(obj)
46 obj.v = rv.time.rungekuttaRV(obj.v, obj.t, obj.k, obj.F, obj.RV, obj.DvDt, obj.coeffs);
47 obj.t = obj.t + obj.k;
48 obj.n = obj.n + 1;
49 obj.dvdt = obj.DvDt(obj.v);
50 [obj.viscosity, obj.Df] = obj.RV.evaluate(obj.v,obj.dvdt);
51 obj.residual = obj.dvdt + obj.Df;
52 end
53 end
54 end