Something we must not forget: We discussed that we can, similar to #850, calculate all the constants that only depend on inputs and the reference particle in compute_constants, instead of doing them per particle in an element's spin_and_phasespace_push method.
