MathOverflow is a question and answer site for professional mathematicians. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

In usual formulation of conjugate gradient algorithm initial search direction is taken to be the residual (so residual and search direction spans Krylov subspace). However, in cases where inexact Krylov is used, search direction is not anymore in span of Krylov subspace. I am in process of analyzing convergence of one such version of CG where search directions are not anymore in span of Krylov subspace of initial residual. I wonder if any such method has been already discussed in literature. In particular, what will be condition of any instance of search direction (more than just being a descent direction) for CG to successfully converge. How does the convergence rate depends on initial search direction?

share|cite|improve this question
up vote 1 down vote accepted

I recommend Simoncini and Szyld, Theory of inexact Krylov subspace methods and applications to scientific computing, SIAM J. Sci. Comput., 25 (2003) pp. 454-477. It looks like it discusses exactly what you want.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.