#!/usr/bin/env looper # Hubbard model # Rok Zitko, zitko@theorie.physik.uni-goettingen.de, Oct 2008 # #PRELUDE: $Nz=4; $DELTA=`cat param.delta`; #LOOP: mkdir -p $cnt ; cd $cnt ; cat - >param ; odesolvasymif1 ; nrginit ; nrgrun ; touch DONE ; cd .. [dmft] Nz=4 scale=4 clip=1e-5 offset=1 alpha=1 [param] xmax=30 adapt=false tri=cpp preccpp=2000 prec=100 symtype=QS model=SIAM U=0.6 delta=$DELTA Gamma=-999 numthreads=1 options=WRITE READHAM READBASIS writedir=../ Lambda=2. Tmin=1e-10 keepmin=500 keepenergy=10.0 keep=10000 safeguard=0.001 discretization=Z band=asymode dos=../Delta.dat @$z = 1.0/$Nz; $z <= 1.00001; $z += 1/$Nz z=$z !print "$cnt $z\n"; strategy=kept ops=A_d specd=A_d-A_d self=true specb=A_d-A_d prec_xy=12 verbose=true mmadebug=2 dmnrg=true goodE=2.3 NN2avg=true smooth=new omega0=1e-99 alpha=0.2 bins=1000 T=1e-99 broaden_max=4 broaden_ratio=1.01 broaden_min=1e-8