Numerics


Best option is a re-write of the integration and Dirac ODE methods.

particularly, additivePIntegral(), required for homogeneous, is very sensitive to method, and not 100% sure method is optimal!

Affects TDHF the most, but also just HF!