From a288b5b4139acc33e59d99182021c6e445981dd9 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Mon, 27 Feb 2012 15:15:21 +0100
Subject: [PATCH] Fix RuinaFunction::evaluate()
---
dune/tectonic/nicefunction.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index fe2008d9..83ae7def 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -58,7 +58,7 @@ class RuinaFunction : public NiceFunction {
void virtual evaluate(double const &x, double &y) const {
double const arg = std::max(eta * x / h, rho);
double const r = arg * (a * (std::log(arg) - 1) + mu + compound_state);
- y = 1 / eta * normalStress * r + a * rho + mu + compound_state;
+ y = 1 / eta * normalStress * (r + a * rho + mu + compound_state);
y *= coefficient;
y *= h;
}
--
GitLab