Mercurial > repos > public > sbplib
changeset 1172:9dd10dcff128 feature/rv
Dont return penalty terms from operators on coarse grid. One can use the penalty terms constructed for the fine grid
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Fri, 28 Jun 2019 13:22:03 +0200 |
parents | 393ed30866a1 |
children | dde29b865244 |
files | +rv/+diffops/constructDiffOpsMultiGrid.m |
diffstat | 1 files changed, 2 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
--- a/+rv/+diffops/constructDiffOpsMultiGrid.m Fri Jun 28 13:18:58 2019 +0200 +++ b/+rv/+diffops/constructDiffOpsMultiGrid.m Fri Jun 28 13:22:03 2019 +0200 @@ -14,14 +14,13 @@ 'penalties_res', penalties_res); end -% TODO: Not clear if this is needed % TODO: Only works for equidistant grids -function [D_coarse, penalties_coarse] = constructCoarserGridDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs) +function D_coarse = constructCoarserGridDiffOp(scheme, g, schemeOrder, schemeParams, opSet, BCs) m = g.m(); lim = g.lim(); m_coarse = (m-1)/2 + 1; g_coarse = grid.equidistant(m_coarse, lim{1}); - [D_coarse, penalties_coarse] = rv.diffops.constructFluxDiffOpWithClosures(scheme, g_coarse, schemeOrder, schemeParams, opSet, BCs); + D_coarse = rv.diffops.constructFluxDiffOpWithClosures(scheme, g_coarse, schemeOrder, schemeParams, opSet, BCs); x = g.x{1}; x_coarse = x(1:2:end); % TODO: Fix interpolation