<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Mar 24, 2017 at 8:56 AM, alexandre this <span dir="ltr"><<a href="mailto:alexandre.this@gmail.com" target="_blank">alexandre.this@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">Dear all,<div><br></div><div>I'm in the configuration where I need to solve a linear system Ax = b where b is very small although not equal to zero.</div><div><br></div><div>It appears that, in this configuration, the initial residual of the first iteration of the solver is completely out of reach and the KSP_DIVERGED_DTOL is raised.</div></div></blockquote><div><br></div><div>For any convergence question, we need the output from</div><div><br></div><div>  -ksp_monitor_true_residual -ksp_converged_reason -ksp_view</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>In particular, when using the "-ksp_monitor_true_residual flag", I get the following :</div><div>0 KSP preconditioned resid norm 2.907111674781e+00 true resid norm 1.955854211540e+02 ||r(i)||/||b|| 1.765818923254e+09</div></div></blockquote><div><br></div><div>We need to see the error, not just the first line.</div><div><br></div><div>You can always raise the tolerance, but divergence means that the norm is increasing.</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 dir="ltr"><div>What is the best course of action in this case ? </div><div><br></div><div>Best,</div><div>Alexandre</div><span class="HOEnZb"><font color="#888888"><div><br></div><div><br>-- <br><div class="m_6334818203009088624gmail_signature"><div dir="ltr"><div><div dir="ltr"><div style="color:rgb(34,34,34);font-family:"times new roman","new york",times,serif;font-size:16px;background-color:rgb(255,255,255)"><div>Alexandre This</div><div><br></div><div>Doctorant CIFRE (Philips Healthcare - INRIA) : <span style="color:rgb(0,0,0);font-family:"times new roman";font-size:medium">Fusion Image / Modèles Numériques pour la </span><span style="color:rgb(0,0,0);font-family:"times new roman";font-size:medium">quantification de la sévérité de la régurgitation m</span><span style="color:rgb(0,0,0);font-family:"times new roman";font-size:medium">itrale</span></div><div><br></div><div>Master 2 Recherche Mathématiques, Vision, Apprentissage - ENS Cachan</div><div>Ingénieur spécialisation Santé&Technologie - ECE Paris<br></div><div>Bachelor of Science : Informatics - Aalborg University, Danemark</div><div><br></div></div><div style="color:rgb(34,34,34);font-family:"times new roman","new york",times,serif;font-size:16px;background-color:rgb(255,255,255)">tel : 06.32.57.12.44</div></div></div></div></div>
</div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" 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>
</div></div>