issue with matrix inversion about basiclinearalgebra HOT 4 CLOSED

Sounds good !

from basiclinearalgebra.

Hey @JM-FRANCE, thanks for the info, you're quite right. I've pushed a fix it's included in v3.4.

Regarding your second point:

Also I'm not sure why you say the matrix is singular only of all the elements are null.

We actually only mark the matrix as singular if there exists a row whose elements are all zero. In that case the matrix would be rank deficient and would have a determinant of zero.

Thanks again for the info, feel free to reopen this issue if you have any further problems

from basiclinearalgebra.

Hi

OK - yes indeed you check for a null row but that’s not necessary (just sufficient) to have a singular matrix. In order to return a correct status you would need to calculate the determinant - or update the documentation on this function to make it the responsibility of the caller to check this separately

from basiclinearalgebra.

Ah ok I see your point, indeed having a null row doesn't cover all singular matrices, thanks for the clarification.

I just pushed a small change (1c63223) which should correct for this. Basically in addition to the previous check we now also flag the matrix as singular if any of the diagonal elements in the LU decomposition are zero. This should be sufficient to detect singular matrices since the determinant of a matrix is the product of those diagonal elements.

As a result the Invert function now also correctly returns false for singular matrices so thanks again for the info.

from basiclinearalgebra.