<div dir="ltr">Dear Gromacs developers,<div>I am trying to hack Gromacs to include a new 4-body short-ranged potential which Gromacs has not implemented yet. My strategy is to define the potential as a new dihedral types. I will have 2 input parameters, so I expected that I will do something very similar to what Gromacs does for improper dihedral angle. Therefore by grepping for F_IDIHS and eNR_IMPROPER, I think I need to modify the following files:<div><div>src/gromacs/listed-forces/bonded.h</div><div>src/gromacs/listed-forces/bonded.cpp</div><div>src/gromacs/gmxlib/ifunc.c<br></div><div>src/gromacs/gmxlib/nrnb.c</div><div>src/gromacs/legacyheaders/types/nrnb.h</div><div>src/gromacs/topology/idef.h</div><div>src/gromacs/gmxpreprocess/topdirs.c<br></div><div>src/gromacs/gmxpreprocess/convparm.c</div><div>src/gromacs/gmxpreprocess/toppush.c<br></div><div>src/gromacs/fileio/tpxio.c</div><div><br></div><div>However, when I run test, I get the following fatal error at grompp step:</div><div><div>&quot;Incorrect number of parameters - found 2, expected 32621 or -26576307 for. &quot;</div></div><div><br></div><div>I think there is something wrong related to toppush.c but I am really clueless. The file is too complicated for me to understand. In your opinion, what is wrong with my modification? Besides the following file, is there any other files which I need to modify also?</div><div><br></div><div>I am looking forward to hearing from you.</div><div>Thank you very much in advance.</div><div>Best regards, </div><div><div><br></div></div><div><br></div><div><div class="gmail_signature"><div dir="ltr"><i>Lan Hoa</i></div></div></div>
</div></div></div>