diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 7ef0bc72e2727bd7c43a3affec5b692667ac6890..2ab1da2be6666d31eb60a48fc701a04c077db62b 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -331,7 +331,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
             CorrectionType gradient = rhs_global;
             for (size_t j=0; j<gradient.size(); j++)
               for (int k=0; k<gradient[j].size(); k++)
-                if ((*ignoreNodes_)[j][k])
+                if ((*mgStep->ignoreNodes_)[j][k])  // global Dirichlet nodes set
                   gradient[j][k] = 0;
 
             if (this->verbosity_ == Solver::FULL and rank==0)