<div class="gmail_quote">On Sun, Oct 9, 2011 at 22:54, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">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;">
<br>
   Jed,<br>
<br>
       It appears your  some of your many<br>
<br>
if (s_t >= x*y*z) {s_t = bases[n26] + (i-1)*x_t;} /* 2d case */<br>
<br>
things are breaking the communication with a stencil width of 2<br></blockquote><div><br></div><div>da3.c:980</div><div>M-x xhg-ann<RET></div><div><br></div><div>Looks like that line comes from</div><div><br></div>
<div><div>changeset:   18790:4310a5aad280</div><div>user:        Ethan Coon <<a href="mailto:ecoon@lanl.gov">ecoon@lanl.gov</a>></div><div>date:        Wed Apr 20 00:57:49 2011 +0200</div><div>summary:     cruft to handle the special case of solving a 2D problem in a 3D algorithm with stencil size > 1.  This only works in the 3D case, z-dimension</div>
</div><div><br></div><div><a href="http://petsc.cs.iit.edu/petsc/petsc-dev/rev/4310a5a">http://petsc.cs.iit.edu/petsc/petsc-dev/rev/4310a5a</a></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

the second layer of ghost points uses the first layer instead of getting the second layer.<br>
<br>
Since I don't understand your /* 2d case */ business I leave it to you to fix<br></blockquote><div><br></div><div>I don't understand it either and haven't tried. Maybe Ethan can take a look.</div><div><br></div>
</div>