<div dir="ltr"><span style="font-size:12.7272720336914px">Hi,</span><br style="font-size:12.7272720336914px"><br style="font-size:12.7272720336914px"><span style="font-size:12.7272720336914px">I would like to call the integrator function do_steep() from within do_md(). I modified some constraints parameters in t_inputrec before passing all the arguments supplied to do_md() by the user to do_steep(). This simple approach, however, does not work and GROMACS crashes complaining about a charge group moving too far. The coordinates seem to change by several nano meters over a single SD step. What is this happening? and Do you have any suggestions on how to implement this right?</span><br style="font-size:12.7272720336914px"><br style="font-size:12.7272720336914px"><span style="font-size:12.7272720336914px">Thanks</span><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><font face="georgia, serif"> <br></font><font color="#888888"><span style="font-size:13px;background-color:rgb(255,255,255)"><div style="font-family:arial,sans-serif"><b style="color:rgb(7,55,99);font-family:&#39;trebuchet ms&#39;,sans-serif;font-size:small">Andrew Abi-Mansour</b><font face="georgia, serif"><br></font></div></span></font></div></div></div></div>
</div>