implement multiprocessing for the smooth version of the hmf as well
add convergence stats
add simplistic nipals code