From 315471dcc62f8ba6e9f189b7b38c6415ebc67cf6 Mon Sep 17 00:00:00 2001
From: podlesny <podlesny@zedat.fu-berlin.de>
Date: Sat, 6 Mar 2021 13:27:06 +0100
Subject: [PATCH] add regularity check for hessian
---
dune/tectonic/data-structures/friction/localfriction.hh | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/dune/tectonic/data-structures/friction/localfriction.hh b/dune/tectonic/data-structures/friction/localfriction.hh
index 6ef937e0..7a477b5c 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;
--
GitLab