From b00c6092dd7c894b9bbb63e189bfa152e653484e Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 2 Nov 2011 01:53:53 +0100
Subject: [PATCH] Handle ignored nodes
---
dune/tectonic/myblockproblem.hh | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index 5ff1cc75..8fd4671a 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]);
--
GitLab