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

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, diffOp] = constructFluxDiffOp(scheme, g, order, schemeParams, opSet, BCs)
    diffOp = scheme(g, order, schemeParams{:}, opSet);
    if ~isa(diffOp.D, 'function_handle')
        D = @(v) diffOp.D*v;
    else
        D = diffOp.D;
    end
end