MC33: Sparse unsymmetric matrix: permute to bordered block triangular form

This subroutine finds row and column permutations that reorder an \(m\) by \(n\) sparse matrix to a bordered block triangular form with full diagonal blocks, or a nested bordered block triangular form with each block itself in bordered block triangular form.