1 |
#ifndef __SIMSETUP_H__ |
2 |
#define __SIMSETUP_H__ |
3 |
|
4 |
#include "MakeStamps.hpp" |
5 |
#include "Globals.hpp" |
6 |
#include "ForceFields.hpp" |
7 |
#include "SimInfo.hpp" |
8 |
#include "ReadWrite.hpp" |
9 |
|
10 |
// this routine is defined in BASS_interface.cpp |
11 |
void set_interface_stamps( MakeStamps* ms, Globals* g ); |
12 |
|
13 |
class SimSetup{ |
14 |
|
15 |
public: |
16 |
SimSetup(); |
17 |
~SimSetup(); |
18 |
|
19 |
void setSimInfo( SimInfo* the_simnfo ) { simnfo = the_simnfo; } |
20 |
void parseFile( char* fileName ); |
21 |
void createSim( void ); |
22 |
|
23 |
|
24 |
private: |
25 |
|
26 |
#ifdef IS_MPI |
27 |
void receiveParse(void); |
28 |
#endif |
29 |
|
30 |
void makeAtoms( void ); |
31 |
void makeBonds( void ); |
32 |
void makeBends( void ); |
33 |
void makeTorsions( void ); |
34 |
|
35 |
void initFromBass( void ); |
36 |
void makeElement( double x, double y, double z ); |
37 |
|
38 |
MakeStamps* stamps; |
39 |
Globals* globals; |
40 |
char* inFileName; |
41 |
|
42 |
SimInfo* simnfo; |
43 |
|
44 |
int n_components; |
45 |
|
46 |
char force_field[100]; |
47 |
char ensemble[100]; |
48 |
Component** the_components; |
49 |
|
50 |
int* components_nmol; |
51 |
MoleculeStamp** comp_stamps; //the stamps matching the components |
52 |
int tot_nmol; |
53 |
int tot_atoms; |
54 |
int tot_bonds; |
55 |
int tot_bends; |
56 |
int tot_torsions; |
57 |
int tot_SRI; |
58 |
|
59 |
Atom** the_atoms; |
60 |
SRI** the_sris; |
61 |
ex_pair *the_excludes; |
62 |
Molecule* the_molecules; |
63 |
ForceFields* the_ff; |
64 |
|
65 |
// needed by makeElement |
66 |
|
67 |
int current_mol; |
68 |
int current_comp_mol; |
69 |
int current_comp; |
70 |
int current_atom_ndx; |
71 |
|
72 |
#ifdef IS_MPI |
73 |
int* globalIndex; |
74 |
#endif //is_mpi |
75 |
|
76 |
}; |
77 |
#endif |