Mercurial > repos > public > sbplib

changeset 831:d0934d1143b7 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 5573913a0949
files +scheme/Burgers1D.m +time/Rungekutta4RV.m
diffstat 2 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Burgers1D.m	Thu Sep 06 12:43:51 2018 +0200
+++ b/+scheme/Burgers1D.m	Tue Sep 11 13:24:08 2018 +0200
@@ -33,9 +33,9 @@
 
             switch pde_form
                 case 'skew-symmetric'
-                    D = @(v, viscosity) -1/3*v.*D1*v - 1/3*D1*v.^2 + D2(obj.params.eps + viscosity)*v;
+                    D = @(v, viscosity) -1/3*v.*D1*v - 1/3*D1*v.^2 + D2(params.eps + viscosity)*v;
                 case 'conservative'
-                    D = @(v, viscosity) -1/2*D1*v.^2 + D2(obj.params.eps + viscosity)*v;
+                    D = @(v, viscosity) -1/2*D1*v.^2 + D2(params.eps + viscosity)*v;
             end
 
             obj.m = m;
--- a/+time/Rungekutta4RV.m	Thu Sep 06 12:43:51 2018 +0200
+++ b/+time/Rungekutta4RV.m	Tue Sep 11 13:24:08 2018 +0200
@@ -38,9 +38,9 @@
             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.t = obj.t + obj.k;
             obj.n = obj.n + 1;
-            obj.RV.update(obj.v, v_prev, obj.k);
         end
     end