<div dir="ltr">hi,<div><br></div><div>I am using KSPBCGSL method and HYPRE boomeramg precondtioner passed into KSPsolve for poisson equation as part of incompressible NS.</div><div><br></div><div>I have two questions</div><div>Q1 is number of iterations for the same solver. max_iter = 20 for the same poisson function </div><div><br></div><div>KSPGetIterationNumber() is used to set up max_iter = 20</div><div><br></div><div><br></div><div><div> 0 KSP Residual norm 4.512245447770e-04 </div><div> 2 KSP Residual norm 6.396709069731e-07 </div><div> 4 KSP Residual norm 1.757784220489e-10 </div><div> 6 KSP Residual norm 1.501464709144e-14 </div><div> 8 KSP Residual norm 1.376381429122e-18 </div><div> 10 KSP Residual norm 1.072650548945e-19 </div><div>In poisson_func(): num_iter = 10, rel_residual = 1.072651e-19 </div></div><div><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. but the solver seems to keep iterating util hitting the max_iter.</div><div><div> 0 KSP Residual norm 5.107005838093e-16 </div><div> 2 KSP Residual norm 3.353634821198e-20 </div><div> 4 KSP Residual norm 1.835096039266e-23 </div><div> 6 KSP Residual norm 1.645496102409e-23 </div><div> 8 KSP Residual norm 1.645496099837e-23 </div><div> 10 KSP Residual norm 1.645496099836e-23 </div><div> 12 KSP Residual norm 1.645496099530e-23 </div><div> 14 KSP Residual norm 1.645496034473e-23 </div><div> 16 KSP Residual norm 1.645461130961e-23 </div><div> 18 KSP Residual norm 1.645460451075e-23 </div><div>In poisson_func(): num_iter = 18, rel_residual = 1.645460e-23 </div></div><div><br></div><div>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.</div><div>is there a way to understand what's happening?<br clear="all"><div><br></div>-- <br><div class="gmail_signature">Hao Zhang<br>Dept. of Applid Mathematics and Statistics,<br>Stony Brook University, <br>Stony Brook, New York, 11790</div>
</div></div>