-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathjob_dos.bash
103 lines (94 loc) · 2.52 KB
/
job_dos.bash
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
#rm -fr `find .|grep -v "\<job.bash\>"`
job_id=2
U=0.1
#Ed_up=-0.076
#Ed_down=-0.076
Lambda=2
alpha=0.69
num_kept=512
smear=0
unsmear=1
dim_imp=4
dim_dot=4
Beta_bar=0.6
Q=0
Q_Sz=0
N_up_N_down=1
occupation=1
N_max=0
Nz=4
dir=`pwd`
Vg=0
h=0
for Gamma0_U in 7 #3 4 5 6 7 8 9
do
Gamma0=`awk "BEGIN {print $Gamma0_U*$U}"`
#Gamma0=`echo "${Gamma0_U}*$U"|bc -l`
for delta_U in -0.2 -0.22 -0.24 -0.26 -0.28 -0.3 -0.32 -0.34 -0.36 -0.38 -0.4 -0.42
#for delta_U in -0.2 -0.22 -0.24 -0.25 -0.27 -0.28 -0.29 -0.3 -0.32 -0.34 -0.36 -0.38 -0.4 -0.42
do
Ed_up=`awk "BEGIN {print ${delta_U}*$U-$U/2.0}"`
Ed_down=${Ed_up}
echo $Gamma0 $Ed_up
mkdir Gamma0_U${Gamma0_U}_delta_U$delta_U
cd Gamma0_U${Gamma0_U}_delta_U$delta_U
dir1=`pwd`
for temperature in 1e-9
do
omega0=$temperature
#rm -fr $temperature
mkdir $temperature
cd $temperature
for z in 0 0.25 0.5 0.75
#for z in 0.5 0.75
do
#rm -fr $z
mkdir $z
cd $z
cp /home/ligy/NRG/fdmnrg.x .
echo $job_id > job_id
echo $U >> input_total
echo $Ed_up >> input_total
echo $Ed_down >> input_total
echo $temperature >> input_total
echo $Lambda >> input_total
echo $alpha >> input_total
echo $num_kept >> input_total
echo $smear >> input_total
echo $unsmear >> input_total
echo $omega0 >> input_total
echo $dim_imp >> input_total
echo $dim_dot >> input_total
echo $Beta_bar >> input_total
echo $Q >> input_total
echo $Q_Sz >> input_total
echo $N_up_N_down >> input_total
echo $occupation >> input_total
echo $N_max >> input_total
echo $temperature >> input_band
echo $Lambda >> input_band
echo $num_kept >> input_band
echo $dim_imp >> input_band
echo $dim_dot >> input_band
echo $Beta_bar >> input_band
echo $Q >> input_band
echo $Q_Sz >> input_band
echo $N_up_N_down >> input_band
echo $N_max >> input_band
cp /home/ligy/NRG/makefreq.cpp .
cp /home/ligy/NRG/work_dir/DOS_U=0.1/test.submit .
icpc makefreq.cpp -o makefreq.x
./makefreq.x
cp ~/data_chain_fitting/chain_Lambda${Lambda}_r1_k${Gamma0}_Vg${Vg}_h${h}_z${z}.dat chain_band.dat
cp ~/data_chain_fitting/chain_Lambda${Lambda}_r1_k${Gamma0}_Vg${Vg}_h${h}_z${z}.dat chain_total.dat
job_name=dos-${delta_U}_${Gamma0_U}-$z
qsub -N ${job_name} test.submit
cd ..
done
cd ..
done
cd ..
done
cd ..
done
exit 0