The current active_particles and eachparticle indexing approach is confusing.
Currently:
eachparticle contains both active and inactive particles
active_particles returns a view of only the active particles from eachparticle
We need to rethink this design:
If we don't need buffer particles, we shouldn't have to deal with the active/inactive logic at all.
Thus, eachparticle should only contain active particles by default
Add a separate iterator, e.g. each_active_and_inactive_particle(), for the rare cases where you need all particles. (is this ever happening?)
The current
active_particlesandeachparticleindexing approach is confusing.Currently:
eachparticlecontains both active and inactive particlesactive_particlesreturns a view of only the active particles fromeachparticleWe need to rethink this design:
If we don't need buffer particles, we shouldn't have to deal with the active/inactive logic at all.
Thus,
eachparticleshould only contain active particles by defaultAdd a separate iterator, e.g.
each_active_and_inactive_particle(), for the rare cases where you need all particles. (is this ever happening?)