Nothing Special   »   [go: up one dir, main page]

skip to main content
article
Free access

PSBLAS: a library for parallel linear algebra computation on sparse matrices

Published: 01 December 2000 Publication History

Abstract

Many computationally intensive problems in engineering and science give rise to the solution of large, sparse, linear systems of equations. Fast and efficient methods for their soltion are very important because these systems usually occur in the innermost loop of the computational scheme. Parallelization is often necessary to achieve an acceptable level of performance. This paper presents the design, implementation, and interface of a library of Basic Linear Algebra Subroutines for sparse matrices (PSBLAS) which is specifically tailored to distributed-memory computers. PSBLAS enables easy, efficient, and portable implementations of parallel iterative solvers for linear systems. The interface keeps in view a Single Program Multiple Data programming model on distributed-memory machines. However, the architecture of the library does not exclude an implementation in different paradigms, such as those based on the shared-memory model.

References

[1]
ARIOLI, M., DUFF, I., AND RUIZ, D. 1992. Stopping criteria for iterative solvers. SIAM J. Matrix Anal. Appl. 13, 1 (Jan.), 138-144.
[2]
BALAY, S., GROPP, W., MCINNES,L.C.,AND SMITH, B. 1995. PETSc 2.0 user manual. ANL-95/11 - Revision 2.0.22. Argonne National Laboratory, Argonne, IL.
[3]
BARRETT, R., BERRY, M., CHAN, T., DEMMEL, J., DONAT, J., DONGARRA, J., EIJKHOUT, V., POZO, R., ROMINE, C., AND VORST,H.V. D. 1994. Templates for the Solution of Linear Systems. SIAM, Philadelphia, PA.
[4]
BRUASET,A.M.AND LANGTANGEN, H. P. 1997. Object-oriented design of preconditioned iterative methods in diffpack. ACM Trans. Math. Softw. 23,1,50-80.
[5]
CARNEY, S., HEROUX,M.A.,LI, G., AND WU, K. 1994. A revised proposal for a sparse BLAS toolkit. Tech. Rep. 94-034. Army High Performance Computing Research Center, Minneapolis, MN.
[6]
CERIONI, F., COLAJANNI, M., FILIPPONE, S., AND MAIOLATESI, S. 1996. PSBLAS user's guide. RI.96.11 (April). Univ. of Rome-Tor Vergata, Rome, Italy.
[7]
CHOI, J., DEMMEL, J., DHILLON, I., DONGARRA,J.J.,OSTROUCHOV, S., PETITET, A., STANLEY, K., WALKER,D.W.,AND WHALEY, R. C. 1995. ScaLAPACK: A portable linear algebra library for distributed memory computers-design issues and performance. In Proceedings of the second international workshop on Applied Parallel Computing (PARA '95, Lyngby, Den-mark), J. J. Dongarra, K. Masden, and J. Wasniewski, Eds. Springer-Verlag, New York, NY, 95-106.
[8]
DONGARRA,J.J.AND WHALEY, R. C. 1995. A user's guide to the BLACS v1.0. LAPACK Working Note 94 Technical Report CS-95-281. Department of Computer Science, University of Tennessee, Knoxville, TN. http://www.netlib.org/lapack/lawns
[9]
DONGARRA,J.J.,DU CROZ,J.J.,HAMMARLING, S., AND DUFF, I. S. 1990. A set of level 3 basic linear algebra subprograms. ACM Trans. Math. Softw. 16, 1 (Mar.), 1-17.
[10]
DONGARRA,J.J.,DU CROZ, J., HAMMARLING, S., AND HANSON, R. J. 1988. An extended set of FORTRAN basic linear algebra subprograms. ACM Trans. Math. Softw. 14, 1 (Mar.), 1-17.
[11]
DUFF,I.S.,MARRONE, M., RADICATI, G., AND VITTOLI, C. 1997. Level 3 basic linear algebra subprograms for sparse matrices: A user-level interface. ACM Trans. Math. Softw. 23,3, 379-401.
[12]
FILIPPONE, S., COLAJANNI, M., AND PASCUCCI, D. 1999. An object-oriented environment for sparse parallel computation on adaptive grids. In Proceedings of IPPS/SPDP (San Juan, Puerto Rico, Apr.).
[13]
FILIPPONE, S., MARRONE, M., AND RADICATI DI BROZOLO, G. 1992. Parallel preconditioned conjugate-gradient type algorithms for general sparsity structures. Inter. J. Comput. Math. 40, 159-167.
[14]
GAREY,M.AND JOHNSON, D. 1979. Computers and Intractability: A Guide to the Theory of NP-Completeness. W. H. Freeman and Co., New York, NY.
[15]
GREENBAUM, A. 1997. Iterative Methods for Solving Linear Systems. SIAM Frontiers in Applied Mathematics Series. SIAM, Philadelphia, PA.
[16]
HENDRICKSON,B.AND LELAND, R. 1995. An improved spectral graph partitioning algorithm for mapping parallel computations. SIAM J. Sci. Comput. 16, 2 (Mar.), 452-469.
[17]
HUTCHINSON, S., PREVOST, L., SHADID, J., TONG, C., AND TUMINARO, R. 1998. Aztec user's guide: Version 2.0. Sandia National Laboratories, Livermore, CA.
[18]
KARYPIS,G.AND KUMAR, V. 1995. METIS: Unstructured graph partitioning and sparse matrix ordering system. Computer Science Department, Univ. of Minnesota, Minneapolis, MN. http://www.cs.umn.edu/ karypis
[19]
KELLEY, C. T. 1995. Iterative Methods for Linear and Nonlinear Equations. SIAM Frontiers in Applied Mathematics Series. SIAM, Philadelphia, PA.
[20]
LAWSON,C.L.,HANSON,R.J.,KINCAID,D.R.,AND KROGH, F. T. 1979. Basic linear algebra subprograms for Fortran usage. ACM Trans. Math. Softw. 5, 3, 308-323.
[21]
MACHIELS,L.AND DEVILLE, M. O. 1997. Fortran 90: An entry to object-oriented programming for the solution of partial differential equations. ACM Trans. Math. Softw. 23, 1, 32-49.
[22]
METCALF,M.AND REID, J. 1990. Fortran 90 Explained. Oxford University Press, Inc., New York, NY.
[23]
POTHEN, A., SIMON,H.D.,AND LIOU, K.-P. 1990. Partitioning sparse matrices with eigenvectors of graphs. SIAM J. Matrix Anal. Appl. 11, 3 (July), 430-452.

Cited By

View all

Recommendations

Reviews

Maurice W. Benson

A Fortran 77 Parallel Sparse Basic Linear Algebra Subroutine (PSBLAS) library for sparse matrix operations on distributed-memory computers is presented along with a Fortran 90 object oriented user interface to this library. The focus is on operations that support iterative methods for the solution of large sparse linear systems such as those arising from the numerical solution of partial differential equations with an emphasis on preconditioned methods where the preconditioner has at most the same communication needs as matrix vector multiplication. Each processor works with a set of rows from the full system and serial sparse Basic Linear Algebra Subroutines (BLAS) are used by the processors for their serial computations. A substantial part of the development involves a clear summary of communications, data structures, PSBLAS operations and the object oriented interface to these operations. The advantages of the object oriented approach are stressed and well illustrated with a couple of examples. Efficiency is illustrated through a set of experimental results and a summary of related software is included. This work should be of interest to researchers working with parallel iterative methods as well as software developers needing parallel sparse matrix operations for an application.

Access critical reviews of Computing literature here

Become a reviewer for Computing Reviews.

Comments

Please enable JavaScript to view thecomments powered by Disqus.

Information & Contributors

Information

Published In

cover image ACM Transactions on Mathematical Software
ACM Transactions on Mathematical Software  Volume 26, Issue 4
Dec. 2000
155 pages
ISSN:0098-3500
EISSN:1557-7295
DOI:10.1145/365723
Issue’s Table of Contents

Publisher

Association for Computing Machinery

New York, NY, United States

Publication History

Published: 01 December 2000
Published in TOMS Volume 26, Issue 4

Permissions

Request permissions for this article.

Check for updates

Author Tag

  1. basic linear algebra subprograms

Qualifiers

  • Article

Contributors

Other Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

  • Downloads (Last 12 months)111
  • Downloads (Last 6 weeks)17
Reflects downloads up to 14 Jan 2025

Other Metrics

Citations

Cited By

View all
  • (2024)Alya toward exascale: algorithmic scalability using PSCToolkitThe Journal of Supercomputing10.1007/s11227-024-05989-y80:10(13533-13556)Online publication date: 1-Jul-2024
  • (2024)Software for Numerical Linear AlgebraMatrix Algebra10.1007/978-3-031-42144-0_12(607-650)Online publication date: 2024
  • (2023)FEOTS v0.0.0: a new offline code for the fast equilibration of tracers in the oceanGeoscientific Model Development10.5194/gmd-16-2795-202316:10(2795-2809)Online publication date: 24-May-2023
  • (2023)Parallel Sparse Computation ToolkitSoftware Impacts10.1016/j.simpa.2022.100463(100463)Online publication date: Jan-2023
  • (2023)Why diffusion‐based preconditioning of Richards equation works: Spectral analysis and computational experiments at very large scaleNumerical Linear Algebra with Applications10.1002/nla.252331:1Online publication date: 19-Jul-2023
  • (2022)The Linear Algebra Mapping Problem. Current State of Linear Algebra Languages and LibrariesACM Transactions on Mathematical Software10.1145/354993548:3(1-30)Online publication date: 10-Sep-2022
  • (2021)AMG Preconditioners for Linear Solvers towards Extreme ScaleSIAM Journal on Scientific Computing10.1137/20M134914X(S679-S703)Online publication date: 26-Aug-2021
  • (2020)A high-order non field-aligned approach for the discretization of strongly anisotropic diffusion operators in magnetic fusionComputer Physics Communications10.1016/j.cpc.2020.107375254(107375)Online publication date: Sep-2020
  • (2020)A parallel multithreaded sparse triangular linear system solverComputers & Mathematics with Applications10.1016/j.camwa.2019.09.01280:2(371-385)Online publication date: Jul-2020
  • (2019)A Fourier domain acceleration framework for convolutional neural networksNeurocomputing10.1016/j.neucom.2019.06.080Online publication date: Jul-2019
  • Show More Cited By

View Options

View options

PDF

View or Download as a PDF file.

PDF

eReader

View online with eReader.

eReader

Login options

Full Access

Media

Figures

Other

Tables

Share

Share

Share this Publication link

Share on social media