Calculate the structural rank of the sparse matrix S.
Note that only the structure of the matrix is used in this calculation based
on a Dulmage-Mendelsohn permutation to block triangular form. As
such the numerical rank of the matrix S is bounded by
sprank (S) >= rank (S). Ignoring floating point errors
sprank (S) == rank (S).
See also: dmperm.
Package: octave