diff --git a/dune/gfe/geodesicfefunctionadaptor.hh b/dune/gfe/geodesicfefunctionadaptor.hh
index 48a26f3f04da9daca1630946868a9648582602b4..cb192d0fa0f27c9d87bb9d639a4bd6a091b77252 100644
--- a/dune/gfe/geodesicfefunctionadaptor.hh
+++ b/dune/gfe/geodesicfefunctionadaptor.hh
@@ -55,7 +55,7 @@ void geodesicFEFunctionAdaptor(GridType& grid, std::vector<TargetSpace>& x)
     for (; eIt!=eEndIt; ++eIt) {
 
         // Set up a local gfe function on the father element
-        Dune::array<TargetSpace,dim+1> coefficients;
+        std::vector<TargetSpace> coefficients(dim+1);
 
         for (int i=0; i<eIt->father()->template count<dim>(); i++)
             coefficients[i] = dofMap.find(idSet.subId(*eIt->father(),i,dim))->second;