<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Nov 9, 2017 at 12:14 AM, zakaryah . <span dir="ltr"><<a href="mailto:zakaryah@gmail.com" target="_blank">zakaryah@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 dir="ltr"><div style="font-size:small">Well the saga of my problem continues.  As I described previously in an epic thread, I'm using the SNES to solve problems involving an elastic material on a rectangular grid, subjected to external forces.  In any case, I'm occasionally getting poor convergence using Newton's method with line search.  In troubleshooting by visualizing the residual, I saw that in data sets which had good convergence, the residual was nevertheless significantly larger along the boundary between different processors.  Likewise, in data sets with poor convergence, the residual became very large on the boundary between different processors.  The residual is not significantly larger on the physical boundary, i.e. the global boundary.  When I run on a single process, convergence seems to be good on all data sets.</div><div style="font-size:small"><br></div><div style="font-size:small">Any clues to fix this?</div></div>
</blockquote></div><br>It sounds like something is wrong with communication across domains:</div><div class="gmail_extra"><br></div><div class="gmail_extra"> - If this is FEM, it sounds like you are not adding contributions from the other domain to shared vertices/edges/faces</div><div class="gmail_extra"><br></div><div class="gmail_extra"> - If this is FDM/FVM, maybe the ghosts are not updated</div><div class="gmail_extra"><br></div><div class="gmail_extra">What DM are you using? Are you using the Local assembly functions (FormFunctionLocal), or just FormFunction()?</div><div class="gmail_extra"><br></div><div class="gmail_extra">  Thanks,</div><div class="gmail_extra"><br></div><div class="gmail_extra">     Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>