Kill member mu

dune/tectonic/nicefunction.hh

@@ -37,20 +37,19 @@ class RuinaFunction : public NiceFunction {
   RuinaFunction(double coefficient, double a, double mu, double eta,
                 double normalStress, double b, double state, double L, double h)
       : a(a),
-        mu(mu),
         eta(eta),
         h(h),
         coefficientProduct(coefficient * normalStress),
-        K(b *
-          (state - std::log(eta * L))), // state is assumed to be logarithmic
-        rho(std::exp(-(mu + K) / a)) {}
+        K(mu + b * (state - std::log(eta * L))), // state is assumed to be
+                                                 // logarithmic
+        rho(std::exp(-K / a)) {}
 
   double virtual leftDifferential(double x) const {
     double const arg = x / h * eta;
     if (arg <= rho)
       return 0;
 
-    double const ret = a * std::log(arg) + mu + K;
+    double const ret = a * std::log(arg) + K;
     return coefficientProduct * ret;
   }
 
@@ -77,7 +76,6 @@ class RuinaFunction : public NiceFunction {
 
 private:
   double const a;
-  double const mu;
   double const eta;
   double const h;
   double const coefficientProduct;