temporary fix for T+P bandpower grad being zero

src/CMBLensing.jl

@@ -79,7 +79,7 @@ import Base: +, -, *, \, /, ^, ~, ≈, <, <=, |, &, ==, !,
 import Base.Broadcast: materialize, preprocess, broadcasted
 import Zygote.ChainRules: rrule
 import LinearAlgebra: checksquare, diag, dot, isnan, ldiv!, logdet, mul!, norm,
-    pinv, StructuredMatrixStyle, structured_broadcast_alloc, tr
+    pinv, StructuredMatrixStyle, structured_broadcast_alloc, tr, det
 import Measurements: ±
 import Statistics: std
 import ChainRules: ProjectTo

src/field_vectors.jl

@@ -66,19 +66,19 @@ promote_rule(::Type{F}, ::Type{<:Scalar}) where {F<:Field} = F
 end
 
 @auto_adjoint function sqrt(A::SA) where {SA<:StaticMatrix{2,2,<:DiagOp}}
-    a,c,b,d = A
+    a, c, b, d = A[1,1], A[2,1], A[2,1], A[2,2]
     s = sqrt(a*d-b*c)
     t = pinv(sqrt(a+(d+2s)))
     SA([t*(a+s) t*b; t*c t*(d+s)])
 end
 
 @auto_adjoint function det(A::StaticMatrix{2,2,<:DiagOp})
-    a,c,b,d = A
+    a, c, b, d = A[1,1], A[2,1], A[2,1], A[2,2]
     a*d-b*c
 end
 
 @auto_adjoint function pinv(A::SA) where {SA<:StaticMatrix{2,2,<:DiagOp}}
-    a,c,b,d = A
+    a, c, b, d = A[1,1], A[2,1], A[2,1], A[2,2]
     idet = pinv(a*d-b*c)
     SA([d*idet -(b*idet); -(c*idet) a*idet])
 end