def calc_log_normalizer(mu, sigma, l, h): return np.log( norm.cdf(h, loc=mu, scale=sigma) - norm.cdf(l, loc=mu, scale=sigma))