From 7a5920c39d1d6913755afbf6c2f9ffe688f17616 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Thu, 16 Aug 2012 09:14:57 +0200
Subject: [PATCH] RuinaFunction::evaluate()
---
dune/tectonic/nicefunction.hh | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index 7acf4bec..10b7417d 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -51,6 +51,15 @@ class RuinaFunction : public NiceFunction {
// logarithmic
rho(std::exp(-K / a)) {}
+ double virtual evaluate(double x) const {
+ double const arg = x / h * eta;
+ if (arg <= rho)
+ return 0;
+
+ return coefficientProduct *
+ (+a * arg * (std::log(arg) - 1) + K * arg + a * rho);
+ }
+
double virtual leftDifferential(double x) const {
double const arg = x / h * eta;
if (arg <= rho)
--
GitLab