#!/usr/bin/env looper
#PRELUDE: $Nz=2;
#AUTOLOOP: nrginit ; nrgrun
#OVERWRITE

[sweep]
Nz=2

[extra]
# We need to define epsilon here. For historical reasons, only
# U, Gamma, delta, and t may appear in the main [param] section,
# all other parameters need to be specified in the section [extra].
# epsilon = delta-U/2
epsilon=-0.15

[param]
symtype=QS
Lambda=4.0
Tmin=1e-8
keepmin=200
keepenergy=10.0
keep=10000

discretization=Z
@$z = 1/$Nz; $z <= 1.00001; $z += 1/$Nz
z=$z

# Alternative definition of SIAM in terms of epsilon, U, Gamma
model=../model.m
U=0.5
Gamma=0.03
#delta=0.1  -- we use epsilon instead!

ops=A_d self_d
specd=A_d-A_d self_d-A_d

fdm=true

# Broadening is performed by an external tool
broaden_max=2
broaden_ratio=1.01
broaden_min=1e-6
bins=300
broaden=false
savebins=true

T=1e-10
