inducer/pymetis

Wrap partMeshNodal

stoiver opened this issue · 1 comments

Is your feature request related to a problem? Please describe.
@inducer thanks for this great work of wrapping metis. I am developer of the flood and tsunami package https://github.com/anuga-community/anuga_core and we use metis (indeed your pymetis wrapper) to partition our triangular mesh for mpi runs.

In the past for our python 2 version we used metis 4 and our own wrapper of the partMeshNodal procedure. When we moved to python 3 we decided to upgrade to using metis 5. Instead of writing our own wrapper we decided to use your wrapper, and your pymetis.part_graph procedure. Worked a charm.

But on closer examination the quality of our partitions (using pymetis.part_graph) and the scaling of our parallel code, it seems that the partitions produced by pymetis.part_graph are a little inferior to the partitions created by partMeshNodal. Our communication regions seem to be about 20% larger with pymetis.part_graph

Describe the solution you'd like
So I would like to request that you also add a wrapper to partMeshNodal

Describe alternatives you've considered
I did have a very short look at your wrapper code. I would be willing to try to code up the extension myself, but would appreciate any help from you.

I'd be happy to consider a PR.

but would appreciate any help from you.

What type of help are you looking for?