def check_pixangvecpix(vlen, ntry, nside, isnest):
base = ph.Healpix_Base(nside, "NEST" if isnest else "RING")
for cnt in range(ntry):
inp = rng.integers(low=0, high=12*nside*nside-1, size=vlen)
out = base.vec2pix(ph.ang2vec(base.pix2ang(inp)))
if not np.array_equal(inp, out):
raise ValueError("Test failed")
def check_vecpixvec(vlen, ntry, nside, isnest):
base = ph.Healpix_Base(nside, "NEST" if isnest else "RING")
for cnt in range(ntry):
inp = ph.ang2vec(random_ptg(vlen))
out = base.pix2vec(base.vec2pix(inp))
if np.any(ph.v_angle(inp, out) > base.max_pixrad()):
raise ValueError("Test failed")
def test_vecangvec(vlen):
rng = np.random.default_rng(42)
inp = random_ptg(rng, vlen)
out = ph.vec2ang(ph.ang2vec(inp))
assert_equal(np.all(np.abs(out - inp) < 1e-14), True)
def test_vecpixvec_ring(vlen, nside_ring):
base = ph.Healpix_Base(nside_ring, "RING")
rng = np.random.default_rng(42)
inp = ph.ang2vec(random_ptg(rng, vlen))
out = base.pix2vec(base.vec2pix(inp))
assert_equal(np.all(ph.v_angle(inp, out) < base.max_pixrad()), True)
def check_vecangvec(vlen, ntry):
for cnt in range(ntry):
inp = random_ptg(vlen)
out = ph.vec2ang(ph.ang2vec(inp))
if np.any(np.greater(np.abs(inp-out), 1e-10)):
raise ValueError("Test failed")
def perf_vec2pix(vlen, ntry, nside, isnest):
inp = ph.ang2vec(random_ptg(vlen))
base = ph.Healpix_Base(nside, "NEST" if isnest else "RING")
genperf(base.vec2pix, "vec2pix", inp, vlen, ntry, nside, isnest)