diff --git a/dune/tectonic/circularconvexfunction.hh b/dune/tectonic/circularconvexfunction.hh index ab72140568fa50da130ca7e206ed6185a8c3cf03..65f79b5d65971fdfa84e868d3975a8ea7657236f 100644 --- a/dune/tectonic/circularconvexfunction.hh +++ b/dune/tectonic/circularconvexfunction.hh @@ -27,10 +27,9 @@ template <class MatrixType, class VectorType> class CircularConvexFunction { tangent(m, t); VectorType tmp; - A.mv(x, tmp); // Ax - tmp -= b; // Ax - b - double const dotp = tmp * t; // <Ax - b,t> - D[0] = D[1] = dotp; + A.mv(x, tmp); // Ax + tmp -= b; // Ax - b + D[0] = D[1] = tmp * t; // <Ax - b,t> } void domain(Interval<double> &domain) const {