From ab5bba608e3af386aec40b975e28cb9ef881c8ed Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Mon, 27 Feb 2012 13:49:27 +0100
Subject: [PATCH] Fix definition of rho

---
 dune/tectonic/nicefunction.hh | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index cc779fc5..fe2008d9 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -35,11 +35,11 @@ class RuinaFunction : public NiceFunction {
       : coefficient(coefficient),
         a(a),
         mu(mu),
-        rho(exp(-mu / a)),
         eta(eta),
         normalStress(normalStress),
         compound_state(b * (state - log(eta * L))),
-        h(h) {}
+        h(h),
+        rho(exp(-(mu + compound_state) / a)) {}
 
   /*
     If mu and sigma_n denote the coefficient of friction and the
-- 
GitLab