This example fails with an assertion in the devel branch.

Code: Select all

bool test(int m=100, int n=10)
{
   MatrixXd A = MatrixXd::Random(m, n);
   VectorXd b(m), x(n);
   return A.svd().solve(b, &x);
}

Is there any workaround for this?

Thanks for the report, i committed a tentative fix, can you try (in the default branch)?

Works like a charm

Thanks for the feedback. Hint: do check that the result is sane (A*x approximately equal to b) because obviously our unit-test here is lacking, otherwise this bug wouldn't have existed in the first place.