<div dir="ltr">Thanks, Patrick. I will try your suggestion 

Richardson/SOR. <div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Xiaodong Liu, PhD<br>X: Computational Physics Division<br>Los Alamos National Laboratory<br>P.O. Box 1663, <br>Los Alamos, NM 87544<br>505-709-0534<br></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 18, 2020 at 10:34 AM Patrick Sanan <<a href="mailto:patrick.sanan@gmail.com">patrick.sanan@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 dir="ltr">I may not be interpreting the question correctly, but you can use SOR as the only preconditioner, e.g. <div><br></div><div>    -ksp_type gmres -pc_type sor</div><div><br></div><div>or you can use Richardson/SOR as your multigrid smoother, which might be what you're asking for</div><div><br></div><div><div>     -ksp_type gmres -pc_type mg -pc_mg_levels 3</div><div>    -mg_levels_ksp_type richardson -mg_levels_pc_type SOR</div></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Am Mi., 18. März 2020 um 18:19 Uhr schrieb Xiaodong Liu <<a href="mailto:xliu29@ncsu.edu" target="_blank">xliu29@ncsu.edu</a>>:<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 dir="ltr"><div>Hi, Petsc team, <br></div><div><br></div><div>I am doing something using multigrid as a preconditioner for GMRES. From the official case, </div><div><br></div><div>most use KSP (Chebyshev) and PC (SOR) for solving the preconditioning system for every level of multigrid except the coarse one. </div><div><br></div><div>-ksp_type gmres -pc_type mg -pc_mg_levels 3</div><div>-mg_levels_ksp_type Chebyshev -mg_levels_pc_type SOR</div><div><br></div><div> I would like to  use SOR for solving the preconditioning system directly. I can not find that Petsc supports non-KSP iterative solver for a linear system, right? </div><div><br></div><div>If I am wrong, please let me know.</div><div><br clear="all"><div><div dir="ltr"><div dir="ltr"><div>Xiaodong Liu, PhD<br>X: Computational Physics Division<br>Los Alamos National Laboratory<br>P.O. Box 1663, <br>Los Alamos, NM 87544<br>505-709-0534<br></div></div></div></div></div></div>
</blockquote></div>
</blockquote></div>