Pixell: A Powerful Library for Rectangular Pixel Map Analysis and Visualization
The output is in the shape one would expect from numpy broadcasting. For example,in the last example, the TE power spectrum would be found in cl[0,1], and theET power spectrum (which is different for the cross-spectrum case) is in cl[1,0].If a Healpix-style compressed spectrum is desired, use pixell.powspec.sym_compress.