<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">On Wed, Oct 31, 2018 at 3:43 PM Manav Bhatia <<a href="mailto:bhatiamanav@gmail.com">bhatiamanav@gmail.com</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"><div style="word-wrap:break-word"><div>Here are the updates. I did not find the options to make much difference in the results. </div><div><br></div><div>I noticed this message in the GAMG output for cases 2, 3:  <span style="color:rgb(242,242,242);font-family:Monaco;font-size:10px;background-color:rgba(0,0,0,0.85)">HARD stop of coarsening on level 3.</span><span style="color:rgb(242,242,242);font-family:Monaco;font-size:10px;background-color:rgba(0,0,0,0.85)">  </span><span style="color:rgb(242,242,242);font-family:Monaco;font-size:10px;background-color:rgba(0,0,0,0.85)">Grid too small: 1 block nodes</span></div></div></blockquote><div><br></div><div>Yea, this is what it looks like but I can see that there are 5 levels below.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word"><div><br></div><div>Is this implying that the mesh on level 3 could not be coarsened towards levels 4/5? </div><div><br></div><div>In one of your earlier emails you had mentioned about reducing the rate of coarsening. What is the option to do that. </div></div></blockquote><div><br></div><div>The coarsening rate is pretty slow but you could try reducing it more with a finite -pc_gamg_threshold val (you could try val=0.01, 0.02, 0.05).</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word"><div><br></div><div>Common options: </div><div><div><div style="margin:0px;font-stretch:normal;font-size:10px;line-height:normal;font-family:Monaco;color:rgb(242,242,242);background-color:rgba(0,0,0,0.85)"><div style="margin:0px;font-stretch:normal;line-height:normal"><span style="font-variant-ligatures:no-common-ligatures">-pc_mg_levels 5 -mg_levels_ksp_max_it 4 -pc_gamg_square_graph 0 -pc_gamg_threshold 0. -mg_levels_ksp_type richardson -gamg_est_ksp_type cg</span></div></div></div><div><br></div><div>Case#  |  #levels   |    #KSP Iters  |  Extra Options </div><div>————————————————————————————</div><div>     1      |     5        |       67            |       </div><div>     2      |     5        |       67            |    -pc_gamg_agg_nsmooths 2 </div><div><div>     3      |     5        |       67            |    -pc_gamg_agg_nsmooths 2  -mg_levels_esteig_ksp_type cg </div></div></div></div></blockquote><div><br></div><div>I have done plates a long time ago and they worked fine, but there is something going wrong here. Perhaps my plate formulations were different but I don't think so (what we called Reissner-Mindlin).</div><div><br></div><div>You might try an easier test case if this is not a super simple test (eg, well shaped elements).</div><div><br></div><div>Otherwise you can limit the number of levels that you use and use a parallel direct solver for the (large) coarse grid. To do that use: "-pc_gamg_use_parallel_coarse_grid_solver -mg_coarse_pc_type lu" and configure your PETSc with a parallel coarse grid solver like MUMPS or SuperLU.</div><div><br></div><div>Then set -pc_gamg_coarse_eq_limit N with N being large. Increase it and minimize the solve time. You will find that parallel lu is faster than MG for "small" problems (note, small could be pretty large here). So you will want to run LU first (ie, no AMG) to get a baseline and then try GAMG.</div><div><br></div><div>Oh, and you could also try hypre (=pc_type hypre -pc_hypre_type boomeramg). Hypre should not be great in theory but it is worth a try. Maybe there's just something going wrong with GAMG.</div><div><br></div><div> </div></div></div></div></div></div></div></div>