kern_nx, kern_ny = np.round(scotts_factor * 2 * np.pi * std_devs)
inv_cov = np.linalg.inv(cov * scotts_factor ** 2)
Implementing a 2D, FFT-based Kernel Density Estimator in python, and comparing it to the SciPy implimentation