diff --git a/dune/gfe/functions/localgeodesicfefunction.hh b/dune/gfe/functions/localgeodesicfefunction.hh
index ea5342c160bf09e7d6547ce00c10fbf1a1c33407..e8b29687e7c1df8d43b5c270d32811d335a7df6b 100644
--- a/dune/gfe/functions/localgeodesicfefunction.hh
+++ b/dune/gfe/functions/localgeodesicfefunction.hh
@@ -251,17 +251,6 @@ namespace Dune::GFE
   {
     Dune::FieldMatrix<RT, embeddedDim, dim> result;
 
-#if 0  // this is probably faster than the general implementation, but we leave it out for testing purposes
-    if (dim==1) {
-
-      EmbeddedTangentVector tmp = TargetSpace::interpolateDerivative(coefficients_[0], coefficients_[1], local[0]);
-
-      for (int i=0; i<embeddedDim; i++)
-        result[i][0] = tmp[i];
-
-    }
-#endif
-
     // ////////////////////////////////////////////////////////////////////////
     //  The derivative is evaluated using the implicit function theorem.
     //  Hence we need to solve a small system of linear equations.