From ac950fb2ae0515cdbe390418c426f806c05b68ea Mon Sep 17 00:00:00 2001 From: Elias Pipping <elias.pipping@fu-berlin.de> Date: Tue, 28 Feb 2012 14:31:47 +0100 Subject: [PATCH] Do not implement differential twice --- dune/tectonic/nicefunction.hh | 7 +------ 1 file changed, 1 insertion(+), 6 deletions(-) diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh index 83ae7def..626c59e6 100644 --- a/dune/tectonic/nicefunction.hh +++ b/dune/tectonic/nicefunction.hh @@ -81,12 +81,7 @@ class RuinaFunction : public NiceFunction { /* see above */ double virtual rightDifferential(double s) const { - double const arg = eta * s / h; - if (arg <= rho) - return 0; - - return coefficient * normalStress * - (a * std::log(arg) + mu + compound_state); + return leftDifferential(s); } /* -- GitLab