<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Mar 14, 2018 at 11:21 PM, Sonia Pozzi <span dir="ltr"><<a href="mailto:sonia.pozzi@usi.ch" target="_blank">sonia.pozzi@usi.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Barry,<div><br></div><div>thank you for the answer. That helped a lot. Just a second curiosity. </div><div>I’m setting A00 to be solved with preonly+lu. I obtain the following</div><div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">ksp_0  KSPGetTotalIterations: 26</span></div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">ksp_1  KSPGetTotalIterations: 22</span></div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Residual ksp_0: 0 Reason ksp_0: 4</span></div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Solution ksp_0 :  Convergence in 1 iterations.</span></div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Residual ksp_1: 5.23074e-14 Reason ksp_1: 2</span></div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Solution ksp_1 :  Convergence in 11 iterations.</span></div><div style="margin:0px;font-size:13px;line-height:normal;font-family:'Andale Mono';color:rgb(41,249,20);background-color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Initial system linear residual:  9.904804e-09 </span></div></div><div><br></div><div>To me is clear that kps_0 iterations number is equal to 22*1+2*2. </div><div>But why is the iterations number for the schur complement multiplied by two?</div><div>I’m giving the mass matrix as a preconditioner for the schur complement.</div></div></blockquote><div><br></div><div>Possible cause:</div><div><br></div><div>  The preconditioner is applied once before the solve happens in order to compute the preconditioned initial residual. Then</div><div>  if you do just one iterate in the solve, you would have two applications of the preconditioner.</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 style="word-wrap:break-word"><div>Sonia</div><div><span class=""><br><div>
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)">Ph.D. Student</font></div><div style="margin:0px"><font color="#6f6f6f">Group of Prof. Rolf Krause</font></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)"><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)">Institute of Computational Science</font></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)">Università della Svizzera italiana</font></div></font></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)">Center for Computation Medicine in Cardiology</font></div><div style="color:rgb(0,0,0);margin:0px"><span style="color:rgb(111,111,111)"><a href="https://maps.google.com/?q=Via+Giuseppe+Buffi,+13+CH-6900+Lugano+Switzerland&entry=gmail&source=g">Via Giuseppe Buffi, 13</a></span></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)"><a href="https://maps.google.com/?q=Via+Giuseppe+Buffi,+13+CH-6900+Lugano+Switzerland&entry=gmail&source=g">CH-6900 Lugano</a></font></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)"><a href="https://maps.google.com/?q=Via+Giuseppe+Buffi,+13+CH-6900+Lugano+Switzerland&entry=gmail&source=g">Switzerland</a></font></div><div style="color:rgb(111,111,111);margin:0px;line-height:normal;min-height:12px"><br></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)">Email:<span class="m_-3896794662654634086Apple-converted-space"> </span><a href="mailto:pozzis@usi.ch" target="_blank">pozzis@usi.ch</a> </font></div><div style="color:rgb(0,0,0);margin:0px"><font face="Helvetica" color="#6f6f6f" style="line-height:normal;color:rgb(111,111,111)">Phone: </font><span style="color:rgb(111,111,111)"><a href="tel:+41%2058%20666%2049%2072" value="+41586664972" target="_blank">+41 58 666 4972</a></span></div></div></div>
</div>
<br></span><div><blockquote type="cite"><div>On 14 Mar 2018, at 14:52, Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:</div><br class="m_-3896794662654634086Apple-interchange-newline"><div><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">KSPGetIterationNumber</span></div></blockquote></div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><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.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>