def test_ring_phi0(): counts = get_ring_pixel_counts(Nside) p0_arr = np.r_[0, np.cumsum(counts)[:-1]] from cmb.mapdatautils import array_vec2ang phi_fid = array_vec2ang( healpix.pix2vec_ring(Nside, np.arange(12 * Nside**2)))[:, 1] phi0 = get_ring_phi0(Nside) yield assert_almost_equal, phi0, phi_fid[p0_arr]
def test_ring_thetas(): theta = get_ring_thetas(Nside) eq_(theta.shape, (4 * Nside - 1,)) hz = healpix.pix2vec_ring(Nside, np.arange(12 * Nside**2))[:, 2] hz = np.unique(hz) htheta = np.arccos(hz)[::-1] assert_almost_equal(theta, htheta) thetapos = get_ring_thetas(Nside, positive_only=True) eq_(thetapos[0], np.arccos(0)) assert_almost_equal(thetapos, theta[:2 * Nside][::-1])