<div class="gmail_quote">On Tue, May 22, 2012 at 9:52 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
How about a default nonlinear Krylov accelerator? fgrmes is not much more expensive than gmres under causal circumstances.</blockquote></div><br><div>I could handle that. The extra storage concerns me more than time.</div>
<div><br></div><div>I'd eventually like to make PCMG linear by default, especially considering that the convergence with GMRES(1) smoothers tends to be worse in my experience (except when coarse level operators are poorly scaled/unstable) and adds a ton of synchronization points.</div>
<div><br></div><div>Oh, FGMRES by default would necessarily switch from left to right preconditioning. I'm not sure that's a good thing since it requires much tighter tolerances for ill-conditioned systems to reach a given error (especially penalty boundaries). There are advantages both ways and I keep waffling. I expect that changing the default will cause some surprises.</div>