Mercurial > repos > public > sbplib

diff +rv/constructDiffOps.m @ 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
line wrap: on
line diff
--- 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)