def test_coadd_unit():
geom = WcsGeom.create(npix=(10, 10), binsz=1, proj="CAR", frame="galactic")
m1 = WcsNDMap(geom, data=np.ones((10, 10)), unit="m2")
m2 = WcsNDMap(geom, data=np.ones((10, 10)), unit="cm2")
m1.coadd(m2)
assert_allclose(m1.data, 1.0001)
def test_wcsndmap_coadd(npix, binsz, frame, proj, skydir, axes):
geom = WcsGeom.create(
npix=npix, binsz=binsz, skydir=skydir, proj=proj, frame=frame, axes=axes
)
m0 = WcsNDMap(geom)
m1 = WcsNDMap(geom.upsample(2))
coords = m0.geom.get_coord()
m1.fill_by_coord(
tuple([qconcatenate(t, t) for t in coords]),
qconcatenate(coords[1], coords[1]).value,
)
m0.coadd(m1)
assert_allclose(np.nansum(m0.data), np.nansum(m1.data), rtol=1e-4)
