From 510440f1b97b5ceaaccf5225042689877aa0e6fe Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 29 Jun 2012 13:01:01 +0200
Subject: [PATCH] RuinaFunction(|.|) is smooth at zero!
---
dune/tectonic/nicefunction.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index 07f6e606..66808943 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -68,12 +68,14 @@ class RuinaFunction : public NiceFunction {
double virtual regularity(double x) const {
double const arg = x / h * eta;
// TODO: Make this controllable
- if (arg < 1e-14 || std::abs(arg - rho) < 1e-14)
+ if (std::abs(arg - rho) < 1e-14)
return std::numeric_limits<double>::infinity();
return std::abs(second_deriv(x));
}
+ bool virtual smoothesNorm() const { return true; }
+
private:
double const a;
double const eta;
--
GitLab