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.
- Publication:
-
arXiv e-prints
- Pub Date:
- June 2011
- DOI:
- 10.48550/arXiv.1106.3678
- arXiv:
- arXiv:1106.3678
- Bibcode:
- 2011arXiv1106.3678Y
- Keywords:
-
- Mathematics - Numerical Analysis;
- Computer Science - Numerical Analysis;
- Mathematics - Optimization and Control;
- 65D15;
- 65F30;
- 65G20;
- 68W10
- E-Print:
- 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