<div dir="ltr">Hi all<div><br></div><div>I'm a bit embarrassed that after careful investigation I found that I made a wrong configuration in my problem settings. This makes the problem also not converged with MUMPS at the first place. After fixing that problem, the NR iterator converges normally with both MUMPS and Hypre. Although ill-posed, the iterative solver achieves comparable tolerance comparing with MUMPS solution.</div><div><br></div><div>Thanks again for all your helpful suggestions.</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 6:32 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"><br>
  Use -ksp_type gmres to get it to print the residuals. With preonly it doesn't compute or print them.<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
> On Jul 15, 2016, at 11:28 AM, Hoang Giang Bui <<a href="mailto:hgbk2008@gmail.com">hgbk2008@gmail.com</a>> wrote:<br>
><br>
> I used<br>
><br>
> -ksp_monitor_true_residual<br>
> -ksp_monitor_true_solution<br>
> -ksp_converged_reason<br>
><br>
> with MUMPS but it does not compute the true residual. Should I compute that myself?<br>
><br>
> Below is a sample for a full log of MUMPS<br>
> <a href="https://www.dropbox.com/s/fy5uknooxw77r19/log13Jun16_mumps?dl=0" rel="noreferrer" target="_blank">https://www.dropbox.com/s/<wbr>fy5uknooxw77r19/log13Jun16_<wbr>mumps?dl=0</a><br>
><br>
><br>
> Giang<br>
><br>
> On Fri, Jul 15, 2016 at 2:52 AM, Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</a>> wrote:<br>
><br>
><br>
> On Thu, Jul 14, 2016 at 7:29 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
> On Thu, Jul 14, 2016 at 6:27 PM, Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</a>> wrote:<br>
><br>
><br>
> Notice that there are 7 orders of magnitude between the apparent residual (using the preconditioner), and the actual residual, Ax - b.<br>
> You are using Hypre, and this generally means the Hypre coarse grid operator is crap. Please<br>
><br>
><br>
> 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?<br>
><br>
> This is why Barry and I asked for a comparsion with MUMPS. If you are right, and its just the condition number,<br>
><br>
> I said norm not condition number. I trust I'm missing something in this thread.<br>
><br>
> the LU<br>
> will not be any more accurate.<br>
><br>
>    Matt<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>
><br>
<br>
</div></div></blockquote></div><br></div></div>