Skip to content

Instantly share code, notes, and snippets.

@zonca
Created July 13, 2011 18:16
Show Gist options
  • Save zonca/1080918 to your computer and use it in GitHub Desktop.
Save zonca/1080918 to your computer and use it in GitHub Desktop.
healpy anafast issue with unseen points
m = np.ones(healpy.nside2npix(64))
cl = healpy.anafast(m) # array([ 12.56637061, 0. ,0. .... all zeros
m[100:200]=healpy.UNSEEN
clunseen = healpy.anafast(m) # array([ 1.39439151e+56, 1.37741024e+56,
m[100:200]=0
clzero = healpy.anafast(m) # array([ 1.25152959e+01, 5.13689503e-05, 5.01205239e-05,
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment