Mercurial > repos > public > sbplib

diff +rv/+diffops/constructSchemeDiffOps.m @ 1190:79618b58b0a0 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Refactor constructDiffops and remove closures from some of the diffops used to compute the residual.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 08 Jul 2019 15:12:19 +0200
parents +rv/+diffops/constructSchemeDiffOp.m@65a577db5ca0
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/+rv/+diffops/constructSchemeDiffOps.m	Mon Jul 08 15:12:19 2019 +0200
@@ -0,0 +1,7 @@
+function [D_scheme, penalties, D_f] = constructSchemeDiffOps(scheme, g, order, schemeParams, opSet, BCs)
+    %% DiffOps for solution vector
+    [D_f, closures, penalties] = rv.diffops.constructFluxDiffOps(scheme, g, order, schemeParams, opSet, BCs);
+    D_scheme = rv.diffops.addClosuresToDiffOp(D_f, closures);
+    D2 = rv.diffops.constructSymmetricD2(g, order, opSet);
+    D_scheme = @(v,viscosity)(D_scheme(v) + D2(viscosity)*v);
+end
\ No newline at end of file