/usr/include/gromacs/topology/topology.h is in libgromacs-dev 2018.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 | /*
* This file is part of the GROMACS molecular simulation package.
*
* Copyright (c) 1991-2000, University of Groningen, The Netherlands.
* Copyright (c) 2001-2004, The GROMACS development team.
* Copyright (c) 2011,2014,2015,2016, by the GROMACS development team, led by
* Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
* and including many others, as listed in the AUTHORS file in the
* top-level source directory and at http://www.gromacs.org.
*
* GROMACS is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public License
* as published by the Free Software Foundation; either version 2.1
* of the License, or (at your option) any later version.
*
* GROMACS is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with GROMACS; if not, see
* http://www.gnu.org/licenses, or write to the Free Software Foundation,
* Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
*
* If you want to redistribute modifications to GROMACS, please
* consider that scientific software is very special. Version
* control is crucial - bugs must be traceable. We will be happy to
* consider code for inclusion in the official distribution, but
* derived work must not be called official GROMACS. Details are found
* in the README & COPYING files - if they are missing, get the
* official version at http://www.gromacs.org.
*
* To help us fund GROMACS development, we humbly ask that you cite
* the research papers on the package. Check out http://www.gromacs.org.
*/
#ifndef GMX_TOPOLOGY_TOPOLOGY_H
#define GMX_TOPOLOGY_TOPOLOGY_H
#include <cstdio>
#include "gromacs/math/vectypes.h"
#include "gromacs/topology/atoms.h"
#include "gromacs/topology/block.h"
#include "gromacs/topology/idef.h"
#include "gromacs/topology/symtab.h"
enum {
egcTC, egcENER, egcACC, egcFREEZE,
egcUser1, egcUser2, egcVCM, egcCompressedX,
egcORFIT, egcQMMM,
egcNR
};
/* Names corresponding to groups */
extern const char *gtypes[egcNR+1];
typedef struct gmx_moltype_t
{
char **name; /* Name of the molecule type */
t_atoms atoms; /* The atoms in this molecule */
t_ilist ilist[F_NRE]; /* Interaction list with local indices */
t_block cgs; /* The charge groups */
t_blocka excls; /* The exclusions */
} gmx_moltype_t;
/*! \brief Block of molecules of the same type, used in gmx_mtop_t */
typedef struct gmx_molblock_t
{
int type; /**< The molecule type index in mtop.moltype */
int nmol; /**< The number of molecules in this block */
int nposres_xA; /**< The number of posres coords for top A */
rvec *posres_xA; /**< Position restraint coordinates for top A */
int nposres_xB; /**< The number of posres coords for top B */
rvec *posres_xB; /**< Position restraint coordinates for top B */
/* Convenience information, derived from other gmx_mtop_t contents */
int natoms_mol; /**< The number of atoms in one molecule */
int globalAtomStart; /**< Global atom index of the first atom in the block */
int globalAtomEnd; /**< Global atom index + 1 of the last atom in the block */
int globalResidueStart; /**< Global residue index of the first residue in the block */
int residueNumberStart; /**< Residue numbers start from this value if the number of residues per molecule is <= maxres_renum */
} gmx_molblock_t;
typedef struct gmx_groups_t
{
t_grps grps[egcNR]; /* Groups of things */
int ngrpname; /* Number of groupnames */
char ***grpname; /* Names of the groups */
int ngrpnr[egcNR];
unsigned char *grpnr[egcNR]; /* Group numbers or NULL */
} gmx_groups_t;
/* This macro gives the group number of group type egc for atom i.
* This macro is useful, since the grpnr pointers are NULL
* for group types that have all entries 0.
*/
#define ggrpnr(groups, egc, i) ((groups)->grpnr[egc] ? (groups)->grpnr[egc][i] : 0)
/* The global, complete system topology struct, based on molecule types.
This structure should contain no data that is O(natoms) in memory. */
typedef struct gmx_mtop_t
{
char **name; /* Name of the topology */
gmx_ffparams_t ffparams;
int nmoltype;
gmx_moltype_t *moltype;
int nmolblock;
gmx_molblock_t *molblock;
gmx_bool bIntermolecularInteractions; /* Are there intermolecular
* interactions? */
t_ilist *intermolecular_ilist; /* List of intermolecular interactions
* using system wide atom indices,
* either NULL or size F_NRE */
int natoms;
int maxres_renum; /* Parameter for residue numbering */
int maxresnr; /* The maximum residue number in moltype */
t_atomtypes atomtypes; /* Atomtype properties */
t_block mols; /* The molecules */
gmx_groups_t groups;
t_symtab symtab; /* The symbol table */
} gmx_mtop_t;
/* The mdrun node-local topology struct, completely written out */
typedef struct gmx_localtop_t
{
t_idef idef; /* The interaction function definition */
t_atomtypes atomtypes; /* Atomtype properties */
t_block cgs; /* The charge groups */
t_blocka excls; /* The exclusions */
} gmx_localtop_t;
/* The old topology struct, completely written out, used in analysis tools */
typedef struct t_topology
{
char **name; /* Name of the topology */
t_idef idef; /* The interaction function definition */
t_atoms atoms; /* The atoms */
t_atomtypes atomtypes; /* Atomtype properties */
t_block cgs; /* The charge groups */
t_block mols; /* The molecules */
gmx_bool bIntermolecularInteractions; /* Inter.mol. int. ? */
t_blocka excls; /* The exclusions */
t_symtab symtab; /* The symbol table */
} t_topology;
void init_mtop(gmx_mtop_t *mtop);
void init_top(t_topology *top);
void done_moltype(gmx_moltype_t *molt);
void done_molblock(gmx_molblock_t *molb);
void done_gmx_groups_t(gmx_groups_t *g);
void done_mtop(gmx_mtop_t *mtop);
void done_top(t_topology *top);
// Frees both t_topology and gmx_mtop_t when the former has been created from
// the latter.
void done_top_mtop(t_topology *top, gmx_mtop_t *mtop);
bool gmx_mtop_has_masses(const gmx_mtop_t *mtop);
bool gmx_mtop_has_charges(const gmx_mtop_t *mtop);
bool gmx_mtop_has_atomtypes(const gmx_mtop_t *mtop);
bool gmx_mtop_has_pdbinfo(const gmx_mtop_t *mtop);
void pr_mtop(FILE *fp, int indent, const char *title, const gmx_mtop_t *mtop,
gmx_bool bShowNumbers, gmx_bool bShowParameters);
void pr_top(FILE *fp, int indent, const char *title, const t_topology *top,
gmx_bool bShowNumbers, gmx_bool bShowParameters);
void cmp_top(FILE *fp, const t_topology *t1, const t_topology *t2, real ftol, real abstol);
void cmp_groups(FILE *fp, const gmx_groups_t *g0, const gmx_groups_t *g1,
int natoms0, int natoms1);
#endif
|