Mercurial > repos > public > sbplib

view +rv/+diffops/constructDiffOpsBdf.m @ 1172:9dd10dcff128 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Dont return penalty terms from operators on coarse grid. One can use the penalty terms constructed for the fine grid
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 13:22:03 +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