Mercurial > repos > public > sbplib

comparison +rv/+diffops/constructDiffOpsMultiGrid.m @ 1171:393ed30866a1 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Actually return coarse grid operators from constructDiffOpsMultiGrid
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 13:18:58 +0200
parents af3c4eb0cbbd
children 9dd10dcff128
comparison
equal deleted inserted replaced
1170:d3bde8a23e08 1171:393ed30866a1
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.constructFluxDiffOpWithClosures(scheme, g, schemeOrder, schemeParams, opSet, BCs); 5 D_coarse = rv.diffops.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,...