35 |
|
#define BRAINS_SIMCREATOR_HPP |
36 |
|
|
37 |
|
|
38 |
+ |
#include "primitives/Molecule.hpp" |
39 |
+ |
#include "brans/SimInfo.hpp" |
40 |
+ |
#include "types/MakeStamps.hpp" |
41 |
+ |
#include "io/Globals.hpp" |
42 |
+ |
#include "UseTheForce/ForceField.hpp" |
43 |
+ |
|
44 |
|
namespace oopse { |
45 |
|
|
46 |
|
/** |
47 |
|
* @class SimCreator SimCreator.hpp "brains/SimCreator.hpp" |
48 |
< |
* @brief SimCreator is an abstract factory which is responsible for setting up the simulation. |
48 |
> |
* @brief |
49 |
|
*/ |
50 |
|
class SimCreator { |
51 |
|
public: |
52 |
< |
void parseFile(const string& mdfile); |
53 |
< |
void createSim(); |
52 |
> |
|
53 |
> |
/** |
54 |
> |
* Setup Simulation |
55 |
> |
* @return a pointer to SimInfo |
56 |
> |
* @param mdfile the meta-data file name |
57 |
> |
*/ |
58 |
> |
SimInfo* createSim(char* mdfile); |
59 |
|
|
60 |
|
protected: |
61 |
|
|
62 |
< |
virtual ForceFiled* createForceFiled(); |
63 |
< |
virtual Molecule* createMolecule(); |
64 |
< |
virtual Atom* createAtom(); |
65 |
< |
virtual RigidBody* createRigidBody(); |
66 |
< |
virtual Bond* createBond(); |
67 |
< |
virtual Bend* createBend(); |
68 |
< |
virtual Torsion* createTorsion(); |
62 |
> |
/** |
63 |
> |
* Parses the meta-data file |
64 |
> |
* @param mdfile |
65 |
> |
* @param stamps |
66 |
> |
* @param globals |
67 |
> |
*/ |
68 |
> |
void parseFile(const std::string mdFileName, MakeStamps* stamps, Globals* globals); |
69 |
> |
|
70 |
> |
#ifdef IS_MPI |
71 |
> |
/** |
72 |
> |
* Divide the molecules among the processors |
73 |
> |
*/ |
74 |
> |
void divideMolecules(); |
75 |
> |
#endif //is_mpi |
76 |
> |
|
77 |
> |
/** create the molecules belong to current processor*/ |
78 |
> |
virtual void createMolecules(); |
79 |
> |
|
80 |
> |
/** |
81 |
> |
* Sets the global index for atoms, rigidbodies and cutoff groups and fill up |
82 |
> |
* globalGroupMembership and globalMolMembership arrays which map atoms' |
83 |
> |
* global index to the global index of the groups (or molecules) they belong to. |
84 |
> |
* These array are never changed during the simulation. |
85 |
> |
*/ |
86 |
> |
void setGlobalIndex(); |
87 |
> |
|
88 |
> |
private: |
89 |
|
|
90 |
+ |
/** Extracts the molecules stamps and adds them into SimInfo class */ |
91 |
+ |
void compList(MakeStamps* stamps, Globals* globals, |
92 |
+ |
std::vector<std::pair<MoleculeStamp*, int> >& moleculeStamps) ; |
93 |
+ |
|
94 |
+ |
/** Initialize fortran, mainly set the cutoff radius and setup*/ |
95 |
+ |
void initFortran(SimInfo* info); |
96 |
+ |
|
97 |
+ |
/** Figure out the cutoff radius and pass it to fortran */ |
98 |
+ |
void setCutoff(SimInfo* info); |
99 |
|
|
100 |
< |
|
61 |
< |
|
100 |
> |
std::string mdFileName_; //save the meta-data file name which may be used later |
101 |
|
}; |
102 |
|
|
103 |
|
} //end namespace oopse |