diff --git a/dune/tectonic/samplefunctional.hh b/dune/tectonic/samplefunctional.hh index d2779321c26a22d62af3efcc3722262109745f4b..a56354d2bdae7ea49fab0c616f6716f73811ab4b 100644 --- a/dune/tectonic/samplefunctional.hh +++ b/dune/tectonic/samplefunctional.hh @@ -29,10 +29,10 @@ template <int dim> class SampleFunctional { double operator()(SmallVector const &v) const { SmallVector y; - A.mv(v, y); // Av - y /= 2; // 1/2 Av - y -= b; // 1/2 Av - b - return y * v + *phi(v); // <1/2 Av - b,v> + H(|v|) + A.mv(v, y); // Av + y /= 2; // 1/2 Av + y -= b; // 1/2 Av - b + return y * v + (*phi)(v); // <1/2 Av - b,v> + H(|v|) } bool descentDirection(SmallVector const x, SmallVector &ret) const {