avoid changing map in map2alm

carronj · Oct 17, 2023 · 2044a2a · 2044a2a
1 parent 09d715b
commit 2044a2a
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/lenspyx/remapping/utils_geom.py b/lenspyx/remapping/utils_geom.py
@@ -183,7 +183,7 @@ def alm2map_spin(self, gclm:np.ndarray, spin:int, lmax:int, mmax:int, nthreads:i
     def map2alm_spin(self, m:np.ndarray, spin:int, lmax:int, mmax:int, nthreads:int, zbounds=(-1., 1.), **kwargs):
         # FIXME: method only here for backwards compatiblity
         assert zbounds[0] == -1 and zbounds[1] == 1., zbounds
-        return self.adjoint_synthesis(m, spin, lmax, mmax, nthreads, **kwargs)
+        return self.adjoint_synthesis(m.copy(), spin, lmax, mmax, nthreads, **kwargs)
 
     def alm2map(self, gclm:np.ndarray, lmax:int, mmax:int, nthreads:int, zbounds=(-1., 1.), **kwargs):
         # FIXME: method only here for backwards compatiblity
@@ -193,7 +193,7 @@ def alm2map(self, gclm:np.ndarray, lmax:int, mmax:int, nthreads:int, zbounds=(-1
     def map2alm(self, m:np.ndarray, lmax:int, mmax:int, nthreads:int, zbounds=(-1., 1.), **kwargs):
         # FIXME: method only here for backwards compatiblity
         assert zbounds[0] == -1 and zbounds[1] == 1., zbounds
-        return self.adjoint_synthesis(m, 0, lmax, mmax, nthreads, **kwargs).squeeze()
+        return self.adjoint_synthesis(m.copy(), 0, lmax, mmax, nthreads, **kwargs).squeeze()
 
     @staticmethod
     def rings2pix(geom:Geom, rings:np.ndarray[int]):