Partition mesh on single process

I am coding up a parallel algorithm in which each processor has a copy of a mesh. Each processor then identifies a subset of the cells in this coarse mesh as its 'home domain.' The home domains should be distinct and cover the whole coarse mesh. I want to partition the mesh on each processor to identify home domains by leveraging the mesh partitioning already implemented in fenics. It is as if I want to parallelize a local mesh: I want to calculate which cells should be owned by each processor, but don't want to actually redistribute these cells. My understanding from asking this question in the FAQ forum is that this is not easily achievable in the current version of fenics, but is on the radar for further development. So I am raising this issue here so that this feature can hopefully be added to the code. Thanks.

Comments (2)