minor update

src/dcore/impurity_solvers/hphi_spectrum.py

@@ -124,6 +124,14 @@ def calc_one_body_green(one_body_green_core):
             G[i, j] = G_ij[i, j]
             G[j, i] = G_ji[i, j]
 
+    # for i in range(G_shape[0]):
+    #     for j in range(i+1, G_shape[1]):  # i<j
+    #         g_iijj = G_diag[i] - G_diag[j]
+    #         A2 = A[i, j] - g_iijj
+    #         B2 = B[i, j] - g_iijj
+    #         G[i, j] = (A2 - 1J * B2) / 2.
+    #         G[j, i] = (A2 + 1J * B2) / 2.
+
     return G.reshape(n_site, n_sigma, n_site, n_sigma, n_T, n_omega)
 
 

src/dcore/impurity_solvers/pomerol.py

@@ -83,14 +83,16 @@ def assign_from_numpy_array(g_block, data, block_names):
     for b in range(data.shape[0]):
         for o1 in range(data.shape[1]):
             for o2 in range(data.shape[2]):
-                ax.plot(data[b, o1, o2, :30].real, marker='.')
+                ax.plot(data[b, o1, o2, :30].real, marker='.', label=f"{b}, ({o1},{o2})")
+    ax.legend()
     fig.savefig("test_gf_re.pdf")
 
     fig, ax = plt.subplots()
     for b in range(data.shape[0]):
         for o1 in range(data.shape[1]):
             for o2 in range(data.shape[2]):
-                ax.plot(data[b, o1, o2, :30].imag, marker='.')
+                ax.plot(data[b, o1, o2, :30].imag, marker='.', label=f"{b}, ({o1},{o2})")
+    ax.legend()
     fig.savefig("test_gf_im.pdf")