<div dir="ltr"><div dir="ltr">On Mon, Mar 16, 2020 at 4:32 PM Xiaodong Liu <<a href="mailto:xliu29@ncsu.edu">xliu29@ncsu.edu</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi, all, <br></div><div><br></div><div>I am testing KSPSetInitialGuessNonzero using the case <br></div><div><br></div><div><a href="https://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex50.c.html" target="_blank">https://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex50.c.html</a></div><div><br></div><div>This case is special. WIth zero initial guess, 1 iteration can deliver the exact solution. But when I tried true for KSPSetInitialGuessNonzero, it shows the same convergence history as false. <br></div><div><br></div><div>Does the code change KSPSetInitialGuessNonzero to false automatically?</div></div></blockquote><div><br></div><div>For that flag to make any difference, you have to pass in a nonzero vector to KSPSolve(). This example does not do that.</div><div>Were you doing that?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Take care.<br></div><div><br></div><div>Thanks,</div><div><br></div><div><div><div dir="ltr"><div dir="ltr"><div>Xiaodong Liu, PhD<br>X: Computational Physics Division<br>Los Alamos National Laboratory<br>P.O. Box 1663, <br>Los Alamos, NM 87544<br>505-709-0534<br></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>