A New Block Parallel SOR Method and Its Analysis

As the development of the PSOR method (a point parallel SOR method by mesh domain partitioning proposed in [SIAM J. Sci. Comput., 20 (1999), pp. 2261-2281], this paper introduces a new mesh domain partition and ordering (the multitype partition and ordering), and proposes a new block parallel SOR (BPSOR) method for numerically solving 2-dimensional (2D) or three-dimensional (3D) elliptic boundary problems. A general mathematical analysis shows that the BPSOR method can have the same asymptotic convergence rate as the corresponding sequential block SOR method if the coefficient matrix of the block linear system is "consistently ordered." It also shows that the original sequential ordering can be maintained in the parallel implementation of the BPSOR method so that the BPSOR method can be effectively applied to solve complex elliptic boundary problems. Furthermore, three particular multitype orderings are proposed based on strip and block mesh partitions, which lead to three effective BPSOR methods for solving the five-point like linear systems (in two dimensions) and the seven-point like linear systems (in three dimensions). In addition, it is shown that the PSOR method can be generated from BPSOR if each block equation is solved approximately by only one point SOR iteration. Thus, the PSOR method is well defined without involving any interprocessor data communication operations, and extended to solving three-dimensional (3D) problems. Finally, numerical results are presented which confirm the theoretical results and show that BPSOR has a good parallel performance on a parallel MIMD computer.

[1]  I. Duff,et al.  The effect of ordering on preconditioned conjugate gradients , 1989 .

[2]  Junping Wang,et al.  Convergence analysis without regularity assumptions for multigrid algorithms based on SOR smoothing , 1992 .

[3]  B. Roux,et al.  Implicit solvent models. , 1999, Biophysical chemistry.

[4]  Loyce Adams,et al.  New Parallel SOR Method by Domain Partitioning , 1999, SIAM J. Sci. Comput..

[5]  Message P Forum,et al.  MPI: A Message-Passing Interface Standard , 1994 .

[6]  S. A. STOTLAND,et al.  Orderings for Parallel Conjugate Gradient Preconditioners , 1997, SIAM J. Sci. Comput..

[7]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[8]  C. Cramer,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics. , 1999, Chemical reviews.

[9]  Irad Yavneh,et al.  On Red-Black SOR Smoothing in Multigrid , 1996, SIAM J. Sci. Comput..

[10]  Charbel Farhat.  Multiprocessors in computational mechanics , 1987 .

[11]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[12]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[13]  Th. Lippert,et al.  A parallel SSOR preconditioner for lattice QCD , 1996, hep-lat/9608066.

[14]  James M. Ortega,et al.  Multicoloring with lots of colors , 1989, ICS '89.

[15]  David M. Young,et al.  Applied Iterative Methods , 2004 .

[16]  Henk A. van der Vorst,et al.  Large tridiagonal and block tridiagonal linear systems on vector and parallel computers , 1987, Parallel Comput..

[17]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[18]  E. L. Poole,et al.  Multicolor ICCG methods for vector computers , 1987 .

[19]  K. Stüben,et al.  Multigrid methods: Fundamental algorithms, model problem analysis and applications , 1982 .

[20]  Donald G. Truhlar,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics , 1999 .

[21]  Eric F. van de Velde,et al.  Concurrent scientific computing , 1994, Texts in applied mathematics.

[22]  James M. Ortega.  Orderings for Conjugate Gradient Preconditionings , 1991, SIAM J. Optim..

[23]  H. F. Jordan,et al.  Is SOR Color-Blind? , 1986 .

[24]  William Gropp,et al.  The MPI Message-Passing Interface Standard: Overview and Status , 1995 .

[25]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[26]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[27]  R. Grimes,et al.  On vectorizing incomplete factorization and SSOR preconditioners , 1988 .

[28]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[29]  Thomas Lippert,et al.  Parallel SSOR preconditioning for lattice QCD , 1999, Parallel Comput..

[30]  Message Passing Interface Forum.  MPI: A message - passing interface standard , 1994 .

[31]  J. Ortega,et al.  A multi-color SOR method for parallel computation , 1982, ICPP.