From 044d511a57a08e3148f1d236788731034219cadd Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sun, 5 Jan 2014 20:18:32 +0100
Subject: [PATCH] [Extend] Handle the nonsmooth case x=0
---
dune/tectonic/localfriction.hh | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/dune/tectonic/localfriction.hh b/dune/tectonic/localfriction.hh
index c6c8900a..baac0904 100644
--- a/dune/tectonic/localfriction.hh
+++ b/dune/tectonic/localfriction.hh
@@ -57,8 +57,7 @@ template <size_t dimension> class LocalFriction {
Interval<double> &D) const {
double const xnorm = x.two_norm();
if (xnorm <= 0.0)
- D[0] = D[1] = 0.0; // WARNING: we assume that the outer
- // function approaches zero superlinearly
+ D[0] = D[1] = func->differential(0.0) * v.two_norm();
else
D[0] = D[1] = func->differential(xnorm) * dotFirstNormalised(x, v);
}
--
GitLab