Premium
An evaluation of limited‐memory sparse linear solvers for thermo‐mechanical applications
Author(s) -
Teranishi Keita,
Raghavan Padma,
Sun Jun,
Michaleris Panagiotis
Publication year - 2007
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.2239
Subject(s) - preconditioner , cholesky decomposition , computer science , incomplete lu factorization , sparse matrix , mathematical optimization , incomplete cholesky factorization , linear system , iterative method , coefficient matrix , factorization , dimension (graph theory) , algorithm , parallel computing , matrix decomposition , mathematics , eigenvalues and eigenvectors , mathematical analysis , physics , quantum mechanics , pure mathematics , gaussian
We consider the performance of sparse linear solvers for problems that arise from thermo‐mechanical applications. Such problems have been solved using sparse direct schemes that enable robust solution at the expense of memory requirements that grow non‐linearly with the dimension of the coefficient matrix. In this paper, we consider a class of preconditioned iterative solvers as a limited‐memory alternative to direct solution schemes. However, such preconditioned iterative solvers typically exhibit complex trade‐offs between reliability and performance. We therefore characterize such trade‐offs for systems from thermo‐mechanical problems by considering several preconditioning schemes including multilevel methods and those based on sparse approximate inversion and incomplete matrix factorization. We provide an analysis of computational costs and memory requirements for model thermo‐mechanical problems, indicating that certain incomplete factorization schemes can achieve good performance. We also provide empirical evaluations that corroborate our analysis and indicate the relative effectiveness of different solution schemes. Our results indicate that our drop‐threshold incomplete Cholesky preconditioning is more robust, efficient and flexible than other popular preconditioning schemes. In addition, we propose preconditioner reuse to amortize preconditioner construction cost over a sequence of linear systems that arise from non‐linear solutions in a plastic regime. Copyright © 2007 John Wiley & Sons, Ltd.