diff --git a/dune/tectonic/ellipticenergy.hh b/dune/tectonic/ellipticenergy.hh index 8359136afe258d1bb6ce1cd928edb89d39c1beb2..c30626f733dfb37cb51cd9f8181c98f1be682586 100644 --- a/dune/tectonic/ellipticenergy.hh +++ b/dune/tectonic/ellipticenergy.hh @@ -294,7 +294,8 @@ void minimisationInitialiser(Functional const &J, Bisection const &bisection, tangentialMinimisation(J, x1, descDir1, bisection); double const Jx1 = J(x1); - SmallVector x2 = { -x1[0], -x1[1] }; + SmallVector x2(0); + x2.axpy(-1, x1); SmallVector const descDir2 = { x2[1], -x2[0] }; tangentialMinimisation(J, x2, descDir2, bisection); double const Jx2 = J(x2);