diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc index ea77c53deeeaf6900e80a876eac0d866d98c9e16..2f0063490136944e60043e9a93677bcd32d9e572 100644 --- a/dune/gfe/riemanniantrsolver.cc +++ b/dune/gfe/riemanniantrsolver.cc @@ -213,7 +213,8 @@ setup(const GridType& grid, if (rank==0) { hasObstacle_.resize(numLevels); - hasObstacle_.back().resize(dynamic_cast<TruncatedCompressedMGTransfer<CorrectionType>* >(mmgStep->mgTransfer_.back())->getMatrix().N(), true); + hasObstacle_.back().resize(basis.size(), true); + // If there is more than one level use the transfer operators to determine the correct coarse level system sizes for (int i=0; i<hasObstacle_.size()-1; i++) hasObstacle_[i].resize(dynamic_cast<TruncatedCompressedMGTransfer<CorrectionType>* >(mmgStep->mgTransfer_[i])->getMatrix().M(),true); }