Added a smart cylinder calculation

This is a 10000 particle beast which took 20 minutes to compute one final frame after settling.

Added numpy because calculations are reduced from hours to minutes