diff --git a/dune/gfe/globalgeodesicfefunction.hh b/dune/gfe/globalgeodesicfefunction.hh
index 914c376dd86d3845c414bd26898f18bb5d2da5a4..26603090b5f2db9d408f7d1b2dd44c09d168ce21 100644
--- a/dune/gfe/globalgeodesicfefunction.hh
+++ b/dune/gfe/globalgeodesicfefunction.hh
@@ -8,6 +8,7 @@
 
 #include <dune/gfe/localgeodesicfefunction.hh>
 
+#include <dune/fufem/functions/virtualgridfunction.hh>
 
 /** \brief Global geodesic finite element function. 
  *
@@ -16,7 +17,7 @@
  *  \tparam CoefficientType - The coefficient vector type.
  */
 template<class Basis, class TargetSpace, class CoefficientType>
-GlobalGeodesicFEFunction {
+GlobalGeodesicFEFunction : public VirtualGridFunction<typename Basis::GridView::Grid, TargetSpace> {
 
 
     typedef typename Basis::LocalFiniteElement LocalFiniteElement;
diff --git a/dune/gfe/globalgfetestfunction.hh b/dune/gfe/globalgfetestfunction.hh
index f026b6f5629e82a9350e891af5ef9cc50504dfad..453b97379057e503f7ddbd86221e7b01d7537b84 100644
--- a/dune/gfe/globalgfetestfunction.hh
+++ b/dune/gfe/globalgfetestfunction.hh
@@ -10,6 +10,7 @@
 
 #include <dune/gfe/localgfetestfunction.hh>
 
+#include <dune/fufem/functions/virtualgridfunction.hh>
 
 /** \brief Global geodesic finite element test function. 
  *
@@ -18,7 +19,7 @@
  *  \tparam CoefficientType - The coefficient vector type.
  */
 template<class Basis, class TargetSpace, class CoefficientType>
-GlobalGFETestFunction {
+GlobalGFETestFunction : public VirtualGridFunction<typename Basis::GridView::Grid, typename TargetSpace::EmbeddedTangentVector>  {
 
     typedef typename Basis::LocalFiniteElement LocalFiniteElement;
     typedef typename Basis::GridView GridView;