#EXAMPLE# Bethe lattice 2-band Hubbard model, single-site DMFT, Bethe lattice

#####################################
#       To load Sigma(w) 
# put the correct path or '' to avoid loading
# It must be a correct path from the directory 
# where the lisa script is executed
#####################################

Reload_Sigma = False

#####################################
#   General parameters
#####################################

# Maximum number of iterations of the main loop
N_Iter_SelfCons_Max           = 3                          

# Temperature
Beta                          = 6.0                          # 1/Temperature
is_SUPRA = False

#DONT_USE_SHIFT_IN_G0 = True

SCHEME = "DMFT_DOS"
# The number of orbitals is given by the size of this array
Half_Bandwidths = [sqrt(2),sqrt(2)]
AF_Computation = True

#####################################
# Forcing terms
#####################################

Forcing_SC      = [0.0, 0.0]
Forcing_AF      = [0.0, 0.0]
Forcing_F       = [0.0, 0.0]
Forcing_DDW     = [0.0, 0.0]
Number_Iterations_With_Forcing =1

#####################################
# Should we compute at fixed density ?
#####################################

Fix_Density = False	
Chemical_potential_init    = 0.0 + 3*1.5

# Not used if fixdoping=False but required for compilation
Density_Required            =  0.8
Precision_Density           =  0.01
Delta_Chemical_Potential   =  0.1   # MUST be non 0 or it crashes the code !!
Max_loop_mu                =  20
Nbre_Iter_wait_adjust_dop  =   0

 
#####################################

#   Interactions
#   Must be a matrix .....
   
#####################################


import utils
U_interact         = 3.0          # The Coulomb interaction

U_matrix_anti = [[U_interact,U_interact],[U_interact,U_interact]] # utils.ScalarMatrix(2,U_interact)
U_matrix_para = [[0,U_interact],[U_interact,0]] # utils.ScalarMatrix(2,0.0)


	
#####################################

# Choose the solver 

#####################################

Solver = "QMC_Hirsh_Fye"

QMC_Number_Time_Slices            =  16 #int(1.0*Beta* 2.0 )
QMC_N_cycles_MAX              = 100000 #32*1000                       # Maximum number of Monte Carlo iterations
QMC_Compute_Susceptibility        = True

#####################################

# Where do we put the data ?

#####################################

Data_Directory = 'ess'
Data_OverWrite = False