RuinaFunction::evaluate()

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)