diff --git a/dune/tectonic/data-structures/friction/localfriction.hh b/dune/tectonic/data-structures/friction/localfriction.hh
index bbd78ce36f7c84172fbb4f817a28db5c6ca9f16c..6ef937e024b62325fa052954f835cfcbb034c7b8 100644
--- a/dune/tectonic/data-structures/friction/localfriction.hh
+++ b/dune/tectonic/data-structures/friction/localfriction.hh
@@ -62,7 +62,7 @@ class WrappedScalarFriction : public LocalFriction<dimension> {
   double regularity(VectorType const &x) const override {
     auto tangential_x = removeNormal(x);
     double const xnorm = tangential_x.two_norm();
-    if (xnorm <= 0.0)
+    if (xnorm < 0.0)
       return std::numeric_limits<double>::infinity();
 
     return func_.regularity(xnorm);