Polynomial Based Iteration Methods for Symmetric Linear Systems

Any book on the solution of nonsingular systems of equations is bound to start with Ax= J, but here, A is assumed to be symmetric. These systems arise frequently in scientific computing, for example, from the discretization by finite differences or by finite elements of partial differential equations. Usually, the resulting coefficient matrix A is large, but sparse.