newmatmcx=qml.allmat2bigmat(bigmatmcx) clf() imshow(np.log10(np.abs(newmatmcx)),interpolation='nearest') colorbar() #### calcul des matrices de cov theoriques I,Q,U reload(qml) ip=np.arange(12*nside**2) ipok=ip[~mask] lmax=3*nside maskl=ell<(lmax+1) bl=np.ones(maskl.size) #spectra=[ell,ctt,cee,cbb,cte] #### This is the new ordering allcovth=qml.covth(nside,ipok,lmax,bl,spectra,polar=True,temp=True,allinone=False) newmatth=qml.allmat2bigmat(allcovth) mini = 0 maxi = 1e4 clf() subplot(2,3,1) imshow(covii,interpolation='nearest', vmin=mini, vmax=maxi) colorbar() title('Auto') subplot(2,3,2) imshow(coviix,interpolation='nearest', vmin=mini, vmax=maxi) colorbar() title('Cross') subplot(2,3,3)