MC37: Sparse symmetric matrix: represent as a sum of element matrices

Given a sparse symmetric matrix \(\mathbf{A}\), this subroutine computes a set of element matrices that, if assembled, would yield the same matrix. Note that this set of elements is not unique. The matrix can be input by the user either in compressed column format (column pointer/row index scheme) or by row and column index pairs in any order.