Premium
Evaluating the LCD algorithm for solving linear systems of equations arising from implicit SUPG formulation of compressible flows
Author(s) -
Catabriga Lucia,
Coutinho Alvaro L. G. A.,
Franca Leopoldo P.
Publication year - 2004
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1012
Subject(s) - generalized minimal residual method , inviscid flow , compressibility , conjugate gradient method , finite element method , system of linear equations , mathematics , linear system , algorithm , computer science , element (criminal law) , mathematical optimization , geometry , mathematical analysis , mechanics , engineering , structural engineering , physics , political science , law
Abstract In this work we evaluate the performance of the left conjugate direction method recently introduced by Yuan, Golub, Plemmons and Cecílio for the solution of non‐symmetric systems of linear equations arising from the implicit semi‐discrete SUPG finite element formulation of advective–diffusive and inviscid compressible flows. We extend the original algorithm to accommodate restarts and typical element‐by‐element preconditioners. We also show how to select the first left conjugate vector to start LCD. Several problems are solved, accessing performance parameters such as number of iterations, memory requirements and CPU times, and results are compared with other algorithms, such as GMRES, TFQMR and Bi‐CGSTAB. Copyright © 2004 John Wiley & Sons, Ltd.