Skip to content
Snippets Groups Projects
Commit 31b3f152 authored by graeser's avatar graeser Committed by graeser
Browse files

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]]
parent 239ae461
No related branches found
No related tags found
No related merge requests found
......@@ -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]);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment