From 31b3f15253f3f34a359dde04c17c9c03ada14eb2 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Carsten=20Gr=C3=A4ser?= <graeser@mi.fu-berlin.de>
Date: Tue, 15 Jan 2013 12:52:54 +0000
Subject: [PATCH] Set size of coarseSolution[maxLevel-1]
Otherwise the vector is not initialized correctly
if we only have two levels leading to an invalid
read.
[[Imported from SVN: r7676]]
---
dune/solvers/iterationsteps/obstacletnnmgstep.hh | 1 +
1 file changed, 1 insertion(+)
diff --git a/dune/solvers/iterationsteps/obstacletnnmgstep.hh b/dune/solvers/iterationsteps/obstacletnnmgstep.hh
index 2b36ca51..dec535ca 100644
--- a/dune/solvers/iterationsteps/obstacletnnmgstep.hh
+++ b/dune/solvers/iterationsteps/obstacletnnmgstep.hh
@@ -376,6 +376,7 @@ class ObstacleTNNMGStep
transfer_[maxLevel-1]->restrict(rhs_, coarseRhs[maxLevel-1]);
transfer_[maxLevel-1]->restrictToFathers(*ignoreNodes_, coarseIgnore[maxLevel-1]);
obstacleRestrictor.restrict(obstacles_, coarseObstacle[maxLevel-1], hasObstacle_, hasObstacle_, *(transfer_[maxLevel-1]), critical);
+ coarseSolution[maxLevel-1].resize(coarseMatrix[maxLevel-1].N());
for (int i = maxLevel-2; i>=0; --i)
{
transfer_[i]->galerkinRestrictSetOccupation(coarseMatrix[i+1], coarseMatrix[i]);
--
GitLab