<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Jul 18, 2015 at 9:53 PM, Michele Rosso <span dir="ltr"><<a href="mailto:mrosso@uci.edu" target="_blank">mrosso@uci.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>


  
  

<div>
Barry,<br>
<br>
I tried the master branch of PETSc as you suggested and unfortunately the problem with KSP being set to GMRES on the coarse level when using "-mg_coarse_ksp_type preonly"  is still there. I attached the output from ksp view and the options stack.<br>
Am I missing something?<br></div></blockquote><div><br></div><div>If you show output from -log_summary we can make sure you are running the version we expect.</div><div><br></div><div>  Thanks,</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>
Thanks,<br>
Michele<br>
<br>
<br>
On Thu, 2015-07-16 at 17:30 -0500, Barry Smith wrote:
<blockquote type="CITE">
<pre>   Michel,

    This is a very annoying feature that has been fixed in master <a href="http://www.mcs.anl.gov/petsc/developers/index.html" target="_blank">http://www.mcs.anl.gov/petsc/developers/index.html</a>  I would like to have changed it in maint but Jed would have a shit-fit :-) since it changes behavior.

  Barry

<font color="#737373">> On Jul 16, 2015, at 4:53 PM, Michele Rosso <<a href="mailto:mrosso@uci.edu" target="_blank">mrosso@uci.edu</a>> wrote:</font>
<font color="#737373">> </font>
<font color="#737373">> Hi,</font>
<font color="#737373">> </font>
<font color="#737373">> I am performing a series of solves inside a loop. The matrix for each solve changes but not enough to justify a rebuilt of the PC at each solve.</font>
<font color="#737373">> Therefore I am using  KSPSetReusePreconditioner to avoid rebuilding unless necessary. The solver is CG + MG with a custom  PC at the coarse level.</font>
<font color="#737373">> If KSP is not updated each time, everything works as it is supposed to. </font>
<font color="#737373">> When instead I allow the default PETSc  behavior, i.e. updating PC every time the matrix changes, the coarse level KSP , initially set to PREONLY, is changed into GMRES </font>
<font color="#737373">> after the first solve. I am not sure where the problem lies (my PC or PETSc), so I would like to have your opinion on this.</font>
<font color="#737373">> I attached the ksp_view for the 2 successive solve and the options stack.</font>
<font color="#737373">> </font>
<font color="#737373">> Thanks for your help,</font>
<font color="#737373">> Michel</font>
<font color="#737373">> </font>
<font color="#737373">> </font>
<font color="#737373">> </font>
<font color="#737373">> <ksp_view.txt><petsc_options.txt></font>

</pre>
</blockquote>
<br>
</div>

</blockquote></div><br><br clear="all"><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>