|
@@ -247,7 +247,9 @@ IGL_INLINE bool igl::mosek_quadprog(
|
|
// Q should be square
|
|
// Q should be square
|
|
assert(Q.rows() == Q.cols());
|
|
assert(Q.rows() == Q.cols());
|
|
// Q should be symmetric
|
|
// Q should be symmetric
|
|
|
|
+#ifdef EIGEN_HAS_A_BUG_AND_FAILS_TO_LET_ME_COMPUTE_Q_MINUS_Q_TRANSPOSE
|
|
assert( (Q-Q.transpose()).sum() < FLOAT_EPS);
|
|
assert( (Q-Q.transpose()).sum() < FLOAT_EPS);
|
|
|
|
+#endif
|
|
// Only keep lower triangular part of Q
|
|
// Only keep lower triangular part of Q
|
|
SparseMatrix<double> QL;
|
|
SparseMatrix<double> QL;
|
|
//QL = Q.template triangularView<Lower>();
|
|
//QL = Q.template triangularView<Lower>();
|