From c70e1aa48359ad5b2a6d06406e3027709ae5b8c8 Mon Sep 17 00:00:00 2001
From: 70akaline <36015241+70akaline@users.noreply.github.com>
Date: Thu, 26 Nov 2020 13:02:33 +0000
Subject: [PATCH] wrong order of orbital setting and set Coulomb mat

Construct Coulomb matrix needs prior to defining the orbital
---
 doc/tutorials/images_scripts/dft_dmft_cthyb.py | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/doc/tutorials/images_scripts/dft_dmft_cthyb.py b/doc/tutorials/images_scripts/dft_dmft_cthyb.py
index 9d8ae67b..b2744bbe 100644
--- a/doc/tutorials/images_scripts/dft_dmft_cthyb.py
+++ b/doc/tutorials/images_scripts/dft_dmft_cthyb.py
@@ -19,10 +19,6 @@ h_field = 0.0
 U = 4.0
 J = 0.65
 dc_type = 1                      # DC type: 0 FLL, 1 Held, 2 AMF
-# Construct U matrix for density-density calculations
-Umat, Upmat = U_matrix_kanamori(n_orb=n_orb, U_int=U, J_hund=J)
-# Construct density-density Hamiltonian
-h_int = h_int_density(spin_names, orb_names, map_operator_structure=SK.sumk_to_solver[0], U=Umat, Uprime=Upmat)
 
 ## SLATER HAMILTONIAN
 ## Define interaction paramters, DC and Hamiltonian
@@ -74,6 +70,10 @@ n_orb = SK.corr_shells[0]['dim']
 l = SK.corr_shells[0]['l']
 spin_names = ["up","down"]
 orb_names = [i for i in range(n_orb)]
+# Construct U matrix for density-density calculations
+Umat, Upmat = U_matrix_kanamori(n_orb=n_orb, U_int=U, J_hund=J)
+# Construct density-density Hamiltonian
+h_int = h_int_density(spin_names, orb_names, map_operator_structure=SK.sumk_to_solver[0], U=Umat, Uprime=Upmat)
 
 # Use GF structure determined by DFT blocks
 gf_struct = [(block, indices) for block, indices in SK.gf_struct_solver[0].items()]