<div class="gmail_extra">On Wed, Dec 5, 2012 at 10:36 AM, Anton Popov <span dir="ltr"><<a href="mailto:popov@uni-mainz.de" target="_blank">popov@uni-mainz.de</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">


This is not exactly clear to me. Residuals for Dirichlet DOF should
    be (set to) zero, should they? Because those are not unknowns. What
    do you mean by Dirichlet boundary residuals?</blockquote><div><br></div><div>If you leave the variables in the vector (convenient for visualization, for example) then you have to write residuals for them. With most preconditioners, those residuals will be zero after the first iteration, but that is not required.</div>


<div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"> Is there any SNES
    example demonstrating usage of Newton or Picard for a problem with
    non-zero Dirichlet BC? </blockquote></div><br></div><div class="gmail_extra">The lid-driven cavity examples do the non-symmetric version. snes ex48 uses symmetric version enforcement and could trivially be inhomogeneous. When working with external applications, I almost always advocate doing this instead of removal or other ad-hoc boundary condition implementations.</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Here is a slightly more formal explanation:</div><div class="gmail_extra"><br></div><div class="gmail_extra"><a href="http://scicomp.stackexchange.com/a/3300/119">http://scicomp.stackexchange.com/a/3300/119</a><br>
</div>