The BI-CG method and its variants such as ccs, BI-CGSTAB, and BI-CGSTA
B2 for solving nonsymmetric linear systems of equations can suffer fro
m breakdown. By using the theory of the formal orthogonal polynomials
we present a breakdown-free BI-CGSTAB algorithm and a breakdown-free B
I-CGSTAB2 algorithm. Only exact breakdowns are cured accurately. (C) 1
997 Elsevier Science Inc.