To evaluate the density of data at x, the current interface is density = np.exp(kde.score_samples(data))

(see in action here)

This seems really silly and unnecessary for such a common operation.

How would people feel about adding a density method to KernelDensity and to GMM to make this easier for users to figure out?