<div class="gmail_quote">On Wed, Feb 16, 2011 at 15:17, Juha Jäykkä <span dir="ltr">&lt;<a href="mailto:juhaj@iki.fi">juhaj@iki.fi</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":4im">even<br>
though the eventual solution returned by SNESSolve() is correct.</div></blockquote></div><br><div>SNESSolve uses a Newton method so the linear system is being solving for a defect. If the initial guess is zero, then it would normally pick up your Dirichlet boundary conditions on the first iteration and all subsequent solves would have zero in those locations.</div>