<div dir="ltr">PS: Perhaps we should simply suggest to use no DD and OpenMP parallelization only?</div><div class="gmail_extra"><br clear="all"><div>--<br>Szilárd</div>
<br><br><div class="gmail_quote">On Fri, May 2, 2014 at 5:57 PM, Szilárd Páll <span dir="ltr"><<a href="mailto:pall.szilard@gmail.com" target="_blank">pall.szilard@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>I see. So that means that SHAKE is even more useless and DD will work only if the domains is empty?</div></div><div class="gmail_extra"><br clear="all"><div>--<br>Szilárd</div>
<br><br><div class="gmail_quote"><div><div class="h5">On Fri, May 2, 2014 at 5:37 PM, <a href="mailto:hess@kth.se" target="_blank">hess@kth.se</a> <span dir="ltr"><<a href="mailto:hess@kth.se" target="_blank">hess@kth.se</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
<div style="font-size:12pt;font-family:Calibri,sans-serif">Hi,<div><br></div><div>Since with the Verlet scheme charge groups are single atoms, this error is correct. But some more explanation could be useful.</div><div><br>
</div><div>Cheers,</div><div><br></div><div>Berk</div><div><div><br><div>----- Reply message -----<br>From: "Szilárd Páll" <<a href="mailto:pall.szilard@gmail.com" target="_blank">pall.szilard@gmail.com</a>><br>
To: "Discussion list for GROMACS development" <<a href="mailto:gmx-developers@gromacs.org" target="_blank">gmx-developers@gromacs.org</a>><br>Subject: [gmx-developers] strange SHAKE + Verlet scheme warning?<br>
Date: Fri, May 2, 2014 17:25</div></div></div></div><div><div><br><div dir="ltr"><div><span style="font-family:arial,sans-serif;font-size:13px">Hi,</span><div style="font-family:arial,sans-serif;font-size:13px">
<br></div><div style="font-family:arial,sans-serif;font-size:13px">I've come across a somewhat strange issue and I thought I'd ask before filing a redmine - especially as I'm not sure whether this is only a documentation/code self-documentation bug or a buggy behavior too.</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">I'm running 5.0 RC1 (git), Verlet scheme with a small water box of ~500 molecules and -DFLEXIBLE with SHAKE constraints. </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Trying to start a DD run I get this error:</div><div style="font-family:arial,sans-serif;font-size:13px">
<br></div><div style="font-family:arial,sans-serif;font-size:13px"><div><div>-------------------------------------------------------</div><div>Program mdrun_d, VERSION 5.0-rc1-dev-20140430-81378bd</div><div>Source code file: /home/pszilard/projects/gromacs/gromacs-5.0/src/gromacs/mdlib/constr.c, line: 1266</div>
<div><br></div><div>Fatal error:</div><div>SHAKE is not supported with domain decomposition and constraint that cross charge group boundaries, use LINCS</div><div>For more information and tips for troubleshooting, please check the GROMACS</div>
<div>website at <a href="http://www.gromacs.org/Documentation/Errors" target="_blank">http://www.gromacs.org/Documentation/Errors</a></div><div>-------------------------------------------------------</div></div><div><br>
</div>
</div><div style="font-family:arial,sans-serif;font-size:13px">Now, this is inconsistent at least with the fact that charge groups are removed by grompp with the Verlet scheme. The documentation does point out that "SHAKE does not support constraints between atoms on different nodes", though.</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Looking a the code, the bInterCGcons that triggers it still refers to charge groups, so there is definitely room for improvement both in error message, variable naming, as well as documentation.</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Cheers,</div><div style="font-family:arial,sans-serif;font-size:13px">--<br>Szilárd</div>
</div>
</div>
</div></div><br></div></div><span class="HOEnZb"><font color="#888888">--<br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List</a> before posting!<br>
<br>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request@gromacs.org</a>.<br>
</font></span></blockquote></div><br></div>
</blockquote></div><br></div>