<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, May 15, 2015 at 1:45 PM, Fande Kong <span dir="ltr"><<a href="mailto:fande.kong@colorado.edu" target="_blank">fande.kong@colorado.edu</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">Hi all,<div><br></div><div>I was wondering why we will get a zero residual when solving the linear step of a nonlinear problem? I am solving a fluid-structure interaction problems using inexact Newton method, within which fGMRES is used together with a multilevel Schwarz preconditioner. The algorithm works well when the mesh is small, but it catched a zero residual issue when we turn to use a little bit large mesh.  The output is listed as follows:</div><div><br></div><div><div>0 SNES Function norm 1.291775443677e-02 </div><div>  Jacobian: rows 5575441  Column 5575441</div><div>call DMGetInterpolation_SpmcsDMMesh </div><div>  Jacobian: rows 1298448  Column 1298448</div><div>    0 KSP Residual norm 1.291775443677e-02 </div><div>    1 KSP Residual norm 0.000000000000e+00 </div><div>Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 0</div></div></div></blockquote><div><br></div><div>This looks like a mistake in the Jacobian or preconditioner, since you have a 0 residual, but</div><div>a bad descent direction.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Any suggestions? Thanks.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Fande,</div></font></span></div>
</blockquote></div><div><br></div>-- <br><div class="gmail_signature">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></div>