Polynomial Based Iteration Methods for Symmetric Linear Systems (Paperback)

Bernd Fischer

ISBN 10: 3663111091 / ISBN 13: 9783663111092

Editore: Springer Fachmedien Wiesbaden, Weisbaden, 2013

Riguardo questo articolo

Paperback. 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. In many cases, the need to store the matrix factors rules out the application of direct solvers, such as Gaussian elimination in which case the only alternative is to use iterative methods. A natural way to exploit the sparsity structure of A is to design iterative schemes that involve the coefficient matrix only in the form of matrix-vector products. To achieve this goal, most iterative methods generate iterates Xn by the simple rule Xn = Xo + Qn-l(A)ro, where ro = f-Axo denotes the initial residual and Qn-l is some polynomial of degree n - 1. The idea behind such polynomial based iteration methods is to choose Qn-l such that the scheme converges as fast as possible. 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. In many cases, the need to store the matrix factors rules out the application of direct solvers, such as Gaussian elimination in which case the only alternative is to use iterative methods. A natural way to exploit the sparsity structure of A is to design iterative schemes that involve the coefficient matrix only in the form of matrix-vector products. To achieve this goal, most iterative methods generate iterates Xn by the simple rule Xn = Xo + Qn-l(A)ro, where ro = f-Axo denotes the initial residual and Qn-l is some polynomial of degree n - 1. The idea behind such polynomial based iteration methods is to choose Qn-l such that the scheme converges as fast as possible. Shipping may be from our Sydney, NSW warehouse or from our UK or US warehouse, depending on stock availability.

Codice articolo 9783663111092

Dati bibliografici

• Titolo: Polynomial Based Iteration Methods for Symmetric Linear Systems (Paperback)
• Editore: Springer Fachmedien Wiesbaden, Weisbaden
• Data di pubblicazione: 2013
• Lingua: Tedesco
• ISBN 10: 3663111091
• ISBN 13: 9783663111092
• Rilegatura: Paperback
• Condizione: new

Descrizione articolo

Riassunto

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. In many cases, the need to store the matrix factors rules out the application of direct solvers, such as Gaussian elimination in which case the only alternative is to use iterative methods. A natural way to exploit the sparsity structure of A is to design iterative schemes that involve the coefficient matrix only in the form of matrix-vector products. To achieve this goal, most iterative methods generate iterates Xn by the simple rule Xn = Xo + Qn-l(A)ro, where ro = f-Axo denotes the initial residual and Qn-l is some polynomial of degree n - 1. The idea behind such polynomial based iteration methods is to choose Qn-l such that the scheme converges as fast as possible.

Dalla quarta di copertina

Contents: Introduction - Orthogonal Polynomials - Chebyshev and Optimal Polynomials - Orthogonal Polynomials and Krylow Subspaces - Estimating the Spectrum and the Distribution function - Parameter Free Methods - Parameter Dependent Methods - The Stokes Problem - Approximating the A-Norm - Bibliography - Notation - Index