On Mon, Nov 29, 2010 at 5:04 PM, Verena Kuhlemann <span dir="ltr">&lt;<a href="mailto:vkuhlem@emory.edu">vkuhlem@emory.edu</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hello,<div><br></div><div>I have a (hopefully) quick question about setting the restart for gmres.</div><div><br></div><div>I set the restart to be 50, but when I use KSPView it looks like GMRES</div><div>still uses a restart of 30. It is possible that I overwrite my settings </div>

<div>somehow? Here is the relevant part of my code:</div><div><br></div><div>KSPCreate(PETSC_COMM_WORLD,&amp;ksp);</div><div>KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN);</div><div>KSPGMRESSetRestart(ksp,50);</div>
</blockquote><div><br></div><div>This is ignored because you have not yet set the KSP type.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div>
KSPSetInitialGuessNonzero(ksp,PETSC_TRUE);</div><div>KSPSetTolerances(ksp,1.e-12,1.e-12,PETSC_DEFAULT,100000);</div><div>KSPSetUp(ksp);</div></blockquote><div><br></div><div>If you move this before the restart call, the type wil be set and it will work.</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>KSPGetPC(ksp,&amp;pc);</div><div>PCSetType(pc,PCNONE);</div><div>KSPView(ksp,PETSC_VIEWER_STDOUT_WORLD);</div>

<div>KSPSolve(ksp,b,x);</div><div><br></div><div><br></div><div>Thanks,</div><div>Verena</div>
</blockquote></div><br><br clear="all"><br>-- <br>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<br>