Example #1
0
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)