symamd: Symmetric approximate minimum degree permutation. p = symamd(S) for a symmetric positive definite matrix S, returns the permutation vector p such that S(p,p) tends to have a sparser Cholesky factor than S. To find the ordering for S, symamd constructs a matrix M such that spones(M’*M) = spones (S), and then computes p = colamd(M). The symamd function may also work well for symmetric indefinite matrices.