diff --git a/dune/gfe/localgfetestfunction.hh b/dune/gfe/localgfetestfunction.hh
index b87d031d60b0bf74e74daf5a09a53bb4dfb7d7f4..15ae8617922bbb79a5fdf54a0ad9904eaaee26e6 100644
--- a/dune/gfe/localgfetestfunction.hh
+++ b/dune/gfe/localgfetestfunction.hh
@@ -48,6 +48,12 @@ public:
     /** \brief Evaluate the derivatives of all shape functions function */
     void evaluateJacobian(const Dune::FieldVector<ctype, dim>& local,
                           std::vector<Dune::FieldMatrix<ctype, EmbeddedTangentVector::dimension, dim> >& out) const;
+                          
+    /** \brief Polynomial order */
+    unsigned int order() const
+    {
+        return localGFEFunction_.localFiniteElement_.order();
+    }
 
 private: