Method ProcessPostUpdate

avida-core/source/main/cMultiProcessWorld.cc:274–356

! Process post-update events. This method is called after each update of the local population completes. Here we check to see if this world has received any migrants from other worlds, and if so, we inject them into the local population. Note that this is an unconditional injection -- that is, migrants are "pushed" to this world. Migrants are injected according to BIRTH_METHOD. \todo

 second mode that relaxes consistency in favor of speed.
 */
void cMultiProcessWorld::ProcessPostUpdate(cAvidaContext& ctx) {
	namespace mpi = boost::mpi;
	using namespace boost;
	
	// restart the timer for this method, and get the elapsed time for the past update:
	m_pf[UPDATE] = m_update_timer.elapsed();
	m_post_update_timer.restart();
	
	// wait until we're sure that this process has sent all its messages:
	mpi::wait_all(m_reqs.begin(), m_reqs.end());
	m_reqs.clear();
	
	// at this point, we know that *this* process has sent everything.  but, we don't
	// know if it's *received* everything.  so, we're going to put in a synchronization
	// barrier, which means that every process must reach the barrier before any are allowed
	// to proceed.  since we just finished waiting for all communication to complete,
	// this means that all messages must have been received, too.
	m_mpi_world.barrier();	
	
	// now, receive all the messages, but store them in order by source and tag:
	typedef std::map<int,migration_message> rx_tag_t;
	typedef std::map<int,rx_tag_t> rx_src_t;
	rx_src_t recvd;
	optional<mpi::status> s = m_mpi_world.iprobe(mpi::any_source,mpi::any_tag);
	while(s.is_initialized()) {
		migration_message msg;
		m_mpi_world.recv(s->source(), s->tag(), msg);
		recvd[s->source()][s->tag()] = msg;		
		// any others?
		s = m_mpi_world.iprobe(mpi::any_source,mpi::any_tag);		
	}
	
	// iterate over received messages in-order, injecting genomes into our population:
	for(rx_src_t::iterator i=recvd.begin(); i!=recvd.end(); ++i) {
		for(rx_tag_t::iterator j=i->second.begin(); j!=i->second.end(); ++j) {
			//std::cerr << "**********  INJECTING *********** " << i->first << " " << j->first << std::endl;
			// ok, add this migrant to the current population
			migration_message& migrant = j->second;
			int target_cell=-1;
			
			switch(GetConfig().BIRTH_METHOD.Get()) {
				case POSITION_OFFSPRING_RANDOM: { // spatial
					// invert the orginating cell
					migrant._x = GetConfig().WORLD_X.Get() - migrant._x - 1;
					migrant._y = GetConfig().WORLD_Y.Get() - migrant._y - 1;
					target_cell = GetConfig().WORLD_Y.Get() * migrant._y + migrant._x;
					break;
				}
				case POSITION_OFFSPRING_FULL_SOUP_RANDOM: { // mass action
					target_cell = GetRandom().GetInt(GetPopulation().GetSize());
					break;
				}
				default: {
					GetDriver().RaiseFatalException(-1, "Avida-MP only supports BIRTH_METHODS 0 (POSITION_OFFSPRING_RANDOM) and 4 (POSITION_OFFSPRING_FULL_SOUP_RANDOM).");
				}
			}
			
			GetPopulation().InjectGenome(target_cell,