make sure we do not solve with a null right hand side

This commit is contained in:
Gael Guennebaud 2012-06-06 17:11:50 +02:00
parent fd32697074
commit 05af70a958

View File

@ -200,7 +200,7 @@ template<typename Solver> void check_sparse_spd_solving(Solver& solver)
SpMat B(size,rhsCols);
DenseVector b = DenseVector::Random(size);
DenseMatrix dB(size,rhsCols);
initSparse<Scalar>(density, dB, B);
initSparse<Scalar>(density, dB, B, ForceNonZeroDiag);
for (int i = 0; i < g_repeat; i++) {
check_sparse_solving(solver, A, b, dA, b);