core: fix sign in terminal GF

usadelndsoc/core.py

@@ -56,6 +56,7 @@ class Core:
         for i, j in zip(*np.where(m)):
             M = self._Omega[i, j] + 1j * self._omega * tau3
             Q = signm(-1j * M)
+            Q = Q * np.sign(self._omega)
             g = np.linalg.solve(np.eye(2) + Q[:2, :2], Q[:2, 2:])
             gt = np.linalg.solve(np.eye(2) - Q[2:, 2:], Q[2:, :2])
             Phi[i, j, 0] = g