Mercurial > repos > public > sbplib

comparison +rv/constructDiffOps.m @ 1023:defc9d0cc1f2 feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove incorrect assertion of the number of BC:s
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 07 Jan 2019 12:06:49 +0100
parents cc61dde120cd
children 44c3ea38097e
comparison
equal deleted inserted replaced
1022:234c1c02ea39 1023:defc9d0cc1f2
1 function [D_rv, D_flux, DvDt, solutionPenalties, residualPenalties] = constructDiffOps(scheme, grid, order, opSet, waveSpeed, BCs, fluxSplitting) 1 function [D_rv, D_flux, DvDt, solutionPenalties, residualPenalties] = constructDiffOps(scheme, grid, order, opSet, waveSpeed, BCs, fluxSplitting)
2 default_arg('fluxSplitting',[]); 2 default_arg('fluxSplitting',[]);
3 assert(size(BCs,1) == grid.D());
4 3
5 %% DiffOps for solution vector 4 %% DiffOps for solution vector
6 [D, solutionPenalties] = constructTotalFluxDiffOp(scheme, grid, order, opSet, waveSpeed, BCs, fluxSplitting); 5 [D, solutionPenalties] = constructTotalFluxDiffOp(scheme, grid, order, opSet, waveSpeed, BCs, fluxSplitting);
7 D2 = constructSymmetricD2Operator(grid, order, opSet); 6 D2 = constructSymmetricD2Operator(grid, order, opSet);
8 D_rv = @(v,viscosity)(D + D2(viscosity))*v; 7 D_rv = @(v,viscosity)(D + D2(viscosity))*v;