diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh index 5ff1cc755a874da551e4753f9aef92d297c0bfcc..8fd4671a8a5b3a894dc63be9b0dbb2fc3b702e63 100644 --- a/dune/tectonic/myblockproblem.hh +++ b/dune/tectonic/myblockproblem.hh @@ -79,8 +79,11 @@ class MyBlockProblem<MyConvexProblemTypeTEMPLATE>::IterateObject { void solveLocalProblem( LocalVectorType& ui, int m, const typename Dune::BitSetVector<block_size>::const_reference ignore) { - // Note: ignore is currently ignored (what's it used for anyway?) { + // TODO: Does it make any sense to ignore single spatial dimensions here? + if (ignore.test(0)) + return; + LocalMatrixType const* localA = NULL; LocalVectorType localb(problem.f[m]);