Mercurial > repos > public > sbplib

comparison +rv/+diffops/constructDiffOpsBdf.m @ 1167:66c0fbbc406f feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make constructDiffOps return a struct with the operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 12:53:47 +0200
parents 65a577db5ca0
children
comparison
equal deleted inserted replaced
1166:17bb3c46472c 1167:66c0fbbc406f
1 function [D_scheme, D_flux, penalties_scheme] = constructDiffOpsBdf(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs) 1 function diffOpsStruct = constructDiffOpsBdf(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs)
2 %% DiffOps for solution vector 2 %% DiffOps for solution vector
3 [D_scheme, penalties_scheme, ~] = rv.diffops.constructSchemeDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs) 3 [D_scheme, penalties_scheme, ~] = rv.diffops.constructSchemeDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs);
4 %% DiffOps for residual viscosity 4 %% DiffOps for residual viscosity
5 [D_flux, ~] = rv.diffops.constructFluxDiffOp(scheme, g, order, schemeParams, opSet, BCs) 5 [D_flux, ~] = rv.diffops.constructFluxDiffOp(scheme, g, residualOrder, schemeParams, opSet, BCs);
6 D_flux = @(v) -D_flux(v); 6 D_flux = @(v) -D_flux(v);
7 diffOpsStruct = struct('D_scheme', D_scheme,...
8 'D_flux', D_flux,...
9 'penalties_scheme', penalties_scheme);
7 end 10 end