view +rv/+diffops/constructDiffOpsStandard.m @ 1171:393ed30866a1 feature/rv

Actually return coarse grid operators from constructDiffOpsMultiGrid
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 13:18:58 +0200
parents 65a577db5ca0
children
line wrap: on
line source

function [D_scheme, D_flux, D_t, penalties_scheme, penalties_res] = constructDiffOpsStandard(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs)
    % DiffOps for solution vector
    [D_scheme, D, penalties_scheme] = rv.diffops.constructSchemeDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs)
    %% DiffOps for residual viscosity
    [D_flux, penalties_res] = rv.diffops.constructFluxDiffOpWithClosures(scheme, g, residualOrder, schemeParams, opSet, BCs);
    D_flux = @(v) -D_flux(v);
    D_t = D;
end