MyConvexProblem::operator()

dune/tectonic/myconvexproblem.hh

@@ -27,6 +27,14 @@ class MyConvexProblem {
                   VectorType const &f)
       : A(A), phi(phi), f(f) {};
 
+  /* Just for debugging */
+  double operator()(VectorType const &x) const {
+    double ret = phi(x) - (f * x);
+    VectorType tmp(x.size());
+    A.mv(x, tmp);
+    return ret + 0.5 * (tmp * x);
+  }
+
   MatrixType const &A;
   Dune::GlobalNonlinearity<MatrixType, VectorType> const &phi;