Fixed MPI by changing the indexing scheme for reading from the source bank. Previously, it was doing index_source -1, with index source starting at 0, which was causing undefined behavior. Now, it uses index_source properly. This may have been an artifact from Fortran days?

This commit is contained in:
John Tramm 2019-11-20 20:07:27 +00:00
parent d99e2d27c1
commit c7399d0d1d

View file

@ -1106,8 +1106,10 @@ void finalize_generation()
void initialize_history(Particle* p, int64_t index_source)
{
//assert(index_source - 1 >= 0 );
// set defaults
p->from_source(&simulation::source_bank[index_source - 1]);
//p->from_source(&simulation::source_bank[index_source - 1]);
p->from_source(&simulation::source_bank[index_source]);
// set identifier for particle
p->id_ = simulation::work_index[mpi::rank] + index_source;