<div dir="ltr">Hi.<br><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">One the possible solutions (there are atomic numbers stored in the gmx_mtop_t)</blockquote><div> </div><div>I think we can work with that. Thanks!<br><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">But that would make the algorithm work only for proteins. Why not do the numbering based on bond connectivity?</blockquote><br>We have two versions of the algorithm, a sequential one that works with any molecule and a parallel one that, as you say, only works with proteins. The parallelization of the algorithm is based on the structure of the peptide chain. The numbering we use is extremely specific; each amino acid bond must always have the same index. The reordering algorithm needs to identify each of the molecule bonds precisely. We have not found a way to achieve that without using the additional information that the atoms' element provides.<br><br>Regards,<br>Lorién<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue, 11 feb 2021 a las 17:30, Berk Hess (<<a href="mailto:hess@kth.se" target="_blank">hess@kth.se</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div>On 2/11/21 4:22 PM, Lorién López
Villellas wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">Hi.<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Why
is the element number relevant for constraints? The element
number is never used in MD calculations.</blockquote>
<br>
We have performed a heavy code optimization based on the
structural patterns of the peptide chains. In order to apply
this optimization, we need a very specific bond numbering. To
get to this bond numbering, we use a bond reordering phase. The
reordering algorithm needs to know the elements of the
molecule's atoms to locate itself.<br>
<br>
</div>
</blockquote>
But that would make the algorithm work only for proteins. Why not do
the numbering based on bond connectivity?<br>
<br>
Cheers,<br>
<br>
Berk<br>
<br>
</div>
-- <br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" rel="noreferrer" 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" rel="noreferrer" 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" rel="noreferrer" 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>.</blockquote></div>