def test_renorm_demote_symmetric():
"""Test that Region._renorm and Region._demote are mutual inverses"""
ra = 13.5
dec = -90
radius = 0.1
region = Region(maxdepth=11)
region.add_circles(np.radians(ra), np.radians(dec), np.radians(radius))
region._demote_all()
start_dict = region.pixeldict.copy()
region._renorm()
region._demote_all()
end_dict = region.pixeldict.copy()
test = True
for i in range(1, region.maxdepth+1):
if len(end_dict[i].difference(start_dict[i])) > 0:
test = False
if not (test): raise AssertionError('renorm and demote are not symmetric')
def test_renorm_demote_symmetric():
"""Test that Region._renorm and Region._demote are mutual inverses"""
ra = 13.5
dec = -90
radius = 0.1
region = Region(maxdepth=11)
region.add_circles(np.radians(ra), np.radians(dec), np.radians(radius))
region._demote_all()
start_dict = region.pixeldict.copy()
region._renorm()
region._demote_all()
end_dict = region.pixeldict.copy()
test = True
for i in range(1, region.maxdepth + 1):
if len(end_dict[i].difference(start_dict[i])) > 0:
test = False
if not (test): raise AssertionError('renorm and demote are not symmetric')
def test_add_circles_list_scalar():
"""Test that Region.add_circles works for vector inputs"""
ra_list = np.radians([13.5, 13.5])
dec_list = np.radians([-90, -90])
radius_list = np.radians([0.1, 0.01])
ra = ra_list[0]
dec = dec_list[0]
radius = radius_list[0]
region1 = Region(maxdepth=11)
region2 = Region(maxdepth=11)
region1.add_circles(ra_list, dec_list, radius_list)
region1._demote_all()
region2.add_circles(ra, dec, radius)
region2._demote_all()
test = True
for i in range(1, region1.maxdepth+1):
if len(region1.pixeldict[i].difference(region2.pixeldict[i])) > 0:
test = False
if not (test): raise AssertionError('add_circles gives different results for lists and scalars')
def test_add_circles_list_scalar():
"""Test that Region.add_circles works for vector inputs"""
ra_list = np.radians([13.5, 13.5])
dec_list = np.radians([-90, -90])
radius_list = np.radians([0.1, 0.01])
ra = ra_list[0]
dec = dec_list[0]
radius = radius_list[0]
region1 = Region(maxdepth=11)
region2 = Region(maxdepth=11)
region1.add_circles(ra_list, dec_list, radius_list)
region1._demote_all()
region2.add_circles(ra, dec, radius)
region2._demote_all()
test = True
for i in range(1, region1.maxdepth + 1):
if len(region1.pixeldict[i].difference(region2.pixeldict[i])) > 0:
test = False
if not (test):
raise AssertionError(
'add_circles gives different results for lists and scalars')
