3395f156e75ae5b7a00d8e218073234828a7624c
[alexxy/gromacs.git] / src / kernel / topexcl.h
1 /*
2  * 
3  *                This source code is part of
4  * 
5  *                 G   R   O   M   A   C   S
6  * 
7  *          GROningen MAchine for Chemical Simulations
8  * 
9  *                        VERSION 3.2.0
10  * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
11  * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
12  * Copyright (c) 2001-2004, The GROMACS development team,
13  * check out http://www.gromacs.org for more information.
14
15  * This program is free software; you can redistribute it and/or
16  * modify it under the terms of the GNU General Public License
17  * as published by the Free Software Foundation; either version 2
18  * of the License, or (at your option) any later version.
19  * 
20  * If you want to redistribute modifications, please consider that
21  * scientific software is very special. Version control is crucial -
22  * bugs must be traceable. We will be happy to consider code for
23  * inclusion in the official distribution, but derived work must not
24  * be called official GROMACS. Details are found in the README & COPYING
25  * files - if they are missing, get the official version at www.gromacs.org.
26  * 
27  * To help us fund GROMACS development, we humbly ask that you cite
28  * the papers on the package - you can find them in the top README file.
29  * 
30  * For more info, check our website at http://www.gromacs.org
31  * 
32  * And Hey:
33  * Gallium Rubidium Oxygen Manganese Argon Carbon Silicon
34  */
35
36 #ifndef _topexcl_h
37 #define _topexcl_h
38 #include "visibility.h"
39 #include "topio.h"
40
41 typedef struct {
42   int nr;               /* nr atoms (0 <= i < nr) (atoms->nr)           */
43   int nrex;             /* with nrex lists of neighbours                */
44                         /* respectively containing zeroth, first        */
45                         /* second etc. neigbours (0 <= nre < nrex)      */
46   int **nrexcl;         /* with (0 <= nrx < nrexcl[i][nre]) neigbours   */ 
47                         /* per list stored in one 2d array of lists     */ 
48   int ***a;             /* like this: a[i][nre][nrx]                    */
49 } t_nextnb;
50
51 GMX_LIBGMXPREPROCESS_EXPORT
52 extern void init_nnb(t_nextnb *nnb, int nr, int nrex);
53 /* Initiate the arrays for nnb (see above) */
54
55 GMX_LIBGMXPREPROCESS_EXPORT
56 extern void done_nnb(t_nextnb *nnb);
57 /* Cleanup the nnb struct */
58
59 #ifdef DEBUG_NNB
60 #define print_nnb(nnb, s) __print_nnb(nnb, s)
61 extern void print_nnb(t_nextnb *nnb, char *s);
62 /* Print the nnb struct */
63 #else
64 #define print_nnb(nnb, s)
65 #endif
66
67 GMX_LIBGMXPREPROCESS_EXPORT
68 extern void gen_nnb(t_nextnb *nnb,t_params plist[]);
69 /* Generate a t_nextnb structure from bond information. 
70  * With the structure you can either generate exclusions
71  * or generate angles and dihedrals. The structure must be
72  * initiated using init_nnb.
73  */
74
75 extern void nnb2excl (t_nextnb *nnb, t_blocka *excl);
76 /* generate exclusions from nnb */
77
78 extern void generate_excl (int nrexcl, int nratoms,
79                            t_params plist[],t_blocka *excl);
80 /* Generate an exclusion block from bonds and constraints in
81  * plist.
82  */
83
84 #endif  /* _topexcl_h */