[Libmesh-devel] Possible bug in Partitioner

[Vijay S. M. <vi...@gm...>]

Hi,

In a simple exercise of building FE from discretizing a cube
(MeshGeneration::build_cube), I have encountered assert failures in
partitioner.C when running in more than 1 processors;

The exact error message is shown below

src/partitioning/partitioner.C:189: static void
Partitioner::partition_unpartitioned_elements(MeshBase&, unsigned
int): Assertion `cnt < global_indices.size()' failed.

I think part of the reason is that the call to MeshCommunication
probably is strongly dependent on LibHilbert to be present. I
currently do not have Hilbert as part of my libMesh build and it is
possible that the global_indices are not populated leading to the
assertion failure.

Here's the call :

MeshCommunication().find_global_indices
(MeshTools::bounding_box(mesh), it, end, global_indices);

I currently do not have access to my libMesh build. When i get time, I
will try to build it again with Hilbert and see if it resolves the
issue. Just thought i should let the developers know about this and if
there is something more to it other than just including Hilbert
library.

Cheers,
Vijay