An element-by-element preconditioned conjugate gradient method implemented on a vector computer

✍ Scribed by Jocelyne Erhel; Alice Traynard; Marina Vidrascu

• Publisher: Elsevier Science
• Year: 1991
• Tongue: English
• Weight: 444 KB
• Volume: 17
• Category: Article
• ISSN: 0167-8191
• DOI: 10.1016/s0167-8191(05)80050-x

✦ Synopsis

We consider the linear equation Ax = b where A is a sparse symmetric positive definite matrix arising from a finite element discretisation. We use the preconditioned conjugate gradient method to solve this equation, introducing an element-by-element preconditioner which is based on a Crout's decomposition of the element matrices and an element-by-element product of them. When the mesh is coloured, this preconditioner is largely vectorizable. We implement this method on a CRAY-2, and test it on 2D and 3D elastic and thermal problems and compare it to other classical preconditioners.