PolyChordLite
PolyChordLite copied to clipboard
Pass custom mpi4py communicator
Hi, Thanks for this sampler! Here are two suggestions:
- it'd be nice to be able to provide, through the Python wrapper, the proper mpi4py communicator in case it is not COMM_WORLD (e.g. to have several runs in parallel).
- it may be worth avoiding MPI_ANY_TAG, e.g. in https://github.com/PolyChord/PolyChordLite/blob/e00386ba83766ea04365c055b5f17d6f62c0ee04/src/polychord/mpi_utils.F90#L492, which may conflict with potential MPI communication between e.g. the user-provided logposterior and dumper.
I have a branch that addresses the above suggestions I think, here: https://github.com/adematti/PolyChordLite/tree/mpi4py (I can file a PR if desired) Best
Hi @adematti. Both of these sound like a great idea. Feel free to file a PR.