[petsc-users] DIVERGED_INDEFINITE_PC in algebraic multigrid

Michele Rosso mrosso at uci.edu
Mon May 26 11:20:25 CDT 2014


Mark,

thank you for your input and sorry my late reply: I saw your email only now.
By setting up the solver each time step you mean re-defining the KSP 
context every time? Why should this help?
I will definitely try that as well as the hypre solution and report back.
Again, thank you.

Michele

On 05/22/2014 09:34 AM, Mark Adams wrote:
> If the solver is degrading as the coefficients change, and I would 
> assume get more nasty, you can try deleting the solver at each time 
> step.  This will be about 2x more expensive, because it does the setup 
> each solve, but it might fix your problem.
>
> You also might try:
>
> -pc_type hypre
> -pc_hypre_type boomeramg
>
>
>
>
> On Mon, May 19, 2014 at 6:49 PM, Jed Brown <jed at jedbrown.org 
> <mailto:jed at jedbrown.org>> wrote:
>
>     Michele Rosso <mrosso at uci.edu <mailto:mrosso at uci.edu>> writes:
>
>     > Jed,
>     >
>     > thank you very much!
>     > I will try with ///-mg_levels_ksp_type chebyshev -mg_levels_pc_type
>     > sor/   and report back.
>     > Yes, I removed the nullspace from both the system matrix and the
>     rhs.
>     > Is there a way to have something similar to Dendy's multigrid or the
>     > deflated conjugate gradient method with PETSc?
>
>     Dendy's MG needs geometry.  The algorithm to produce the interpolation
>     operators is not terribly complicated so it could be done, though DMDA
>     support for cell-centered is a somewhat awkward.  "Deflated CG"
>     can mean
>     lots of things so you'll have to be more precise.  (Most everything in
>     the "deflation" world has a clear analogue in the MG world, but the
>     deflation community doesn't have a precise language to talk about
>     their
>     methods so you always have to read the paper carefully to find out if
>     it's completely standard or if there is something new.)
>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140526/82644d03/attachment.html>


More information about the petsc-users mailing list