fix Model::fix did not take into account quadratic terms

hlefebvr · Dec 18, 2023 · 4b345ff · 4b345ff
1 parent 694cc7a
commit 4b345ff
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/lib/src/modeling/matrix/Row.cpp b/lib/src/modeling/matrix/Row.cpp
@@ -16,6 +16,10 @@ idol::Row idol::Row::fix(const Solution::Primal &t_primals) const {
         result.linear() += constant.fix(t_primals) * var;
     }
 
+    for (const auto& [var1, var2, constant] : quadratic()) {
+        result.quadratic() += constant.fix(t_primals) * var1 * var2;
+    }
+
     result.rhs() = rhs().fix(t_primals);
 
     return result;