Hello all,

It is possible to have a code snippet on how and when to call the parmetis package to partition the a mesh from a libmesh application code.

I've done like this after building a mesh object:

    // Partition the mesh with ParMetis package.
    ParmetisPartitioner pmetis;
    pmetis.partition(mesh, 4); // Partition the mesh on 4 processors
 
    // Print information about the mesh to the screen.
    mesh.print_info();

I get the following output which shows that I have only 1 subdomain:
 
Mesh Information:
  mesh_dimension()=3
  spatial_dimension()=3
  n_nodes()=29791
  n_elem()=3375
   n_local_elem()=843
   n_active_elem()=3375
  n_subdomains()=1
  n_processors()=4
  processor_id()=0


So, what's wrong ?

Thanks !
-as