<div dir="ltr">thanks.<div><br></div><div><span style="font-size:12.8px">Yes. KSPSetTolerances is used.</span><br></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">tolerance is set to be 1.e-14</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">I was thinking if there is a debug option I could use from PETSc to get some insight?  GDB was used.</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 17, 2017 at 12:25 AM, 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 Sep 16, 2017, at 10:39 PM, Hao Zhang <<a href="mailto:hbcbh1999@gmail.com">hbcbh1999@gmail.com</a>> wrote:<br>
><br>
> hi,<br>
><br>
> I am using KSPBCGSL method and HYPRE boomeramg precondtioner passed into KSPsolve for poisson equation as part of incompressible NS.<br>
><br>
> I have two questions<br>
> Q1 is number of iterations for the same solver. max_iter = 20 for the same poisson function<br>
><br>
> KSPGetIterationNumber() is used to set up max_iter = 20<br>
<br>
</span>  What do you mean, KSPGetIterationNumber() tells you the current number of iterations. You cannot use it to set the maximum number of iterations. Use KSPSetTolerances() or -ksp_max_it <its> to set the number<br>
<span class="">><br>
><br>
>  0 KSP Residual norm 4.512245447770e-04<br>
>   2 KSP Residual norm 6.396709069731e-07<br>
>   4 KSP Residual norm 1.757784220489e-10<br>
>   6 KSP Residual norm 1.501464709144e-14<br>
>   8 KSP Residual norm 1.376381429122e-18<br>
>  10 KSP Residual norm 1.072650548945e-19<br>
> In poisson_func(): num_iter = 10, rel_residual = 1.072651e-19<br>
><br>
> and the very same function at the next time step gives: remind you guys max_iter = 20. I assume the number of iterations should stop at num_iter = 2.<br>
<br>
</span>   Why would it stop at num_iter 2. The number of iterations allowed is the same for each solve, it does not accumulate over all solves.<br>
<span class=""><br>
>  but the solver seems to keep iterating util hitting the max_iter.<br>
>  0 KSP Residual norm 5.107005838093e-16<br>
>   2 KSP Residual norm 3.353634821198e-20<br>
>   4 KSP Residual norm 1.835096039266e-23<br>
>   6 KSP Residual norm 1.645496102409e-23<br>
>   8 KSP Residual norm 1.645496099837e-23<br>
>  10 KSP Residual norm 1.645496099836e-23<br>
>  12 KSP Residual norm 1.645496099530e-23<br>
>  14 KSP Residual norm 1.645496034473e-23<br>
>  16 KSP Residual norm 1.645461130961e-23<br>
>  18 KSP Residual norm 1.645460451075e-23<br>
> In poisson_func(): num_iter = 18, rel_residual = 1.645460e-23<br>
<br>
</span>   But the way it is almost never realistic to try to get residual norms to be this small. You might consider using something like -ksp_atol 1.e-16 or something<br>
<span class=""><br>
><br>
> Q2 is: it seems HYPRE with KSPsolve works with coarse mesh but when I double the mesh size. the solve halted indefinitely. ksp_monitor didn't produce any results.<br>
> is there a way to understand what's happening?<br>
<br>
</span>   You could run in the debugger and use control c to interrupt the code when it is "hanging" and type bt to see the stack trace for the hanging.<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
><br>
> --<br>
> Hao Zhang<br>
> Dept. of Applid Mathematics and Statistics,<br>
> Stony Brook University,<br>
> Stony Brook, New York, 11790<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Hao Zhang<br>Dept. of Applid Mathematics and Statistics,<br>Stony Brook University, <br>Stony Brook, New York, 11790</div>
</div>