Permutation Pr represented as an array of indices.
The row permutation matrix can be reconstructed via:
>>> Pr = np.zeros((n, n)) >>> Pr[perm_r, np.arange(n)] = 1
scipy.sparse.linalg.SuperLU.perm_c
scipy.sparse.linalg.SuperLU.L