Anonymous

void cell::update( ){	particles_t::iterator end = m_particles.end();	for (particles_t::iterator part_it = m_particles.begin(); part_it != end; ++part_it) {		particle& p = **part_it;		p.r = p.r1;		p.v = p.v1;		p.r1 = 0.0;		p.v1 = 0.0;		if (m_volume->outside_cell(p.r,m_i,m_j,m_k)) {			m_volume->lost(p);			m_particles.erase(part_it);		}	}}