solvers: assert prim/dual res or duality not nan

Author: fabinsch

include/proxsuite/proxqp/dense/solver.hpp

@@ -1317,6 +1317,11 @@ qp_solve( //
               << std::endl;
   }
   qpwork.dirty = true;
+
+  assert(!std::isnan(qpresults.info.pri_res));
+  assert(!std::isnan(qpresults.info.dua_res));
+  assert(!std::isnan(qpresults.info.duality_gap));
+
   PROXSUITE_EIGEN_MALLOC_ALLOWED();
 }
 

include/proxsuite/proxqp/sparse/solver.hpp

@@ -7,6 +7,8 @@
 #define PROXSUITE_PROXQP_SPARSE_SOLVER_HPP
 
 #include <chrono>
+#include <cmath>
+
 #include <proxsuite/linalg/dense/core.hpp>
 #include <proxsuite/linalg/sparse/core.hpp>
 #include <proxsuite/linalg/sparse/factorize.hpp>
@@ -1409,6 +1411,10 @@ qp_solve(Results<T>& results,
               << std::endl;
   }
 
+  assert(!std::isnan(results.info.pri_res));
+  assert(!std::isnan(results.info.dua_res));
+  assert(!std::isnan(results.info.duality_gap));
+
   work.set_dirty();
 }
 } // namespace sparse