298 |
|
simnfo->n_torsions = tot_torsions; |
299 |
|
simnfo->n_SRI = tot_SRI; |
300 |
|
simnfo->n_mol = tot_nmol; |
301 |
+ |
|
302 |
|
simnfo->molMembershipArray = new int[tot_atoms]; |
303 |
|
|
304 |
|
#ifdef IS_MPI |
348 |
|
} |
349 |
|
local_SRI = local_bonds + local_bends + local_torsions; |
350 |
|
|
350 |
– |
if (worldRank != 0) { |
351 |
– |
for (i =0; i < tot_atoms; i++){ |
352 |
– |
std::cerr << "i = " << i << " molMembershipArray[i] = " << simnfo->molMembershipArray[i] << "\n"; |
353 |
– |
} |
354 |
– |
} |
355 |
– |
|
351 |
|
simnfo->n_atoms = mpiSim->getMyNlocal(); |
352 |
|
|
353 |
|
if( local_atoms != simnfo->n_atoms ){ |
745 |
|
else if( !strcmp( force_field, "LJ" ) ){ |
746 |
|
new Verlet( *simnfo, the_ff, the_extendedsystem ); |
747 |
|
} |
748 |
< |
else { |
754 |
< |
std::cerr << "I'm a bug.\n"; |
755 |
< |
fprintf( stderr, "Ima bug. stderr %s\n", force_field); |
756 |
< |
} |
748 |
> |
|
749 |
|
#ifdef IS_MPI |
750 |
|
mpiSim->mpiRefresh(); |
751 |
|
#endif |