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