Mercurial > repos > public > sbplib

changeset 1161:856bd6291d17 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use same order for both time derivative and flux term when computing the residual
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 25 Jun 2019 17:00:27 +0200
parents 76e3bb7836cf
children 0ec06ca3fc36
files +rv/ResidualViscosity.m +rv/constructDiffOps.m
diffstat 2 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+rv/ResidualViscosity.m	Tue Jun 25 16:52:54 2019 +0200
+++ b/+rv/ResidualViscosity.m	Tue Jun 25 17:00:27 2019 +0200
@@ -18,7 +18,7 @@
         %       wrapping it in a function.
         function obj = ResidualViscosity(grid, Df, waveSpeed, Cmax, Cres, h, normalization, postProcess)
             default_arg('normalization',@(v)abs(obj.minmaxDiffNeighborhood1d(v)-norm(v-mean(v),inf)));
-            default_arg('postProcess','none');
+            default_arg('postProcess','maximum neighbors');
             obj.Df = Df;
             obj.waveSpeed = waveSpeed;
             obj.h = h;
--- a/+rv/constructDiffOps.m	Tue Jun 25 16:52:54 2019 +0200
+++ b/+rv/constructDiffOps.m	Tue Jun 25 17:00:27 2019 +0200
@@ -18,7 +18,7 @@
     % change the sign.
     D_flux = @(v) -D_res(v);
     % DiffOp for time derivative in residual viscosity
-    DvDt = D;
+    DvDt = D_res;
 end
 
 function [D, penalties] = constructFluxDiffOp(scheme, g, order, schemeParams, opSet, BCs)