IPM: fixed computation of quadratic term in evaluate_subproblem_objec…

…tive()
cvanaret · Dec 12, 2024 · 035dd45 · 035dd45
1 parent d2ffd7c
commit 035dd45
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/...ents/inequality_handling_methods/interior_point_methods/PrimalDualInteriorPointMethod.cpp b/...ents/inequality_handling_methods/interior_point_methods/PrimalDualInteriorPointMethod.cpp
@@ -348,7 +348,7 @@ namespace uno {
 
    double PrimalDualInteriorPointMethod::evaluate_subproblem_objective(const Direction& direction) const {
       const double linear_term = dot(direction.primals, this->objective_gradient);
-      const double quadratic_term = 0.; // TODO this->sol->hessian.quadratic_product(direction.primals, direction.primals) / 2.;
+      const double quadratic_term = this->hessian.quadratic_product(direction.primals, direction.primals) / 2.;
       return linear_term + quadratic_term;
    }