Add another test

test/test_operators/test_correlated_fields.py

@@ -36,7 +36,6 @@ def _stats(op, samples):
                 ift.HPSpace(8), ift.GLSpace(4)])
 @pmp('N', [0, 2])
 def testAmplitudesInvariants(sspace, N):
-
     fsspace = ift.RGSpace((12,), (0.4,))
     dofdex1, dofdex2, dofdex3 = None, None, None
     if N == 2:
@@ -55,6 +54,8 @@ def testAmplitudesInvariants(sspace, N):
     for ampl in fa.normalized_amplitudes:
         ift.extra.check_jacobian_consistency(ampl, ift.from_random(ampl.domain),
                                              ntries=10)
+    ift.extra.check_jacobian_consistency(op, ift.from_random(op.domain),
+                                         ntries=10)
 
     samples = [ift.from_random(op.domain) for _ in range(100)]
     tot_flm, _ = _stats(fa.total_fluctuation, samples)