<div dir="ltr"><div dir="ltr">* Do you set:<div><br></div><div>    PetscCall(MatSetOption(Amat, MAT_SPD, PETSC_TRUE));<br>    PetscCall(MatSetOption(Amat, MAT_SPD_ETERNAL, PETSC_TRUE));<br></div><div><br></div><div>Do that to get CG Eigen estimates. Outright failure is usually caused by a bad Eigen estimate.</div><div>-pc_gamg_esteig_ksp_monitor_singular_value<br></div><div>Will print out the estimates as its iterating. You can look at that to check that the max has converged.</div><div><br></div><div>*  -pc_gamg_aggressive_coarsening 0<br></div><div><br></div><div>will slow coarsening as well as threshold.</div><div><br></div><div>* you can run with '-info :pc' and send me the output (grep on GAMG)</div><div><br></div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Mar 24, 2023 at 2:47 PM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">You can -pc_gamg_threshold .02 to slow the coarsening and either stronger smoother or increase number of iterations used for estimation (or increase tolerance). I assume your system is SPD and you've set the near-null space.<br>
<br>
Blaise Bourdin <<a href="mailto:bourdin@mcmaster.ca" target="_blank">bourdin@mcmaster.ca</a>> writes:<br>
<br>
> Hi,<br>
><br>
> I am having issue with GAMG for some very ill-conditioned 2D linearized elasticity problems (sharp variation of elastic moduli with thin  regions of nearly incompressible material). I use snes_type newtonls, linesearch_type cp, and pc_type gamg without any further options. pc_type Jacobi converges fine (although slowly of course).<br>
><br>
><br>
> I am not really surprised that gamg would not converge out of the box, but don’t know where to start to investigate the convergence failure. Can anybody help?<br>
><br>
> Blaise<br>
><br>
> — <br>
> Canada Research Chair in Mathematical and Computational Aspects of Solid Mechanics (Tier 1)<br>
> Professor, Department of Mathematics & Statistics<br>
> Hamilton Hall room 409A, McMaster University<br>
> 1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada <br>
> <a href="https://www.math.mcmaster.ca/bourdin" rel="noreferrer" target="_blank">https://www.math.mcmaster.ca/bourdin</a> | +1 (905) 525 9140 ext. 27243<br>
</blockquote></div></div>