From 350119b12ec8fe7c9ef6ac3c3636138f67875f2e Mon Sep 17 00:00:00 2001
From: Pauli Virtanen <pauli.t.virtanen@jyu.fi>
Date: Thu, 25 Apr 2024 12:33:06 +0300
Subject: [PATCH] tests: test self-consistency with h
---
tests/test_solver.py | 20 ++++++++++++++++++++
1 file changed, 20 insertions(+)
diff --git a/tests/test_solver.py b/tests/test_solver.py
index ab8ec4b..0568812 100644
--- a/tests/test_solver.py
+++ b/tests/test_solver.py
@@ -347,3 +347,23 @@ def test_soc_analytic(T, n, alpha_soc):
Js.append(J)
Jans.append(Jy_an)
+
+def test_selfcons_h():
+ T = 0.4
+ h = 0.5
+
+ solver = Solver(nx=5, ny=2)
+ solver.mask[...] = MASK_NONE
+ solver.Lx = 5
+ solver.Ly = 5
+ solver.eta = 0
+ solver.Delta[...] = 1.0
+ solver.Omega[...] += h * np.kron(S_z, S_y)
+
+ T_c0 = np.exp(np.euler_gamma) / np.pi
+ Delta, I0, _, success = solver.self_consistency(T=T, T_c0=T_c0, workers=1)
+
+ Delta0 = BCS_Delta(T, T_c0, h=h)
+ Deltam = np.broadcast_to((Delta0 * S_0)[None, None, :, :], Delta.shape)
+
+ assert_allclose(Delta.copy(), Deltam, rtol=1e-3)
--
GitLab