1
1
# pylint: disable=E1101, C0103, R0912, R0913, R0914, R0915, W0212
2
2
3
- #Copyright © 2016 United States Government as represented by the Administrator
3
+ #Copyright 2016 United States Government as represented by the Administrator
4
4
#of the National Aeronautics and Space Administration. All Rights Reserved.
5
5
6
6
'''
@@ -164,6 +164,7 @@ def apply_coupling(cls_in, Mscal, Mpol, Mcross, inverse=True):
164
164
TT , EE , BB , TE , EB , TB = 0 , 1 , 2 , 3 , 4 , 5
165
165
doV = False
166
166
elif ncls == 4 :
167
+ doV = False
167
168
TT , EE , BB , TE = 0 , 1 , 2 , 3
168
169
169
170
cls_out = np .empty ([ncls , nell ])
@@ -202,7 +203,9 @@ def apply_coupling(cls_in, Mscal, Mpol, Mcross, inverse=True):
202
203
203
204
cls_out [EE , :] = polvect_out [0 :nell ]
204
205
cls_out [BB , :] = polvect_out [nell :2 * nell ]
205
- cls_out [EB , :] = polvect_out [2 * nell :]
206
+
207
+ if ncls > 4 :
208
+ cls_out [EB , :] = polvect_out [2 * nell :]
206
209
207
210
crossvect = np .zeros (2 * nell )
208
211
crossvect [:nell ] = cls_in [TE ]
@@ -212,7 +215,9 @@ def apply_coupling(cls_in, Mscal, Mpol, Mcross, inverse=True):
212
215
crossvect_out = np .dot (B_cross , crossvect )
213
216
214
217
cls_out [TE , :] = crossvect_out [:nell ]
215
- cls_out [TB , :] = crossvect_out [nell :]
218
+
219
+ if ncls > 4 :
220
+ cls_out [TB , :] = crossvect_out [nell :]
216
221
217
222
if doV :
218
223
crossvect [:nell ] = cls_in [EV ]
0 commit comments