diff --git a/dune/tectonic/samplefunctional.hh b/dune/tectonic/samplefunctional.hh index 50d8f7850650f4bca3bc4eb92592f8532a333bd3..6ca943e90d19268dad6516ebc8aa731c9bb831ec 100644 --- a/dune/tectonic/samplefunctional.hh +++ b/dune/tectonic/samplefunctional.hh @@ -29,9 +29,9 @@ template <int dim> class SampleFunctional { double operator()(SmallVector const &v) const { SmallVector y; - A.mv(v, y); // y = Av - y /= 2; // y = 1/2 Av - y -= b; // y = 1/2 Av - b + 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|) }