diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index d8439808fe3bc02f8d78fb2db2a8ba2361b323aa..fe4261b7f2fd5cee4f3e439ba1073abcb51b3295 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -44,16 +44,6 @@ Dune::FieldMatrix<K,m,n> operator- ( const Dune::FieldMatrix<K, m, n> &A, const
     return ret;
 }
 
-#if 0
-template< class K, int m, int n>
-void transpose(Dune::FieldMatrix<K, m, n> &A)
-{
-   for( size_type i = 0; i < m; ++i )
-        for( size_type j = 0; j < i; ++j )
-			std::swap(A[i][j], A[j][i]);
-}
-#endif
-
 
 /** \brief A function defined by simplicial geodesic interpolation 
            from the reference element to a Riemannian manifold.