In this paper, the block generalized product-type bi-conjugate gradient (GPBi-CG) method for solving large, sparse nonsymmetric linear systems of equations with multiple right-hand sides is proposed. The new algorithm is based on the block BiCG process.We analyze the convergence behavior of this method and present a bound for the residual norm of block GPBi-CG according to the residual norm of Bl-GMRES method. In addition, we prove that convergence is guaranteed when A is positive real. The numerical experiments show the efficiency of the new method and confirm the theoretical results.

Keywords