Mercurial > repos > public > sbplib

view +rv/+diffops/constructDiffOpsBdf.m @ 1178:e3d8f24b2c1c feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add wrapper function calling constructDiffOps* for a specified RV method
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 15:13:44 +0200
parents 66c0fbbc406f
children
line wrap: on
line source

function diffOpsStruct = constructDiffOpsBdf(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs)
    %% DiffOps for solution vector
    [D_scheme, penalties_scheme, ~] = rv.diffops.constructSchemeDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs);
    %% DiffOps for residual viscosity
    [D_flux, ~] = rv.diffops.constructFluxDiffOp(scheme, g, residualOrder, schemeParams, opSet, BCs);
    D_flux = @(v) -D_flux(v);
    diffOpsStruct = struct('D_scheme', D_scheme,...
                 'D_flux', D_flux,...
                 'penalties_scheme', penalties_scheme);
end