[gmx-developers] How to obtain the element of an atom?

[Berk Hess]

Hi,

Why is the element number relevant for constraints? The element number 
is never used in MD calculations.

Cheers,

Berk

On 2/11/21 3:29 PM, Lorién López Villellas wrote:
> Hi all.
>
> We are working on a new parallel algorithm for imposing constraints 
> based on SHAKE. We are currently trying to integrate it into GROMACS.
>
> We use the iatoms array (topology of triplets of constraint type, and 
> indices of the two atoms that compose a bond) just like SHAKE. We need 
> a way to obtain the element (Carbon, Nitrogen, Hidrogen...) of each of 
> the atom indexes in iatoms. As far as we have researched, there is no 
> way to access that data from src/gromacs/mdlib/constr. Where is that 
> information stored?
>
> If there is no way to obtain an atom element, we could also adapt our 
> code to work with similar information if it exists.
>
> We really appreciate any help you can provide.
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20210211/838a514b/attachment.html>