Fast Krylov Methods for N-Body Learning