--- a/+rv/constructDiffOps.m	Tue Jun 25 15:04:53 2019 +0200
+++ b/+rv/constructDiffOps.m	Tue Jun 25 16:52:54 2019 +0200
@@ -5,7 +5,7 @@
     D_rv = @(v,Viscosity)(D(v) + D2(Viscosity)*v);
 
     %% DiffOps for residual viscosity
-    [D_flux, residualPenalties] = constructFluxDiffOp(scheme, g, residualOrder, schemeParams, opSet, BCs);
+    [D_res, residualPenalties] = constructFluxDiffOp(scheme, g, residualOrder, schemeParams, opSet, BCs);
     % TODO: Construct D_flux without closures when using bdfs.
     % diffOp = scheme(g, residualOrder, schemeParams{:}, opSet);
     % if ~isa(diffOp.D, 'function_handle')
@@ -16,7 +16,7 @@
 
     % DiffOp for flux in residual viscosity. Due to sign conventions of the implemented schemes, we need to
     % change the sign.
-    D_flux = @(v) -D_flux(v);
+    D_flux = @(v) -D_res(v);
     % DiffOp for time derivative in residual viscosity
     DvDt = D;
 end