Parallel implementation of efficient preconditioned linear solver for grid-based applications in chemical physics: I: Block Jacobi diagonalization

  • Authors:
  • Wenwu Chen;Bill Poirier

  • Affiliations:
  • Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, Lubbock, TX;Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, Lubbock, TX

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2006

Quantified Score

Hi-index 31.45

Visualization

Abstract

Linear systems in chemical physics often involve matrices with a certain sparse block structure. These can often be solved very effectively using iterative methods (sequence of matrix-vector products) in conjunction with a block Jacobi preconditioner [Numer. Linear Algebra Appl. 7 (2000) 715]. In a two-part series, we present an efficient parallel implementation, incorporating several additional refinements. The present study (paper I) emphasizes construction of the block Jacobi preconditioner matrices. This is achieved in a preprocessing step, performed prior to the subsequent iterative linear solve step, considered in a companion paper (paper II). Results indicate that the block Jacobi routines scale remarkably well on parallel computing platforms, and should remain effective over tens of thousands of nodes.