From 403aa801b12fcc9a347873c3e3233fceac066b20 Mon Sep 17 00:00:00 2001
From: podlesny <podlesny@zedat.fu-berlin.de>
Date: Mon, 1 Mar 2021 13:00:52 +0100
Subject: [PATCH] fix issue for no friction case
---
dune/tectonic/data-structures/friction/localfriction.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/tectonic/data-structures/friction/localfriction.hh b/dune/tectonic/data-structures/friction/localfriction.hh
index bbd78ce3..6ef937e0 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);
--
GitLab