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)