diff --git a/dune/tectonic/localnonlinearity.hh b/dune/tectonic/localnonlinearity.hh
index 339959ef1ff81153460b756e502503e329ebc1aa..fcf0d7b1252427a540bab1f03f09207663787566 100644
--- a/dune/tectonic/localnonlinearity.hh
+++ b/dune/tectonic/localnonlinearity.hh
@@ -27,7 +27,9 @@ template <int dimension> class LocalNonlinearity {
   }
 
   double regularity(VectorType const &x) const {
-    if (!func_->smoothesNorm() && x.two_norm() < 1e-12) // TODO
+    if (!func_->smoothesNorm() && x.two_norm() < 1e-14) // TODO: Make this
+                                                        // controllable
+                                                        // (truncationRadius?)
       return std::numeric_limits<double>::infinity();
 
     return func_->regularity(x.two_norm());