<div dir="ltr">I used<div><br></div><div><div>-ksp_monitor_true_residual</div><div>-ksp_monitor_true_solution</div><div>-ksp_converged_reason</div></div><div><br></div><div>with MUMPS but it does not compute the true residual. Should I compute that myself?</div><div><br></div><div>Below is a sample for a full log of MUMPS</div><div><a href="https://www.dropbox.com/s/fy5uknooxw77r19/log13Jun16_mumps?dl=0">https://www.dropbox.com/s/fy5uknooxw77r19/log13Jun16_mumps?dl=0</a><br></div><div><br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Giang</div></div></div>
<br><div class="gmail_quote">On Fri, Jul 15, 2016 at 2:52 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Thu, Jul 14, 2016 at 7:29 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div><div><div class="gmail_quote">On Thu, Jul 14, 2016 at 6:27 PM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br><div><br></div><div>Notice that there are 7 orders of magnitude between the apparent residual (using the preconditioner), and the actual residual, Ax - b.</div><div>You are using Hypre, and this generally means the Hypre coarse grid operator is crap. Please</div><div><br></div></div></div></div></div></blockquote><div><br></div></span><div>Huh?, this data looks fine, both the true and preconditioned residual stay separated by about 9 orders of magnitude. This just tells you that the norm of A (or is it A^-1) is 10^9.  Am I misunderstanding this?</div></div></div></div>
</blockquote></div><br></div></div>This is why Barry and I asked for a comparsion with MUMPS. If you are right, and its just the condition number, </div></div></blockquote><div><br></div></span><div>I said norm not condition number. I trust I'm missing something in this thread.</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra">the LU</div><div class="gmail_extra">will not be any more accurate.</div><div class="gmail_extra"><br></div><div class="gmail_extra">   Matt<span><br clear="all"><div><br></div>-- <br><div data-smartmail="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>
</span></div></div>
</blockquote></span></div><br></div></div>
</blockquote></div><br></div></div>