diff --git a/contrib/HighOrderMeshOptimizer/OptHomObjContribMetricMin.h b/contrib/HighOrderMeshOptimizer/OptHomObjContribMetricMin.h
index b3b7bb3975a75455b8b876fbe986a398048fd271..251fbe6c9e9b72e9c834b31c6d400eec59fbe123 100644
--- a/contrib/HighOrderMeshOptimizer/OptHomObjContribMetricMin.h
+++ b/contrib/HighOrderMeshOptimizer/OptHomObjContribMetricMin.h
@@ -65,9 +65,9 @@ bool ObjContribMetricMin<FuncType>::addContrib(double &Obj, alglib::real_1d_arra
     _mesh->metricMinAndGradients(iEl, mM, gMM);
     for (int l = 0; l < _mesh->nBezEl(iEl); l++) {                                  // Add contribution for each Bezier coeff.
       Obj += _weight * FuncType::compute(mM[l]);
+      const double dfact = _weight * FuncType::computeDiff(mM[l]);
       for (int iPC = 0; iPC < _mesh->nPCEl(iEl); iPC++)
-        gradObj[_mesh->indPCEl(iEl,iPC)] +=
-            _weight * FuncType::computeDiff(mM[l], gMM[_mesh->indGSJ(iEl, l, iPC)]);
+        gradObj[_mesh->indPCEl(iEl,iPC)] += dfact *  gMM[_mesh->indGSJ(iEl, l, iPC)];
       _min = std::min(_min, mM[l]);
       _max = std::max(_max, mM[l]);
     }