[petsc-users] Replacement for Euclid ?

[Barry Smith]

> On Feb 7, 2017, at 2:49 AM, Michel Kern <michel.kern at inria.fr> wrote:
> 
> I have a code that uses euclid for ILU preconditioning (as part of a CPR-AMG preconditioner for a mulitphase flow solver). As noted in the release notes, euclid is no longer supported starting from 3.6. What is the recommended replacement ? Is block Jacobi plus a local solver an alternative ? 

   I would use block Jacobi with Gauss-Seidel as the local solver first. Does that not work due to zero pivots?

    Then you can use PCFIELDSPLIT with the option -pc_fieldsplit_detect_saddle_point 

    or you can use block Jacobi with ILU(0) as the local solver.

    There are many combinations to try.

  Barry

> 
> Thanks
> Michel
> 
> -- 
> Michel Kern
> Centre de Recherche Inria de Paris,    Maison de la Simulation
> 2 rue Simone Iff, bureau A413,            Digiteo Labs, bât 565, CEA Saclay 
> 75012 Paris                                          91191 Gif-sur-Yvette cedex 
> tel. : +33 (0)1 80 49 42 36,                  tel. : +33 (0)1 69 08 58 32            
> 
>