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
Attachment:
pgp0DsvnHOCRd.pgp
Description: PGP signature