<div class="gmail_quote">On Fri, May 25, 2012 at 1:07 AM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>And, I've never seen Gauss-Siedel used with Cheby because G-S has the correct damping properties, as is, for the Laplacian.</div></blockquote><div><br></div><div>The point is to have something adequate for things that are not Laplacians. I tried running SOR without Cheby, but it was far less robust.</div>
<div><br></div><div>So I know it looks funny, but I don't have a similarly robust alternative. If we are living in a world where local work is cheap, we might as well do local SOR instead of pbjacobi. (Note that Cheby+pbjacobi is nearly as good as Cheby+SOR in some cases, but much worse in others.)</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br></div><div>Note, G-S is not symmetric and Cheby for unsymmetric is a different can of worms.  So if A is symmetric then maybe try SSOR.</div>
</blockquote></div><br><div>The default SOR is local_symmetric.</div>