[Extend] Handle the nonsmooth case x=0

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);
   }