Fixed call parameters of MfResponse

resolve/likelihood.py

@@ -120,7 +120,7 @@ def ImagingLikelihood(
             R = FullPolResponse(observation, sky_operator.target)
     else:
         if len(sdom.shape) == 3:
-            R = MfResponse(observation, sdom)
+            R = MfResponse(observation, sdom[0], sdom[1])
         else:
             R = StokesIResponse(observation, sdom)
     model_data = R @ sky_operator