diff --git a/dune/tectonic/localfriction.hh b/dune/tectonic/localfriction.hh
index c6c8900aaa5ea0b33f67e3f5f7e272bd71ecb86c..baac090430c2d67e28830c39267c64c072bfe120 100644
--- a/dune/tectonic/localfriction.hh
+++ b/dune/tectonic/localfriction.hh
@@ -57,8 +57,7 @@ template <size_t dimension> class LocalFriction {
                           Interval<double> &D) const {
     double const xnorm = x.two_norm();
     if (xnorm <= 0.0)
-      D[0] = D[1] = 0.0; // WARNING: we assume that the outer
-                         // function approaches zero superlinearly
+      D[0] = D[1] = func->differential(0.0) * v.two_norm();
     else
       D[0] = D[1] = func->differential(xnorm) * dotFirstNormalised(x, v);
   }