Skip to content

Commit

Permalink
help preseve T in bandpower_rescaled Cov
Browse files Browse the repository at this point in the history
  • Loading branch information
marius311 committed Jan 5, 2023
1 parent 2787dc1 commit 9d95c6d
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/proj_lambert.jl
Original file line number Diff line number Diff line change
Expand Up @@ -379,7 +379,7 @@ function Cℓ_to_Cov(::Val{:I}, proj::ProjLambert{T,V}, (Cℓ, ℓedges, θname)
ℓbin_indices = findbin.(Ref(adapt(proj.storage, ℓedges)), proj.ℓmag)
Cov(θ) = Diagonal(LambertFourier(bandpower_rescale(C₀.diag.arr, ℓbin_indices, θ), proj))
ParamDependentOp(@eval Main let Cov=$Cov
(;$θname=$(ones(T,length(ℓedges)-1)), _...) -> Cov($θname)
(;$θname=$(ones(T,length(ℓedges)-1)), _...) -> Cov($T.($θname))
end)
end

Expand All @@ -388,7 +388,7 @@ function Cℓ_to_Cov(::Val{:P}, proj::ProjLambert{T}, (CℓEE, ℓedges, θname)
ℓbin_indices = findbin.(Ref(adapt(proj.storage, ℓedges)), proj.ℓmag)
Cov(θ) = Diagonal(LambertEBFourier(bandpower_rescale(C₀.diag.El, ℓbin_indices, θ), one(eltype(θ)) .* C₀.diag.Bl, proj))
ParamDependentOp(@eval Main let Cov=$Cov
(;$θname=$(ones(T,length(ℓedges)-1)), _...) -> Cov($θname)
(;$θname=$(ones(T,length(ℓedges)-1)), _...) -> Cov($T.($θname))
end)
end

Expand Down

0 comments on commit 9d95c6d

Please sign in to comment.