How to run mdrun with mpiexec?

[dpadula85]

Hello,

I went a bit through past issues and the code, but I can't seem to figure out how to change the configuration file so that mdrun gets called with an mpiexec -n prefix.

Can someone help me?

Thanks,
D.

[orbeckst]

You can use gromacs.run.MDrunner. This code is not used that often and the implementation and documentation is not great, I am afraid. But perhaps it's useful for you:

You'll have to derive your own class from MDrunner. There isn't a lot of documentation so you have to look at the source code. There are some examples https://gromacswrapper.readthedocs.io/en/latest/gromacs/blocks/simulation.html#example-implementations (again, look at the source), for example MDrunnerOpenMP

GromacsWrapper/gromacs/run.py

Lines 283 to 289 in 76ac509

	class MDrunnerOpenMP(MDrunner):
	"""Manage running :program:`mdrun` as an OpenMP_ multiprocessor job.
	.. _OpenMP: http://openmp.org/wp/
	"""
	mdrun = ("mdrun_openmp", "gmx_openmp mdrun")
	mpiexec = "mpiexec"

You can provide your full path to mpiexec in

GromacsWrapper/gromacs/run.py

Line 289 in 76ac509

mpiexec = "mpiexec"

Put the names of your mpi-enabled mdrun binaries in

GromacsWrapper/gromacs/run.py

Line 288 in 76ac509

mdrun = ("mdrun_openmp", "gmx_openmp mdrun")

If you need to run some code before or after you can use the prehook() and posthook() methods as show in MDrunnerMpich2Smpd

GromacsWrapper/gromacs/run.py

Line 291 in 76ac509

class MDrunnerMpich2Smpd(MDrunner):

— but you might not need such a complicated setup.

Maybe define a class in your own code like this

import gromacs.run
class MDrunnerMPI(gromacs.run.MDrunner):
    """Manage running :program:`mdrun` as an MPI multiprocessor job.
    """
    mdrun = ("mdrun_mpi", "gmx_mpi mdrun")
    mpiexec = "/opt/local/bin/mpiexec"

and then use it as

mdrun_mpi = MDrunnerMPI(s="md.tpr", deffnm="md")     # add mdrun commands here
rc = mdrun_mpi.run(ncores=16)

The above basic MDrunnerMPI only supports running mpiexec -n ncores gmx mdrun ..., i.e., only the -n ncores arguments for mpiexec are supported. If you need more then you need write your own mpicommand() method, which you would add to your own MDrunnerMPI class. If this is the case, I suggest you just ask here if you need help.

[dpadula85]

Thanks,

I had figured out the classes to use to override the default one! Your examples are of great help, and maybe this can be inserted in the docs!

Thanks again,
D.

[orbeckst]

Glad that this was helpful.

I will re-open your issue as a reminder for documentation.