diff --git a/dune/tectonic/data-structures/friction/localfriction.hh b/dune/tectonic/data-structures/friction/localfriction.hh index 6ef937e024b62325fa052954f835cfcbb034c7b8..7a477b5c10adeebaa19526d1c9739d1766e0b475 100644 --- a/dune/tectonic/data-structures/friction/localfriction.hh +++ b/dune/tectonic/data-structures/friction/localfriction.hh @@ -136,7 +136,8 @@ class WrappedScalarFriction : public LocalFriction<dimension> { }*/ auto tangential_x = removeNormal(x); const double xnorm = tangential_x.two_norm(); - if (xnorm <= 0.0) + + if (std::isinf(func_.regularity(xnorm)) or xnorm<=0.0) return; VectorType y = tangential_x;