diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh index 07f6e606f70fba3c89a12b54559bd56169c8abf8..66808943bc30ed429fd8c42564f8e07544c13a25 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;