Mercurial > repos > public > sbplib

diff +time/Rungekutta4RV.m @ 835:008496ca38f3 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Compute the residual in between each runge-kutta stage. Note: It is not clear whether the correct residual is used when computing the stages. Must investigate further.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 13 Sep 2018 18:14:54 +0200
parents d0934d1143b7
children 1e057b0f2fed
line wrap: on
line diff
--- a/+time/Rungekutta4RV.m	Thu Sep 13 14:21:47 2018 +0200
+++ b/+time/Rungekutta4RV.m	Thu Sep 13 18:14:54 2018 +0200
@@ -28,17 +28,13 @@
             t = obj.t;
         end
 
-        function [residual, viscosity, t] = getRV(obj)
-            residual = obj.RV.getResidual();
-            viscosity = obj.RV.getViscosity();
-            t = obj.t;
+        function state = getState(obj)
+            [residual, u_t, f_x] = obj.RV.getResidual();
+            state = struct('v', obj.v, 'residual', residual, 'u_t', u_t, 'f_x', f_x, 'viscosity', obj.RV.getViscosity(), 't', obj.t);
         end
 
         function obj = step(obj)
-            v_prev = obj.v;
-            F = @(v,t) obj.F(v, t, obj.RV.getViscosity());
-            obj.v = time.rk4.rungekutta_4(obj.v, obj.t, obj.k, F);
-            obj.RV.update(obj.v, v_prev, obj.k);
+            obj.v = time.rk4.rungekutta_4RV(obj.v, obj.t, obj.k, obj.F, obj.RV);
             obj.t = obj.t + obj.k;
             obj.n = obj.n + 1;
         end