Fixes and improvements to h_max_active
Optimisations:
- In all the ghost operations, start the counter for h_max_active at 0.
- In the star drift, only collect h_max_active of the stars that are active for feedback (i.e. the stars that will actually loop). This matches the ghost updates.
- In the star recv, only collect h_max_active of the stars that are active for feedback.
Bug fixes:
- In the star ghost, only check h > h_max_active for stars active for feedback. (bug in the debug check)
- In all the ghosts, also update the cell's h_max_active counters for particles at the global h_max or h_min limits.
Edited by Matthieu Schaller