Make sure that target domain is actually the correct one

nifty5/library/correlated_fields.py

@@ -55,7 +55,7 @@ def CorrelatedField(target, amplitude_operator, name='xi', codomain=None):
     if codomain is None:
         codomain = tgt[0].get_default_codomain()
     h_space = codomain
-    ht = HarmonicTransformOperator(h_space, tgt[0])
+    ht = HarmonicTransformOperator(h_space, target=tgt[0])
     p_space = amplitude_operator.target[0]
     power_distributor = PowerDistributor(h_space, p_space)
     A = power_distributor(amplitude_operator)
@@ -92,7 +92,7 @@ def MfCorrelatedField(target, amplitudes, name='xi'):
 
     hsp = DomainTuple.make([tt.get_default_codomain() for tt in tgt])
     ht1 = HarmonicTransformOperator(hsp, target=tgt[0], space=0)
-    ht2 = HarmonicTransformOperator(ht1.target, space=1)
+    ht2 = HarmonicTransformOperator(ht1.target, target=tgt[1], space=1)
     ht = ht2 @ ht1
 
     psp = [aa.target[0] for aa in amplitudes]