diff --git a/dune/gfe/mixedriemanniantrsolver.hh b/dune/gfe/mixedriemanniantrsolver.hh
index bb2ecf4b3d283acfbf0c03d1726cb812add15dec..7e5acd784df9bfdc4b810cc0b3ec55ab2fb83209 100644
--- a/dune/gfe/mixedriemanniantrsolver.hh
+++ b/dune/gfe/mixedriemanniantrsolver.hh
@@ -8,6 +8,8 @@
 #include <dune/istl/bcrsmatrix.hh>
 #include <dune/istl/bvector.hh>
 
+#include <dune/grid/utility/globalindexset.hh>
+
 #include <dune/solvers/common/boxconstraint.hh>
 #include <dune/solvers/norms/h1seminorm.hh>
 #include <dune/solvers/solvers/iterativesolver.hh>
@@ -19,7 +21,6 @@
 #include <dune/fufem/functionspacebases/p3nodalbasis.hh>
 
 #include <dune/gfe/mixedgfeassembler.hh>
-#include <dune/gfe/parallel/globalindex.hh>
 
 /** \brief Riemannian trust-region solver for geodesic finite-element problems */
 template <class GridType,