Mercurial > repos > public > sbplib

view +rv/+diffops/constructSchemeDiffOp.m @ 1188:7173b6fd4063 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Rename class property grid to g, to avoid confusion with grid package
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 08 Jul 2019 14:50:19 +0200
parents 65a577db5ca0
children
line wrap: on
line source

function [D_scheme, penalties_scheme, D] = constructSchemeDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs)
    %% DiffOps for solution vector
    [D, penalties_scheme] = rv.diffops.constructFluxDiffOpWithClosures(scheme, g, schemeOrder, schemeParams, opSet, BCs);
    D2 = rv.diffops.constructSymmetricD2(g, schemeOrder, opSet);
    D_scheme = @(v,Viscosity)(D(v) + D2(Viscosity)*v);
end