Mixture of small improvements:
- Do not call external functions for min and max.
- Do not use double-precision constants.
- Reduce the number of divisions.
I believe this speeds up the Evrad collapse by 15-20%. @jborrow could you confirm this on your accuracy vs. time plot?
@bvandenbroucke I am likely to have made a mistake somewhere. Could you cross-check that it is all fine?