lensing sources in camb

[bin hu]

Dear all,

I find that the anisotropic stress contribution from radiation and neutrino are not included in the lensing source in the subroutine "output" of "equations.f90" even in the latest version

******************************************************************
!phi_lens = Phi - 1/2 kappa (a/k)^2 sum_i rho_i pi_i
!Neglect pi contributions because not accurate at late time anyway
phi = -(dgrho +3*dgq*adotoa/k)/(k2*EV%Kf(1)*2)
! - (grhor_t*pir + grhog_t*pig+ pinu*gpnu_t)/k2

sources(3) = -2*phi*f_K(tau-tau_maxvis)/(f_K(CP%tau0-tau_maxvis)*f_K(CP%tau0-tau))
! sources(3) = -2*phi*(tau-tau_maxvis)/((CP%tau0-tau_maxvis)*(CP%tau0-tau))
!We include the lensing factor of two here
******************************************************************
As far as I know, there is the only point where the lensing sources are computed. Does the massive neutrino contribution to the lensing sources are included at some other places?

Thanks in advance!

[Antony Lewis]

No, that's it. I agree it would probably be a good idea to add those terms back in, esp. for non-vanilla models.