Use FFT to replace inefficient FT algorithm, not perfectly working now, rdm and kappa right, but emb basis probably wrong