Abstract
It is well-known that Conjugate Gradient (CG) method with DIC factorization can be viewed as one of useful solvers for a linear system of equations with a symmetric positive definite matrix. On the other hand, Eisenstat’s trick also can be viewed as an efficient implementation for CG methed without preconditioning only. However, this technique had limited effectiveness because of insufficiency and instability of DIC factorization. In this paper, we will propose a strategy for a more efficient and stable implementation combined with accelerated DIC factorization and Eisenstat’s trick and by numerical experiments we will present that this indeed may lead to faster and stable convergence.
