STFC Website

part of UK Research & Innovation

Version 1.0.0

12th July 2004

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.