P_M_masked tested, should properly apply the mask in Fourier space now

proxtoolbox/ProxOperators/P_M.py

@@ -69,6 +69,7 @@ class P_M_masked(P_M):
         -------
         array_like - p_M: the projection IN THE PHYSICAL (time) DOMAIN
         """
-        constrained = super(P_M_masked, self).work(u)
+        fourier_space_iterate = np.fft.fft2(u, axes=None)
+        constrained = magproj(fourier_space_iterate, self.M)
         update = np.where(self.mask, u, constrained)
-        return update
+        return np.fft.ifft2(update, axes=None)