I am currently using `lufact()`

in a project, but I couldn’t figure out by just reading the docs how one can check if the factorization failed. Normally, one would check if the pivots are zero, but how to do that with the LU object returned by `lufact()`

?

Also, I was reading the Linear Algebra docs and encountered the `bkfact()`

factorization. It is apparently the second option after Cholesky that the operator `\`

picks up when the matrix is not positive definite, but only symmetric. Anyone familiar with this factorization could elaborate on the advantages of using it instead of LU for example? First time reading about it.