6 |
|
#include <vector> |
7 |
|
|
8 |
|
#include "Atom.hpp" |
9 |
+ |
#include "RigidBody.hpp" |
10 |
|
#include "Molecule.hpp" |
11 |
+ |
#include "Exclude.hpp" |
12 |
+ |
#include "SkipList.hpp" |
13 |
|
#include "AbstractClasses.hpp" |
14 |
|
#include "MakeStamps.hpp" |
15 |
|
#include "SimState.hpp" |
20 |
|
#include "GenericData.hpp" |
21 |
|
//#include "Minimizer.hpp" |
22 |
|
//#include "OOPSEMinimizer.hpp" |
23 |
< |
|
23 |
> |
double roundMe( double x ); |
24 |
|
class OOPSEMinimizer; |
25 |
|
class SimInfo{ |
26 |
|
|
31 |
|
|
32 |
|
int n_atoms; // the number of atoms |
33 |
|
Atom **atoms; // the array of atom objects |
34 |
+ |
|
35 |
+ |
vector<RigidBody*> rigidBodies; // A vector of rigid bodies |
36 |
+ |
vector<StuntDouble*> integrableObjects; |
37 |
|
|
38 |
|
double tau[9]; // the stress tensor |
39 |
|
|
51 |
|
|
52 |
|
int n_dipoles; // number of dipoles |
53 |
|
|
54 |
< |
|
55 |
< |
int n_exclude; // the # of pairs excluded from long range forces |
50 |
< |
Exclude** excludes; // the pairs themselves |
51 |
< |
|
54 |
> |
int n_exclude; |
55 |
> |
Exclude* excludes; // the exclude list for ignoring pairs in fortran |
56 |
|
int nGlobalExcludes; |
57 |
|
int* globalExcludes; // same as above, but these guys participate in |
58 |
|
// no long range forces. |
142 |
|
int getNDF(); |
143 |
|
int getNDFraw(); |
144 |
|
int getNDFtranslational(); |
145 |
< |
|
145 |
> |
int getTotIntegrableObjects(); |
146 |
|
void setBox( double newBox[3] ); |
147 |
|
void setBoxM( double newBox[3][3] ); |
148 |
|
void getBoxM( double theBox[3][3] ); |
166 |
|
|
167 |
|
void wrapVector( double thePos[3] ); |
168 |
|
|
165 |
– |
void matMul3(double a[3][3], double b[3][3], double out[3][3]); |
166 |
– |
void matVecMul3(double m[3][3], double inVec[3], double outVec[3]); |
167 |
– |
void invertMat3(double in[3][3], double out[3][3]); |
168 |
– |
void transposeMat3(double in[3][3], double out[3][3]); |
169 |
– |
void printMat3(double A[3][3]); |
170 |
– |
void printMat9(double A[9]); |
171 |
– |
double matDet3(double m[3][3]); |
172 |
– |
double matTrace3(double m[3][3]); |
173 |
– |
|
174 |
– |
void crossProduct3(double a[3],double b[3], double out[3]); |
175 |
– |
double dotProduct3(double a[3], double b[3]); |
176 |
– |
double length3(double a[3]); |
177 |
– |
|
169 |
|
SimState* getConfiguration( void ) { return myConfiguration; } |
170 |
|
|
171 |
|
void addProperty(GenericData* prop); |
172 |
|
GenericData* getProperty(const string& propName); |
173 |
< |
vector<GenericData*> getProperties(); |
173 |
> |
//vector<GenericData*>& getProperties() {return properties;} |
174 |
|
|
175 |
|
int getSeed(void) { return seed; } |
176 |
|
void setSeed(int theSeed) { seed = theSeed;} |
189 |
|
double distXY; |
190 |
|
double distYZ; |
191 |
|
double distZX; |
201 |
– |
|
202 |
– |
|
192 |
|
|
193 |
|
void calcHmatInv( void ); |
194 |
|
void calcBoxL(); |
195 |
|
double calcMaxCutOff(); |
196 |
|
|
208 |
– |
|
197 |
|
// private function to initialize the fortran side of the simulation |
198 |
|
setFortranSim_TD setFsimulation; |
199 |
|
|
206 |
|
|
207 |
|
}; |
208 |
|
|
221 |
– |
|
222 |
– |
|
209 |
|
#endif |