mesh ordering and partition

Jed Brown jed at
Fri May 16 12:58:58 CDT 2008

On Sat 2008-05-17 00:47, Gong Ding wrote:
> I am developing parallel code for 3D semiconductor device simulation.
> From the experience of 2D code, the GMRES solver with ILU works well (the 
> matrix is asymmetric.)
> As a result, I'd like to use GMRES+ILU again for 3D,  in parallel.
> Does   -pc_type ilu -pc_factor_mat_ordering_type rcm still work?
> Since the parallel martrix requires continuous index in subdomain, the 
> matrix ordering seems troublesome.
> maybe only a local ordering can be done... Am I right?

For parallel ILU, you could try -pc_type hypre -pc_hypre_type euclid.
Unfortunately, ILU requires a lot of communication so the parallel scaling tends
to be poor.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 197 bytes
Desc: not available
URL: <>

More information about the petsc-users mailing list