diff --git a/dune/gfe/globalgeodesicfefunction.hh b/dune/gfe/globalgeodesicfefunction.hh
index 26603090b5f2db9d408f7d1b2dd44c09d168ce21..a85e0de5e821d5e7aa60d7e4c1e37f418bb49152 100644
--- a/dune/gfe/globalgeodesicfefunction.hh
+++ b/dune/gfe/globalgeodesicfefunction.hh
@@ -12,13 +12,15 @@
 
 /** \brief Global geodesic finite element function. 
  *
- *  \tparam Basis  - The global basis type.
+ *  \tparam B  - The global basis type.
  *  \tparam TargetSpace - The manifold that this functions takes its values in.
  *  \tparam CoefficientType - The coefficient vector type.
  */
-template<class Basis, class TargetSpace, class CoefficientType>
+template<class B, class TargetSpace, class CoefficientType>
 GlobalGeodesicFEFunction : public VirtualGridFunction<typename Basis::GridView::Grid, TargetSpace> {
 
+public:
+    typedef B Basis;
 
     typedef typename Basis::LocalFiniteElement LocalFiniteElement;
     typedef typename Basis::GridView GridView;
@@ -34,7 +36,6 @@ GlobalGeodesicFEFunction : public VirtualGridFunction<typename Basis::GridView::
     //! Dimension of the embedded tanget space
     enum { embeddedDim = EmbeddedTangentVector::dimension };
 
-public:
 
     //! Create global function by a global basis and the corresponding coefficient vector
     GlobalGeodesicFEFunction(const Basis& basis, const CoefficientType& coefficients) :
@@ -77,6 +78,12 @@ public:
         out = localGFE.evaluateDerivative(local);
     }
 
+    /** \brief Export basis */
+    const Basis& basis() const
+    {
+        return basis_;
+    }
+
 private:
     //! The global basis
     const Basis& basis_;