mesh ordering and partition

[Jed Brown]

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.

Jed
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 197 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080516/b3b6f9c6/attachment.pgp>