The Design of Data-Structure-Neutral Libraries for the Iterative Solution of Sparse Linear Systems

Barry F. Smith, William D. Gropp

Research output: Contribution to journalArticlepeer-review

Abstract

Over the past few years several proposals have been made for the standardization of sparse matrix storage formats in order to allow for the development of portable matrix libraries for the iterative solution of linear systems. We believe that this is the wrong approach. Rather than define one standard (or a small number of standards) for matrix storage, the community should define an interface (i.e., the calling sequences) for the functions that act on the data. In addition, we cannot ignore the interface to the vector operations because, in many applications, vectors may not be stored as consecutive elements in memory. With the acceptance of shared memory, distributed memory, and cluster memory parallel machines, the flexibility of the distribution of the elements of vectors is also extremely important. This issue is ignored in most proposed standards. In this article we demonstrate how such libraries may be written using data encapsulation techniques.

Original languageEnglish (US)
Pages (from-to)329-336
Number of pages8
JournalScientific Programming
Volume5
Issue number4
DOIs
StatePublished - 1996
Externally publishedYes

ASJC Scopus subject areas

  • Software
  • Computer Science Applications

Fingerprint

Dive into the research topics of 'The Design of Data-Structure-Neutral Libraries for the Iterative Solution of Sparse Linear Systems'. Together they form a unique fingerprint.

Cite this