Avoiding breakdown in the CGS algorithm

The conjugate gradient squared algorithm can suffer of similar breakdowns as Lanczos type methods for the same reason that is the non-existence of some formal orthogonal polynomials. Thus curing such breakdowns is possible by jumping over these non-existing polynomials and using only those of them which exist. The technique used is similar to that employed for avoiding breakdowns in Lanczos type methods. The implementation of these new methods is discussed. Numerical examples are given.