An introduction to ML(n)BiCGStab
Abstract
ML(n)BiCGStab is a Krylov subspace method for the solution of large, sparse and non-symmetric linear systems. In theory, it is a method that lies between the well-known BiCGStab and GMRES/FOM. In fact, when n = 1, ML(1)BiCGStab is BiCGStab and when n = N, ML(N)BiCGStab is GMRES/FOM where N is the size of the linear system. Therefore, ML(n)BiCGStab is a bridge that connects the Lanczos-based BiCGStab and the Arnoldi-based GMRES/FOM. In computation, ML(n)BiCGStab can be much more stable and converge much faster than BiCGStab when a problem with ill-condition is solved. We have tested ML(n)BiCGStab on the standard oil reservoir simulation test data called SPE9 and found that ML(n)BiCGStab reduced the total computational time by more than 60% when compared to BiCGStab. Tests made on the data from Matrix Market also support the superiority of ML(n)BiCGStab over BiCGStab. Because of the O(N^2) storage requirement in the full GMRES, one has to adopt a restart strategy to get the storage under control when GMRES is implemented. In comparison, ML(n)BiCGStab is a method with only O(nN) storage requirement and therefore it does not need a restart strategy. In this paper, we introduce ML(n)BiCGStab (in particular, a new algorithm involving A-transpose), its relations to some existing methods and its implementations.
Keywords
Cite
@article{arxiv.1106.3678,
title = {An introduction to ML(n)BiCGStab},
author = {Man-Chung Yeung},
journal= {arXiv preprint arXiv:1106.3678},
year = {2011}
}
Comments
This paper introduces the three algorithms of the ML(n)BiCGStab method for the solution of large, nonsymmetric liear systems. One of the algorithms involves A-transpose in its implmentation. This paper will be presented in the 33rd International Conference on Boundary Elements and other Mesh Reduction Methods, 28-30 June 2011, New Forest, UK