In the paper we consider the iterative solution of linear system Ax=b by the conjugate gradient method using the factorized preconditioner B=(I+LZ)Y(I+ZU), where A=D+L+U is the additive splitting of the coefficient matrix into the strictly lower triangular, the diagonal, and the strictly upper triangular parts. For an arbitrary symmetric positive definite matrix A, the diagonal matrices Y>0 and Z are constructed as the minimizers of a certain upper bound for the K-condition number of the inverse preconditioned matrix. The main advantages of the new method are as follows: wide range of applicability, low operation number count per iteration, good parallelizability for all the stages of computation, and sufficient reduction of the iteration number (for a properly chosen preconditioning parameters). Numerical results are given for several test problems.

Сonjugate gradient method, factorized preconditioner, the K-condition number.

Mathematical problems and theory of numerical methods