<p>It's easy to experiment, but I don't think it pays off.</p>
<div class="gmail_quote">On May 25, 2012 8:41 AM, "Mark F. Adams" <<a href="mailto:mark.adams@columbia.edu">mark.adams@columbia.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">I like to compartmetnalize: deal with the math and computer science separately.<div><br></div><div>And you (or someone) decided to only do one comm for SSOR.  You could, and maybe should, communicate between SOR steps.</div>
<div><div><br></div><div><div><div>On May 25, 2012, at 9:33 AM, Jed Brown wrote:</div><br><blockquote type="cite"><p>Fair point for large subdomains, but two additive have two comm steps instead of one.</p>
<div class="gmail_quote">On May 25, 2012 8:24 AM, "Mark F. Adams" <<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div style="word-wrap:break-word"><br><div><div>On May 25, 2012, at 8:59 AM, Jed Brown wrote:</div><br><blockquote type="cite"><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>
</blockquote><br></div><div>So it does a forward and backward pass.  So it is really two smoothing steps.  One should compare it with two additive smoothers.</div><br></div></blockquote></div>
</blockquote></div><br></div></div></div></blockquote></div>