<div class="gmail_quote">On Tue, Dec 21, 2010 at 19:53, Vijay S. Mahadevan <span dir="ltr"><<a href="mailto:vijay.m@gmail.com">vijay.m@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":25d">I am running a linear problem discretized with FEM on a diffusion<br>
reaction system, with discontinuous source distribution. When I run<br>
FGMRes with geometric multigrid as its preconditioner, I notice that<br>
every time after the restart in fgmres, the new residual is orders of<br>
magnitude higher than the previous iteration. I might be wrong on this<br>
but should the restart not preserve monotonicity in convergence ? Or<br>
am I thinking of a different variant of Gmres here.</div></blockquote></div><br><div>It is not possible to guarantee monotonicity for nonsymmetric matrices without storing the full subspace. There is no variant of GMRES, or any Krylov method for that matter, that can do what you want. You are seeing a particularly large jump, if you actually have a linear preconditioner (if you don't use Krylov cycles inside your smoothers) then you might try using bcgs or some variant thereof which would avoid the high cost of restart. Or you could stop using restarts, it looks like you were getting close to an adequate tolerance. Or find a way to make the preconditioner strong enough to converge in a reasonable number of iterations.</div>