4 * This source code is part of
8 * GROningen MAchine for Chemical Simulations
11 * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
12 * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
13 * Copyright (c) 2001-2004, The GROMACS development team,
14 * check out http://www.gromacs.org for more information.
16 * This program is free software; you can redistribute it and/or
17 * modify it under the terms of the GNU General Public License
18 * as published by the Free Software Foundation; either version 2
19 * of the License, or (at your option) any later version.
21 * If you want to redistribute modifications, please consider that
22 * scientific software is very special. Version control is crucial -
23 * bugs must be traceable. We will be happy to consider code for
24 * inclusion in the official distribution, but derived work must not
25 * be called official GROMACS. Details are found in the README & COPYING
26 * files - if they are missing, get the official version at www.gromacs.org.
28 * To help us fund GROMACS development, we humbly ask that you cite
29 * the papers on the package - you can find them in the top README file.
31 * For more info, check our website at http://www.gromacs.org
34 * Gromacs Runs On Most of All Computer Systems
50 extern void check_index(char *gname
,int n
,atom_id index
[],
51 char *traj
,int natoms
);
52 /* Checks if any index is smaller than zero or larger than natoms,
53 * if so a fatal_error is given with the gname (if gname=NULL, "Index" is used)
54 * and traj (if traj=NULL, "the trajectory" is used).
57 t_block
*init_index(char *gfile
, char ***grpname
);
58 /* Lower level routine than the next */
60 void rd_index(char *statfile
,int ngrps
,int isize
[],
61 atom_id
*index
[],char *grpnames
[]);
62 /* Assume the group file is generated, so the
63 * format need not be user-friendly. The format is:
64 * nr of groups, total nr of atoms
65 * for each group: name nr of element, elements
66 * The function opens a file, reads ngrps groups, puts the
67 * sizes in isize, the atom_id s in index and the names of
68 * the groups in grpnames.
70 * It is also assumed, that when ngrps groups are requested
71 * memory has been allocated for ngrps index arrays, and that
72 * the dimension of the isize and grpnames arrays are ngrps.
75 void rd_index_nrs(char *statfile
,int ngrps
,int isize
[],
76 atom_id
*index
[],char *grpnames
[],int grpnr
[]);
77 /* the same but also reads the number of the selected group*/
79 void get_index(t_atoms
*atoms
, char *fnm
, int ngrps
,
80 int isize
[], atom_id
*index
[],char *grpnames
[]);
81 /* Does the same as rd_index, but if the fnm pointer is NULL it
82 * will not read from fnm, but it will make default index groups
83 * for the atoms in *atoms.
91 extern t_aa_names
*get_aa_names(void);
92 /* Read the database in aminoacids.dat */
94 extern bool is_protein(t_aa_names
*aan
,char *resnm
);
95 /* gives true if resnm occurs in aminoacids.dat */
97 extern void done_aa_names(t_aa_names
**aan
);
98 /* Free memory. Pass address of the pointer youget from get_aa_names */
100 extern t_block
*new_block(void);
101 /* allocate new block */
103 extern void write_index(char *outf
, t_block
*b
,char **gnames
);
104 /* Writes index blocks to outf (writes an indexfile) */
106 void add_grp(t_block
*b
,char ***gnames
,int nra
,atom_id a
[],const char *name
);
107 /* Ads group a with name name to block b and namelist gnames */
109 extern void analyse(t_atoms
*atoms
,t_block
*gb
,char ***gn
,
110 bool bASK
,bool bVerb
);
111 /* Makes index groups gb with names gn for atoms in atoms.
112 * bASK=FALSE gives default groups.
115 extern int find_group(char s
[], int ngrps
, char **grpname
);
122 #endif /* _index_h */