From 05731412a45e8df97af866ac86913d5ea409c5a0 Mon Sep 17 00:00:00 2001 From: ismael2395 Date: Wed, 11 Dec 2024 10:38:55 -0800 Subject: [PATCH] likelihood update comments --- bpd/likelihood.py | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/bpd/likelihood.py b/bpd/likelihood.py index f12ca10..9197155 100644 --- a/bpd/likelihood.py +++ b/bpd/likelihood.py @@ -30,12 +30,11 @@ def shear_loglikelihood( _, _, _ = e_post.shape # (N, K, 2) _prior = partial(prior, sigma=sigma_e) + # denom e_post_mag = norm(e_post, axis=-1) denom = interim_prior(e_post_mag) # (N, K), can ignore angle in prior as uniform - # for num, use trick - # p(w_n' | g, alpha ) = p(w_n' \cross^{-1} g | alpha ) = p(w_n | alpha) * |jac(w_n / w_n')| - + # for num, need to include Jacobian of shear transoformation # shape = (N, K, 2) grad1 = _grad_fnc1(e_post, g) grad2 = _grad_fnc2(e_post, g)