Mercurial > repos > public > sbplib
diff +rv/+diffops/constructDiffOpsBdf.m @ 1167:66c0fbbc406f feature/rv
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 |
line wrap: on
line diff
--- a/+rv/+diffops/constructDiffOpsBdf.m Thu Jun 27 17:15:11 2019 +0200 +++ b/+rv/+diffops/constructDiffOpsBdf.m Fri Jun 28 12:53:47 2019 +0200 @@ -1,7 +1,10 @@ -function [D_scheme, D_flux, penalties_scheme] = constructDiffOpsBdf(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs) +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) + [D_scheme, penalties_scheme, ~] = rv.diffops.constructSchemeDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs); %% DiffOps for residual viscosity - [D_flux, ~] = rv.diffops.constructFluxDiffOp(scheme, g, order, schemeParams, opSet, BCs) + [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 \ No newline at end of file