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. |
91 |
|
int useReactionField; |
92 |
|
int useGB; |
93 |
|
int useEAM; |
94 |
+ |
int useMolecularCutoffs; |
95 |
|
|
96 |
|
bool useInitXSstate; |
97 |
|
double orthoTolerance; |
124 |
|
char statusName[300]; // the name of the stat file to be written |
125 |
|
|
126 |
|
int seed; //seed for random number generator |
127 |
+ |
|
128 |
+ |
|
129 |
+ |
vector<double> mfact; |
130 |
+ |
int ngroup; |
131 |
+ |
vector<int> groupList; |
132 |
+ |
vector<int> groupStart; |
133 |
+ |
|
134 |
|
// refreshes the sim if things get changed (load balanceing, volume |
135 |
|
// adjustment, etc.) |
136 |
|
|
150 |
|
int getNDF(); |
151 |
|
int getNDFraw(); |
152 |
|
int getNDFtranslational(); |
153 |
< |
|
153 |
> |
int getTotIntegrableObjects(); |
154 |
|
void setBox( double newBox[3] ); |
155 |
|
void setBoxM( double newBox[3][3] ); |
156 |
|
void getBoxM( double theBox[3][3] ); |
174 |
|
|
175 |
|
void wrapVector( double thePos[3] ); |
176 |
|
|
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 |
– |
|
177 |
|
SimState* getConfiguration( void ) { return myConfiguration; } |
178 |
|
|
179 |
|
void addProperty(GenericData* prop); |
180 |
|
GenericData* getProperty(const string& propName); |
181 |
< |
vector<GenericData*> getProperties(); |
181 |
> |
//vector<GenericData*>& getProperties() {return properties;} |
182 |
|
|
183 |
|
int getSeed(void) { return seed; } |
184 |
|
void setSeed(int theSeed) { seed = theSeed;} |
197 |
|
double distXY; |
198 |
|
double distYZ; |
199 |
|
double distZX; |
201 |
– |
|
202 |
– |
|
200 |
|
|
201 |
|
void calcHmatInv( void ); |
202 |
|
void calcBoxL(); |
203 |
|
double calcMaxCutOff(); |
204 |
|
|
208 |
– |
|
205 |
|
// private function to initialize the fortran side of the simulation |
206 |
|
setFortranSim_TD setFsimulation; |
207 |
|
|
214 |
|
|
215 |
|
}; |
216 |
|
|
217 |
+ |
void getFortranGroupArray(SimInfo* info, vector<double>& mfact, int& ngroup, |
218 |
+ |
vector<int>& groupList, vector<int>& groupStart); |
219 |
|
|
222 |
– |
|
220 |
|
#endif |