“An unbiased, low-variance estimator for `KL[q,p]` is given by the expression `(r - 1) - log(r)`, where `r = p(x)/q(x)`.”