marius311 / CMBLensing.jl

The automatically differentiable and GPU-compatible toolkit for CMB analysis.

Home Page:https://cosmicmar.com/CMBLensing.jl

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

support non-zero EB/TB covariance

kimmywu opened this issue · comments

CMBLensing currently doesn't support non-zero EB/TB covariance in its signal/noise covariance construction. See

# A Flat TEB covariance of the form:
#
# [ΣTT ΣTE ⋅
# ΣTE ΣEE ⋅
# ⋅ ⋅ ΣBB]
#
# We store the 2x2 block as a 2x2 SMatrix, ΣTE, so that we can easily call sqrt/inv on
# it, and the ΣBB block separately as ΣB.
struct FlatIEBCov{T,F} <: ImplicitOp{IEBFourier,S02,Pix}
ΣTE :: SMatrix{2,2,Diagonal{T,F},4}
ΣB :: Diagonal{T,F}
end

Would be good to add that feature to allow modeling of correlated EB/TB noise/signal. I haven't thought through all the details of how to implement this yet. So this is a reminder.