Clearer

dune/tectonic/ellipticenergy.hh

@@ -39,6 +39,7 @@ template <int dim> class EllipticEnergy {
       // If there is a direction of descent, this is it
       SmallVector d;
       smoothGradient(x, d);
+      d *= -1;
 
       Interval<double> D;
       phi->directionalSubDiff(x, d, D);
@@ -48,7 +49,6 @@ template <int dim> class EllipticEnergy {
 
       if (combinedDecline < 0) {
         ret = d;
-        ret *= -1;
       } else {
         ret = 0;
       }