diff --git a/dune/solvers/iterationsteps/mmgstep.cc b/dune/solvers/iterationsteps/mmgstep.cc index 72deb29fcf3edc9f473989c24dc0cd040d92b7d6..470f5a5987a7ebf33131e05f4693271b0dc271c6 100644 --- a/dune/solvers/iterationsteps/mmgstep.cc +++ b/dune/solvers/iterationsteps/mmgstep.cc @@ -14,6 +14,7 @@ template <class MatrixType, class VectorType> void MonotoneMGStep<MatrixType, VectorType>:: preprocess() { + // Unset the recompute bitfields, so we compute the full stiffness matrix hierarchy at the beginning for (size_t i=0; i<this->mgTransfer_.size(); i++) dynamic_cast<TruncatedMGTransfer<VectorType>*>(this->mgTransfer_[i])->setRecomputeBitField(nullptr);