x = A b;
- Is A square?
no => use QR to solve least squares problem.
- Is A triangular or permuted triangular?
yes => sparse triangular solve
- Is A symmetric with positive diagonal elements?
yes => attempt Cholesky after symmetric minimum degree.
- Otherwise
=> use LU on A (:, colamd(A))
