diff --git a/dune/gfe/localprojectedfefunction.hh b/dune/gfe/localprojectedfefunction.hh
index 846264b78c9319a420f1041d6be6e939ac6c5ba3..20157a4163a4ea8c58704b8c79f0d4a8df1585cb 100644
--- a/dune/gfe/localprojectedfefunction.hh
+++ b/dune/gfe/localprojectedfefunction.hh
@@ -193,17 +193,6 @@ Dune::FieldMatrix< K, m, p > operator* ( const Dune::FieldMatrix< K, m, n > &A,
 
       static const int spaceDim = TargetSpace::TangentVector::dimension;
 
-      static FieldMatrix<field_type,3,3> transpose(const FieldMatrix<field_type,3,3>& matrix)
-      {
-        FieldMatrix<field_type,3,3> result;
-
-        for (int i=0; i<3; i++)
-          for (int j=0; j<3; j++)
-            result[i][j] = matrix[j][i];
-
-        return result;
-      }
-
       static FieldMatrix<field_type,3,3> polarFactor(const FieldMatrix<field_type,3,3>& matrix)
       {
         // Use Higham's method