Parallel iterative solvers for boundary value methods

  • Authors:
  • P. Amodio;F. Mazzia

  • Affiliations:
  • Dipartimento di Matematica, Universitá di Bari Via E. Orabona 4,1-70125 Bari, Italy;Dipartimento di Matematica, Universitá di Bari Via E. Orabona 4,1-70125 Bari, Italy

  • Venue:
  • Mathematical and Computer Modelling: An International Journal
  • Year:
  • 1996

Quantified Score

Hi-index 0.98

Visualization

Abstract

A parallel variant of the block Gauss-Seidel iteration for the solution of block-banded linear systems is presented. The coefficient matrix is partitioned among the processors as in the domain decomposition methods and then it is split so that the resulting iterative method has the same spectral properties of the block Gauss-Seidel iteration. The parallel algorithm is applied to the solution of block-banded linear systems arising from the numerical discretization of initial value problems by means of Boundary Value Methods (BVMs). BVMs define a new approach for the solution of ordinary differential equations and seem to be attractive for their interesting stability properties and a possible parallel implementation. In this paper, we refer to BVMs based on the extended trapezoidal rules.