[petsc-users] Laplace Equation preconditioner
Barry Smith
bsmith at petsc.dev
Wed Oct 12 13:32:43 CDT 2022
What KSP are you using? DIVERGED_BREAKDOWN is very rare for KSPGMRES. If you are using one of its lesser cousins like bcgs you might consider switching to bcgsl or gmres.
I assume because of boundary conditions or the discretization you do not have symmetric positive definite and thus cannot use CG.
Barry
> On Oct 12, 2022, at 2:00 PM, Matthew Knepley <knepley at gmail.com> wrote:
>
> On Wed, Oct 12, 2022 at 1:38 PM Alfredo J Duarte Gomez <aduarteg at utexas.edu <mailto:aduarteg at utexas.edu>> wrote:
> Good morning PETSC users,
>
> I have a current solver that requires the solution of a Laplace equation to be reused for all future time steps.
>
> The configuration is axisymmetric with Dirichlet BCs at the top and bottom boundaries, and Zero Neumman conditions at the axis and far field. The grid is curvilinear and structured.
>
> So far I have been using PCHYPRE boomeramg as the preconditioner, which often works well, but I have also experienced DIVERGED_BREAKDOWN on many occasions. When I use direct solver PCMUMPS it always produces a satisfactory answer, which gives me confidence that the solution exists for the given grid in all these configurations where boomeramg fails.
>
> I do not know why Hypre is breaking down. Did you try ML or GAMG? They are easier to diagnose I think.
>
> Thanks,
>
> Matt
>
> I am looking for recommendations on other preconditioners to try in this problem that can produce a solution faster than PCMUMPS, or recommendations on which parameters to adjust on boomeramg.
>
> Thank you,
>
>
> --
> Alfredo Duarte
> Graduate Research Assistant
> The University of Texas at Austin
>
>
> --
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221012/ebcf1004/attachment.html>
More information about the petsc-users
mailing list