6 |
|
#include "parse_me.h" |
7 |
|
#include "LRI.hpp" |
8 |
|
#include "Integrator.hpp" |
9 |
– |
#include "mpiInterface.h" |
9 |
|
|
10 |
+ |
#ifdef IS_MPI |
11 |
+ |
#include "mpiBASS.h" |
12 |
+ |
#include "bassDiag.hpp" |
13 |
+ |
#endif |
14 |
+ |
|
15 |
|
SimSetup::SimSetup(){ |
16 |
|
stamps = new MakeStamps(); |
17 |
|
globals = new Globals(); |
26 |
|
|
27 |
|
inFileName = fileName; |
28 |
|
set_interface_stamps( stamps, globals ); |
29 |
< |
#ifdef MPI |
29 |
> |
#ifdef IS_MPI |
30 |
|
mpiEventInit(); |
31 |
|
#endif |
32 |
|
yacc_BASS( fileName ); |
33 |
< |
#ifdef MPI |
33 |
> |
#ifdef IS_MPI |
34 |
|
throwMPIEvent(NULL); |
35 |
|
#endif |
36 |
|
|
37 |
|
} |
38 |
|
|
39 |
< |
#ifdef MPI |
39 |
> |
#ifdef IS_MPI |
40 |
|
void SimSetup::receiveParse(void){ |
41 |
|
|
42 |
|
set_interface_stamps( stamps, globals ); |
44 |
|
mpiEventLoop(); |
45 |
|
|
46 |
|
} |
43 |
– |
#endif |
47 |
|
|
48 |
+ |
|
49 |
|
void SimSetup::testMe(void){ |
50 |
|
bassDiag* dumpMe = new bassDiag(globals,stamps); |
51 |
|
dumpMe->dumpStamps(); |
52 |
|
delete dumpMe; |
53 |
|
} |
54 |
< |
|
54 |
> |
#endif |
55 |
|
void SimSetup::createSim( void ){ |
56 |
|
|
57 |
|
MakeStamps *the_stamps; |
156 |
|
// create the atom and short range interaction arrays |
157 |
|
|
158 |
|
the_atoms = new Atom*[tot_atoms]; |
159 |
+ |
Atom::createArrays(tot_atoms); |
160 |
|
the_molecules = new Molecule[tot_nmol]; |
161 |
|
|
162 |
|
|
393 |
|
current_atom = comp_stamps[i]->getAtom( k ); |
394 |
|
if( current_atom->haveOrientation() ){ |
395 |
|
|
396 |
< |
dAtom = new DirectionalAtom; |
396 |
> |
dAtom = new DirectionalAtom(index); |
397 |
|
simnfo->n_oriented++; |
398 |
|
the_atoms[index] = dAtom; |
399 |
|
|
413 |
|
dAtom->setSUz( uz ); |
414 |
|
} |
415 |
|
else{ |
416 |
< |
the_atoms[index] = new GeneralAtom; |
416 |
> |
the_atoms[index] = new GeneralAtom(index); |
417 |
|
} |
418 |
|
the_atoms[index]->setType( current_atom->getType() ); |
419 |
|
the_atoms[index]->setIndex( index ); |