mesh ordering and partition

[Gong Ding]

Hi,
I am studying the parallel programming, using libmesh/petsc as an excellent 
example.

I have some questions about the partition and mesh ordering.



It seems libmesh does not reorder the mesh nodes. It only calls metis to 
partition the mesh, and

uses original node order to build the matrix. I wonder if a bad mesh 
ordering may cause low efficency

of ILU preconditioner. However, If I did RCM ordering to the mesh, the 
node's order may conflict with

contiguous index set int the subdomain partitioned by metis. How should I 
balance the ordering (to reduce filling)

and partition (to reduce communication)? any good ideas?



Regards,

Gong Ding