Premium
Projection methods for the numerical solution of non‐self‐adjoint elliptic partial differential equations
Author(s) -
Kamath C.,
Weeratunga S.
Publication year - 1992
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.1690080104
Subject(s) - conjugate gradient method , mathematics , projection (relational algebra) , relaxation (psychology) , scalar (mathematics) , successive over relaxation , projection method , partial differential equation , iterative method , elliptic partial differential equation , block (permutation group theory) , convergence (economics) , mathematical analysis , algorithm , dykstra's projection algorithm , local convergence , geometry , psychology , social psychology , economics , economic growth
We compare the relative performances of two iterative schemes based on projection techniques for the solution of large sparse nonsymmetric systems of linear equations, encountered in the numerical solution of partial differential equations. The Block–Symmetric Successive Over‐Relaxation (Block‐SSOR) method and the Symmetric–Kaczmarz method are derived from the simplest of projection methods, that is, the Kaczmarz method. These methods are then accelerated using the conjugate gradient method, in order to improve their convergence. We study their behavior on various test problems and comment on the conditions under which one method would be better than the other. We show that while the conjugate‐gradient‐accelerated Block‐SSOR method is more amenable to implementation on vector and parallel computers, the conjugate‐gradient accelerated Symmetric–Kaczmarz method provides a viable alternative for use on a scalar machine.