ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/SimInfo.hpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/SimInfo.hpp (file contents):
Revision 644 by mmeineke, Tue Jul 22 16:41:08 2003 UTC vs.
Revision 1139 by gezelter, Wed Apr 28 22:06:29 2004 UTC

# Line 1 | Line 1
1   #ifndef __SIMINFO_H__
2   #define __SIMINFO_H__
3  
4 + #include <map>
5 + #include <string>
6 + #include <vector>
7  
5 –
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"
16  
17   #define __C
18   #include "fSimulation.h"
19   #include "fortranWrapDefines.hpp"
20 <
21 <
22 <
20 > #include "GenericData.hpp"
21 > //#include "Minimizer.hpp"
22 > //#include "OOPSEMinimizer.hpp"
23 > double roundMe( double x );
24 > class OOPSEMinimizer;
25   class SimInfo{
26  
27   public:
28  
29    SimInfo();
30 <  ~SimInfo(){}
30 >  ~SimInfo();
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  
40 <  unsigned int n_bonds;    // number of bends
41 <  unsigned int n_bends;    // number of bends
42 <  unsigned int n_torsions; // number of torsions
43 <  unsigned int n_oriented; // number of of atoms with orientation
44 <  unsigned int ndf;        // number of actual degrees of freedom
45 <  unsigned int ndfRaw;     // number of settable degrees of freedom
40 >  int n_bonds;    // number of bends
41 >  int n_bends;    // number of bends
42 >  int n_torsions; // number of torsions
43 >  int n_oriented; // number of of atoms with orientation
44 >  int ndf;        // number of actual degrees of freedom
45 >  int ndfRaw;     // number of settable degrees of freedom
46 >  int ndfTrans;   // number of translational degrees of freedom
47 >  int nZconstraints; // the number of zConstraints
48  
49 <  unsigned int setTemp;   // boolean to set the temperature at each sampleTime
49 >  int setTemp;   // boolean to set the temperature at each sampleTime
50 >  int resetIntegrator; // boolean to reset the integrator
51  
52 <  unsigned int n_dipoles; // number of dipoles
52 >  int n_dipoles; // number of dipoles
53  
54 <
55 <  int n_exclude;  // the # of pairs excluded from long range forces
42 <  Exclude** excludes;       // the pairs themselves
43 <
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.
# Line 50 | Line 62 | class SimInfo{ (public)
62  
63    int n_constraints; // the number of constraints on the system
64  
65 <  unsigned int n_SRI;   // the number of short range interactions
65 >  int n_SRI;   // the number of short range interactions
66  
67    double lrPot; // the potential energy from the long range calculations.
68  
# Line 74 | Line 86 | class SimInfo{ (public)
86    int usePBC; // whether we use periodic boundry conditions.
87    int useLJ;
88    int useSticky;
89 <  int useDipole;
89 >  int useCharges;
90 >  int useDipoles;
91    int useReactionField;
92    int useGB;
93    int useEAM;
94 +  int useMolecularCutoffs;
95    
96 +  bool useInitXSstate;
97 +  double orthoTolerance;
98  
99    double dt, run_time;           // the time step and total time
100    double sampleTime, statusTime; // the position and energy dump frequencies
101    double target_temp;            // the target temperature of the system
102    double thermalTime;            // the temp kick interval
103    double currentTime;            // Used primarily for correlation Functions
104 +  double resetTime;              // Use to reset the integrator periodically
105  
106    int n_mol;           // n_molecules;
107    Molecule* molecules; // the array of molecules
108    
109 <  int nComponents;           // the number of componentsin the system
109 >  int nComponents;           // the number of components in the system
110    int* componentsNmol;       // the number of molecules of each component
111    MoleculeStamp** compStamps;// the stamps matching the components
112    LinkedMolStamp* headStamp; // list of stamps used in the simulation
# Line 99 | Line 116 | class SimInfo{ (public)
116    char mixingRule[100]; // the mixing rules for Lennard jones/van der walls
117    BaseIntegrator *the_integrator; // the integrator of the simulation
118  
119 +  OOPSEMinimizer* the_minimizer; // the energy minimizer
120 +  bool has_minimizer;
121 +
122    char finalName[300];  // the name of the eor file to be written
123    char sampleName[300]; // the name of the dump file to be written
124    char statusName[300]; // the name of the stat file to be written
125  
126 <
126 >  int seed;                    //seed for random number generator
127    // refreshes the sim if things get changed (load balanceing, volume
128    // adjustment, etc.)
129  
# Line 112 | Line 132 | class SimInfo{ (public)
132  
133    // sets the internal function pointer to fortran.
134  
135 <  void setInternal( void (*fSetup) setFortranSimList,
136 <                    void (*fBox) setFortranBoxList,
137 <                    void (*fCut) notifyFortranCutOffList ){
135 >  void setInternal( setFortranSim_TD fSetup,
136 >                    setFortranBox_TD fBox,
137 >                    notifyFortranCutOff_TD fCut){
138      setFsimulation = fSetup;
139      setFortranBoxSize = fBox;
140      notifyFortranCutOffs = fCut;
# Line 122 | Line 142 | class SimInfo{ (public)
142  
143    int getNDF();
144    int getNDFraw();
145 <
145 >  int getNDFtranslational();
146 >  int getTotIntegrableObjects();
147    void setBox( double newBox[3] );
148    void setBoxM( double newBox[3][3] );
149    void getBoxM( double theBox[3][3] );
150    void scaleBox( double scale );
151    
152 <  void setRcut( double theRcut );
153 <  void setEcr( double theEcr );
154 <  void setEcr( double theEcr, double theEst );
152 >  void setDefaultRcut( double theRcut );
153 >  void setDefaultEcr( double theEcr );
154 >  void setDefaultEcr( double theEcr, double theEst );
155 >  void checkCutOffs( void );
156  
157    double getRcut( void )  { return rCut; }
158    double getRlist( void ) { return rList; }
159    double getEcr( void )   { return ecr; }
160    double getEst( void )   { return est; }
161 +  double getMaxCutoff( void ) { return maxCutoff; }
162  
163    void setTime( double theTime ) { currentTime = theTime; }
164 <  void incrTime( double dt ) { currentTime += dt; }
165 <  void decrTime( double dt ) { currentTime -= dt; }
164 >  void incrTime( double the_dt ) { currentTime += the_dt; }
165 >  void decrTime( double the_dt ) { currentTime -= the_dt; }
166    double getTime( void ) { return currentTime; }
167  
168    void wrapVector( double thePos[3] );
169  
170 <  void matMul3(double a[3][3], double b[3][3], double out[3][3]);
148 <  void matVecMul3(double m[3][3], double inVec[3], double outVec[3]);
149 <  void invertMat3(double in[3][3], double out[3][3]);
150 <  void transposeMat3(double in[3][3], double out[3][3]);
151 <  void printMat3(double A[3][3]);
152 <  void printMat9(double A[9]);
153 <  double matDet3(double m[3][3]);
170 >  SimState* getConfiguration( void ) { return myConfiguration; }
171    
172 +  void addProperty(GenericData* prop);
173 +  GenericData* getProperty(const string& propName);
174 +  //vector<GenericData*>& getProperties()  {return properties;}    
175 +
176 +  int getSeed(void) {  return seed; }
177 +  void setSeed(int theSeed) {  seed = theSeed;}
178 +
179   private:
180  
181 <  double origRcut, origEcr;
158 <  int boxIsInit, haveOrigRcut, haveOrigEcr;
181 >  SimState* myConfiguration;
182  
183 <  double oldEcr;
161 <  double oldRcut;
183 >  int boxIsInit, haveRcut, haveEcr;
184  
185    double rList, rCut; // variables for the neighborlist
186    double ecr;             // the electrostatic cutoff radius
187    double est;             // the electrostatic skin thickness
188    double maxCutoff;
189 +
190 +  double distXY;
191 +  double distYZ;
192 +  double distZX;
193    
194    void calcHmatInv( void );
195    void calcBoxL();
196 <  void checkCutOffs( void );
196 >  double calcMaxCutOff();
197  
198    // private function to initialize the fortran side of the simulation
199 <  void (*setFsimulation) setFortranSimList;
199 >  setFortranSim_TD setFsimulation;
200  
201 <  void (*setFortranBoxSize) setFortranBoxList;
201 >  setFortranBox_TD setFortranBoxSize;
202    
203 <  void (*notifyFortranCutOffs) notifyFortranCutOffList;
203 >  notifyFortranCutOff_TD notifyFortranCutOffs;
204 >  
205 >  //Addtional Properties of SimInfo
206 >  map<string, GenericData*> properties;
207  
208   };
209  
181 –
182 –
210   #endif

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines