1 |
mmeineke |
377 |
#ifndef _MOLECULE_H_ |
2 |
|
|
#define _MOLECULE_H_ |
3 |
|
|
|
4 |
|
|
#include "Atom.hpp" |
5 |
|
|
#include "SRI.hpp" |
6 |
mmeineke |
407 |
#include "MoleculeStamp.hpp" |
7 |
gezelter |
416 |
#include "Exclude.hpp" |
8 |
mmeineke |
377 |
|
9 |
mmeineke |
407 |
typedef struct{ |
10 |
|
|
|
11 |
|
|
int stampID; // the ID in the BASS component stamp array |
12 |
|
|
int nAtoms; // the number of atoms in the molecule |
13 |
|
|
int nBonds; // ... .. .. . .bonds .. .. . . . . |
14 |
|
|
int nBends; // . . . . .. . .bends . . . . .. . |
15 |
gezelter |
416 |
int nTorsions; // .. . . .. . . torsions . . .. . . |
16 |
|
|
int nOriented; // .. . . . .. . oriented atoms . . . |
17 |
mmeineke |
425 |
int nExcludes; // . . .. . . . .exclude pairs.. . . . |
18 |
gezelter |
446 |
|
19 |
mmeineke |
412 |
Atom** myAtoms; // the array of atoms |
20 |
|
|
Exclude** myExcludes;// the array of Excluded pairs |
21 |
mmeineke |
407 |
Bond** myBonds; // arrays of all the short range interactions |
22 |
|
|
Bend** myBends; |
23 |
|
|
Torsion** myTorsions; |
24 |
|
|
|
25 |
|
|
|
26 |
|
|
} molInit; |
27 |
|
|
|
28 |
mmeineke |
377 |
class Molecule{ |
29 |
|
|
|
30 |
|
|
public: |
31 |
|
|
|
32 |
mmeineke |
407 |
Molecule( void ); |
33 |
|
|
~Molecule( void ); |
34 |
mmeineke |
377 |
|
35 |
mmeineke |
407 |
void initialize( molInit &theInit ); |
36 |
chuckv |
438 |
|
37 |
|
|
void setMyIndex( int theIndex ){ myIndex = theIndex;} |
38 |
mmeineke |
423 |
|
39 |
|
|
int getNAtoms ( void ) {return nAtoms;} |
40 |
|
|
int getNBonds ( void ) {return nBonds;} |
41 |
|
|
int getNBends ( void ) {return nBends;} |
42 |
|
|
int getNTorsions( void ) {return nTorsions;} |
43 |
|
|
int getNOriented( void ) {return nOriented;} |
44 |
|
|
int getNExcludes( void ) {return nExcludes;} |
45 |
|
|
int getNMembers ( void ) {return nMembers;} |
46 |
|
|
int getStampID ( void ) {return stampID;} |
47 |
mmeineke |
407 |
|
48 |
mmeineke |
423 |
Atom** getMyAtoms ( void ) {return myAtoms;} |
49 |
|
|
Bond** getMyBonds ( void ) {return myBonds;} |
50 |
|
|
Bend** getMyBends ( void ) {return myBends;} |
51 |
mmeineke |
424 |
Torsion** getMyTorsions( void ) {return myTorsions;} |
52 |
|
|
Exclude** getMyExcludes( void ) {return myExcludes;} |
53 |
mmeineke |
377 |
|
54 |
mmeineke |
423 |
void setStampID( int info ) {stampID = info;} |
55 |
|
|
|
56 |
|
|
void calcForces( void ); |
57 |
|
|
double getPotential( void ); |
58 |
mmeineke |
377 |
|
59 |
mmeineke |
435 |
void printMe( void ); |
60 |
mmeineke |
377 |
|
61 |
mmeineke |
449 |
void getCOM( double COM[3] ); |
62 |
mmeineke |
452 |
void moveCOM( double delta[3] ); |
63 |
gezelter |
475 |
double getCOMvel( double COMvel[3] ); |
64 |
gezelter |
446 |
|
65 |
mmeineke |
377 |
private: |
66 |
|
|
|
67 |
mmeineke |
407 |
int stampID; // the ID in the BASS component stamp array |
68 |
|
|
int nAtoms; // the number of atoms in the molecule |
69 |
|
|
int nBonds; // ... .. .. . .bonds .. .. . . . . |
70 |
|
|
int nBends; // . . . . .. . .bends . . . . .. . |
71 |
gezelter |
416 |
int nTorsions; // .. . . .. . . torsions . . .. . . |
72 |
|
|
int nOriented; // .. . . . .. . oriented atoms . . . |
73 |
mmeineke |
423 |
int nMembers; // .. . . . . . .atoms (legacy code) . . . |
74 |
|
|
int nExcludes; // . . . . .. .. excludes .. . . |
75 |
mmeineke |
407 |
|
76 |
chuckv |
438 |
int myIndex; // mostly just for debug |
77 |
|
|
|
78 |
mmeineke |
407 |
Atom** myAtoms; // the array of atoms |
79 |
mmeineke |
423 |
Bond** myBonds; // arrays of all the short range interactions |
80 |
mmeineke |
407 |
Bend** myBends; |
81 |
|
|
Torsion** myTorsions; |
82 |
mmeineke |
423 |
Exclude** myExcludes; // array of the excluded pairs of long range forces |
83 |
mmeineke |
407 |
|
84 |
mmeineke |
377 |
}; |
85 |
|
|
|
86 |
|
|
#endif |