<div dir="ltr">Thanks a lot Pierre, I use PCGAMG as I type this answer.  Not only that it works, but converges much faster than PCASM :-)<div><br></div><div>    Have a great day,</div><div><br></div><div>    Bojan</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 20, 2023 at 8:07 AM Pierre Jolivet <<a href="mailto:pierre.jolivet@lip6.fr">pierre.jolivet@lip6.fr</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"><br>
> On 20 Apr 2023, at 7:53 AM, Bojan Niceno <<a href="mailto:bojan.niceno.scientist@gmail.com" target="_blank">bojan.niceno.scientist@gmail.com</a>> wrote:<br>
> <br>
> Dear all,<br>
> <br>
> <br>
> I am solving a Laplace equation with finite volume method on an unstructured grid with a Fortran code I have developed, and PETSc 3.19 library.<br>
> <br>
> I first used cg solver with asm preconditioner, which converges nicely when executed sequentially, but fails for MPI parallel version.  I believed that there must be an error in which I set up and assemble parallel matrices for PETSc, but I soon noticed that if I use bicg with asm, everything works fine, the parallel bicg/asm shows almost the same convergence as sequential version.<br>
<br>
KSPCG requires a symmetric PC.<br>
By default, PCASMType is PC_ASM_RESTRICT, which yields a non-symmetric preconditioner.<br>
With a single process, this does not matter, but with more than one process, it does.<br>
If you switch to -pc_asm_type basic, KSPCG should converge.<br>
That being said, for Laplace equation, there are much faster alternatives than PCASM, e.g., PCGAMG.<br>
<br>
Thanks,<br>
Pierre<br>
<br>
> I could carry on with bicg, but I am still worried a little bit.  To my knowledge of Krylov solvers, which is admittedly basic since I am a physicist only using linear algebra, the convergence of bicg should be very similar to that of cg when symmetric systems are solved.  When I run my cases sequentially, I see that's indeed the case.  But in parallel, bicg converges and cg fails.<br>
> <br>
> Do you see the above issues as an anomaly and if so, could you advise how to search for a cause?  <br>
> <br>
> <br>
>     Kind regards,<br>
> <br>
> <br>
>     Bojan<br>
> <br>
> <br>
> <br>
> <br>
> <br>
<br>
</blockquote></div>