> On Feb 7, 2017, at 2:49 AM, Michel Kern <[email protected]> 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
>
>
>
