Pairwise KL divergence of a batch of Multivariate Gaussian Distributions

Hi,
Can someone please suggest an efficient way to calculate pairwise KL divergence of a batch of multivariate gaussian distributions?

Thanks,
Elita