From e5b1c88dd19133c30fb12080d5cebf7e7c114e8b Mon Sep 17 00:00:00 2001
From: Irene Abril Cabezas <79768254+iabrilcabezas@users.noreply.github.com>
Date: Sun, 10 Aug 2025 16:31:18 +0100
Subject: [PATCH] fix to "okmap becomes complex164 if left unattended"

---
 mnms/utils.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/mnms/utils.py b/mnms/utils.py
index b590029..5cf82fe 100644
--- a/mnms/utils.py
+++ b/mnms/utils.py
@@ -1383,7 +1383,7 @@ def fourier_resample(imap, omap=None, shape=None, wcs=None, dtype=None):
 
     kx = np.fft.rfftfreq(omap.shape[-1]).astype(dtype, copy=False)
     ky = np.fft.fftfreq(omap.shape[-2])[..., None].astype(dtype, copy=False)
-    phase = np.exp(-2j*np.pi*(ky*shift[0] + kx*shift[1]))
+    phase = np.exp(-2j*np.pi*(ky*shift[0] + kx*shift[1])).astype(okmap.dtype, copy=False)
 
     okmap = concurrent_op(np.multiply, okmap, phase)
     okmap = enmap.ndmap(okmap, omap.wcs)