Mercurial > repos > public > sbplib

comparison +rv/+diffops/constructDiffOpsMultiGrid.m @ 1177:270611e08398 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in constructDiffOpsMultiGrid
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 15:11:45 +0200
parents 9dd10dcff128
children
comparison
equal deleted inserted replaced
1176:ebec2b86f539 1177:270611e08398
1 function diffOpsStruct = constructDiffOpsMultiGrid(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs) 1 function diffOpsStruct = constructDiffOpsMultiGrid(scheme, g, schemeOrder, residualOrder, schemeParams, opSet, BCs)
2 % DiffOps for stabilized solution vector 2 % DiffOps for stabilized 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 % DiffOp for unstabilized solution vector 4 % DiffOp for unstabilized solution vector
5 D_coarse = rv.diffops.constructCoarserGridDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs); 5 D_coarse = constructCoarserGridDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs);
6 %% DiffOps for residual viscosity 6 %% DiffOps for residual viscosity
7 [D_t, penalties_res] = rv.diffops.constructFluxDiffOpWithClosures(scheme, g, residualOrder, schemeParams, opSet, BCs); 7 [D_t, penalties_res] = rv.diffops.constructFluxDiffOpWithClosures(scheme, g, residualOrder, schemeParams, opSet, BCs);
8 D_flux = @(v) -D_t(v); 8 D_flux = @(v) -D_t(v);
9 diffOpsStruct = struct('D_scheme', D_scheme,... 9 diffOpsStruct = struct('D_scheme', D_scheme,...
10 'D_coarse', D_coarse,... 10 'D_coarse', D_coarse,...