<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Apr 10, 2015 at 1:05 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
> On Apr 10, 2015, at 12:58 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
><br>
> On Fri, Apr 10, 2015 at 12:57 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
> > On Apr 10, 2015, at 8:07 AM, Klaij, Christiaan <<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a>> wrote:<br>
> ><br>
> > Barry,<br>
> ><br>
> > Sure, I can call PCFieldSplitGetSubKSP() to get the fieldsplit_0<br>
> > ksp and then KSPGetIterationNumber, but what does this number<br>
> > mean?<br>
> ><br>
> > It appears to be the number of iterations of the last time that<br>
> > the subsystem was solved, right? If so, this corresponds to the<br>
> > last iteration of the coupled system, how about all the previous<br>
> > iterations?<br>
><br>
>    If I added a KSPGetTotalIterations() would that help?<br>
><br>
> Its hard to see how to define this consistently.<br>
<br>
</span>  Just count the number of iterations that KSP object has done since it was created. How could that be inconsistent?</blockquote><div><br></div><div>Oh, you want the counts to be localized in the subKSP. This not nearly as useful. For example, we could have left this</div><div>count in the KSP, but we hoist it to the SNES. I am fine totaling the counts in the KSP, but I would propose custom</div><div>monitors as a way of hoisting the info to the next level.</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 class="HOEnZb"><div class="h5">
> I think this is better handled by custom monitors.<br>
> We could provide a default monitor for FieldSplit which does this.<br>
><br>
>   Matt<br>
><br>
><br>
> Barry<br>
><br>
> ><br>
> > Chris<br>
> > ________________________________________<br>
> > From: Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>><br>
> > Sent: Friday, April 10, 2015 2:48 PM<br>
> > To: Klaij, Christiaan<br>
> > Cc: <a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><br>
> > Subject: Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1<br>
> ><br>
> >   Chris,<br>
> ><br>
> >     It appears you should call PCFieldSplitGetSubKSP() and then get the information you want out of the individual KSPs. If this doesn't work please let us know.<br>
> ><br>
> >   Barry<br>
> ><br>
> >> On Apr 10, 2015, at 6:48 AM, Klaij, Christiaan <<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a>> wrote:<br>
> >><br>
> >> A question when using PCFieldSplit: for each linear iteration of<br>
> >> the system, how many iterations for fielsplit 0 and 1?<br>
> >><br>
> >> One way to find out is to run with -ksp_monitor,<br>
> >> -fieldsplit_0_ksp_monitor and -fieldsplit_0_ksp_monitor. This<br>
> >> gives the complete convergence history.<br>
> >><br>
> >> Another way, suggested by Matt, is to use -ksp_monitor,<br>
> >> -fieldsplit_0_ksp_converged_reason and<br>
> >> -fieldsplit_1_ksp_converged_reason. This gives only the totals<br>
> >> for fieldsplit 0 and 1 (but without saying for which one).<br>
> >><br>
> >> Both ways require to somehow process the output, which is a bit<br>
> >> inconvenient. Could KSPGetResidualHistory perhaps return (some)<br>
> >> information on the subsystems' convergence for processing inside<br>
> >> the code?<br>
> >><br>
> >> Chris<br>
> >><br>
> >><br>
> >> dr. ir. Christiaan Klaij<br>
> >> CFD Researcher<br>
> >> Research & Development<br>
> >> E mailto:<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a><br>
> >> T <a href="tel:%2B31%20317%2049%2033%2044" value="+31317493344">+31 317 49 33 44</a><br>
> >><br>
> >><br>
> >> MARIN<br>
> >> 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands<br>
> >> T <a href="tel:%2B31%20317%2049%2039%2011" value="+31317493911">+31 317 49 39 11</a>, F <a href="tel:%2B31%20317%2049%2032%2045" value="+31317493245">+31 317 49 32 45</a>, I <a href="http://www.marin.nl" target="_blank">www.marin.nl</a><br>
> >><br>
> ><br>
><br>
><br>
><br>
><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>
<br>
</div></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>