core: fix prefactors of i in S_H

examples/bulk_j_an.py

@@ -86,17 +86,15 @@ def j_anom(T, h, n=5, alpha_soc=0.01, perturbative=False):
         J[:, :, UP] /= dd[0]
         J[:, :, RIGHT] /= dd[1]
         J[:, :, DOWN] /= dd[0]
-        Jtot += pi * aa * J
+        Jtot += 1j * pi * aa * J
 
         # Analytic
         gp = (ww - 1j * h) / np.sqrt((ww - 1j * h) ** 2 + 1)
         gm = (ww + 1j * h) / np.sqrt((ww + 1j * h) ** 2 + 1)
         fp = 1 / np.sqrt((ww - 1j * h) ** 2 + 1)
         fm = 1 / np.sqrt((ww + 1j * h) ** 2 + 1)
-        # XXX: check the 32*eta
-        Jy_an = (
-            (gp - gm) * (1 + gp * gm + fp * fm) * (-1j * alpha_soc**3) * (32 * eta)
-        )
+        # XXX: check the 8*eta
+        Jy_an = (gp - gm) * (1 + gp * gm + fp * fm) * (-1j * alpha_soc**3) * (8 * eta)
         Jantot += pi * aa * Jy_an
 
     return Jtot, Jantot

src/action.hpp

@@ -251,7 +251,7 @@ inline Scalar S_H(const array::ArrayView<const Mask, Shape0>& mask,
         }
     }
 
-    return (2.0*eta/(dd[0]*dd[1])) * S;
+    return (2.0*eta/(Complex(0,4.0)*dd[0]*dd[1])) * S;
 }
 
 /* Function that has the same value as the AD variable,