Becksteinlab / GromacsWrapper

GromacsWrapper wraps system calls to GROMACS tools into thin Python classes (GROMACS 4.6.5 - 2024 supported).

Home Page:https://gromacswrapper.readthedocs.org

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

setup.py:energy_minimize ignores mdp parameter

luancarvalhomartins opened this issue · comments

In setup.py, line 662:

mdp = deffnm+'.mdp'

This overwrites user-supplied mdp parameter.

Yes – if you want to provide a PR I am happy to merge it. I won't be able to work on it during this week.

Thanks for reporting!

I can do so, but only between Friday and Sunday. Would you rather prefer setup.py:energy_minimize to try to alter user-supplied mdp file through cbook.py:editmdp or to use it as is? I tend to think that the former allows more flexibility.

You are welcome.

I understood this issue that a user-supplied mdp parameter in function gromacs.setup.energy_minimize(mdp="custom.mdp") is ignored. In fact, this is not true and it works just fine with a supplied mdp file. The code is not very clear because the variable mdp changes meaning half-way through (sorry). I added a comment to the code and a test for a custom mdp file in PR #128.

(It is true that a custom em/em.mdp file would be overwritten and this is not checked for.)