Mercurial > repos > public > sbplib

view +rv/+diffops/constructSchemeDiffOp.m @ 1184:ecc605453733 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Move specialized constructDiffOps* into constructDiffOps.m. Add multigrid RV using the standard RungekuttaRV time stepper
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 05 Jul 2019 17:51:11 +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