vectorised gstdint
hn integration only calls erf once if it can
re-introduced the spline approximation for hazard rate models