<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Mar 18, 2022 at 11:44 AM Fabio Durastante <<a href="mailto:fabio.durastante@unipi.it">fabio.durastante@unipi.it</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">Hi everybody,<br>
<br>
I'm trying to run the rotated anisotropy example ex54f using CG and GAMG as preconditioner, I run it with the command:<br>
<br>
mpirun -np 2 ./ex54f -ne 1011 \<br>
-theta 18.0 \<br>
-epsilon 100.0 \<br>
-pc_type gamg \<br>
-pc_gamg_type agg \<br>
-log_view \<br>
-log_trace \<br>
-ksp_view \<br>
-ksp_monitor \<br>
-ksp_type cg \<br>
-mg_levels_pc_type jacobi \<br>
-mg_levels_ksp_type richardson \<br>
-mg_levels_ksp_max_it 4 \<br>
-ksp_atol 1e-9 \<br>
-ksp_rtol 1e-12<br>
<br>
But the KSP CG seems to stop just after two iterations:<br>
<br>
   0 KSP Residual norm 6.666655711717e-02<br>
   1 KSP Residual norm 9.859661350927e-03<br>
<br>
I'm attaching the full log, the problem seems to appear when I modify the value of epsilon, if I leave it to the default (1.0) it prints<br>
<br>
   0 KSP Residual norm 5.862074869050e+00<br>
   1 KSP Residual norm 5.132711016122e-01<br>
   2 KSP Residual norm 1.198566629717e-01<br>
   3 KSP Residual norm 1.992885901625e-02<br>
   4 KSP Residual norm 4.919780086064e-03<br>
   5 KSP Residual norm 1.417045143681e-03<br>
   6 KSP Residual norm 3.559622318760e-04<br>
   7 KSP Residual norm 9.270786187701e-05<br>
   8 KSP Residual norm 1.886403709163e-05<br>
   9 KSP Residual norm 2.940634415714e-06<br>
  10 KSP Residual norm 5.015043022637e-07<br>
  11 KSP Residual norm 9.760219712757e-08<br>
  12 KSP Residual norm 2.320857464659e-08<br>
  13 KSP Residual norm 4.563772507631e-09<br>
  14 KSP Residual norm 8.896675476997e-10<br>
<br>
that is very strange because the case with epsilon 1 should be easier.<br></blockquote><div><br></div><div>As Barry said, it failed (use -ksp_converged_reason). </div><div><br></div><div>You are making it pretty hard with 100. The example uses .1 (ie, 10)</div><div> <br></div><div>You are adding -mg_levels_ksp_type richardson which does work with jacobi. Use sor (instead of jacobi) and it should work.</div><div>The default ksp_type is chebyshev.</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
Any help with this would be great.<br>
<br>
Thank you very much,<br>
<br>
Fabio Durastante<br>
</blockquote></div></div>