def test_log_mean_exp_fail(): X = np.abs(np.random.randn(100)) X[0] = -3000 direct_exp_log_mean_exp = np.exp(log_mean_exp(X)) assert_equal(direct_exp_log_mean_exp, np.inf)
def test_log_mean_exp_equals_avg_prob_of_log_probs(): X = np.abs(np.random.randn(100)) direct_exp_log_mean_exp = np.exp(log_mean_exp(X)) safe = avg_prob_of_log_probs(X) assert_allclose(direct_exp_log_mean_exp, safe)
def test_log_mean_exp(): X = np.abs(np.random.randn(100)) direct = np.log(np.mean(np.exp(X))) indirect = log_mean_exp(X) assert_allclose(direct, indirect)