ann/plot-eps-scale.py at master · abhirathb/ann · GitHub

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
import numpy as np
import matplotlib.pyplot  as plt

prefix = 'result_22-8/eps-scale/'
suffix = '_1000_100'

x = np.array([0,1,2,3,4])
xticks = ['0.1','0.01','0.001','0.0001','0.00001']

h=np.loadtxt(prefix+'ham'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,h,"o")
plt.title('Hamiltonian')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'ham' + suffix +'.png')

pot=np.loadtxt(prefix+'pot'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,pot,"o")
plt.title('Potential')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'pot' + suffix +'.png')

kin=np.loadtxt(prefix+'kin'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,kin,"o")
plt.title('Kinetic Energy')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'kin' + suffix +'.png')

diff=np.loadtxt(prefix+'diff'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,diff,"o")
plt.title('Difference in Hamiltonian')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'diff' + suffix +'.png')


diff_k=np.loadtxt(prefix+'diff-k'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,diff_k,"o")
plt.title('Difference in Kinetic Energy')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'diff-k' + suffix +'.png')

diff_u=np.loadtxt(prefix+'diff-u'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,diff_u,"o")
plt.title('Difference in Potential Energy')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'diff-u' + suffix +'.png')


ratio=np.loadtxt(prefix+'ratio'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,ratio,"o")
plt.title('DelH/H')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'ratio' + suffix +'.png')

ratio_k=np.loadtxt(prefix+'ratio-k'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,ratio_k,"o")
plt.title('DelK/K')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'ratio-k' + suffix +'.png')


ratio_u=np.loadtxt(prefix+'ratio-u'+suffix)
plt.figure()
plt.xticks(x,xticks)
plt.plot(x,ratio_u,"o")
plt.title('DelU/U')
plt.ylabel('Energy')
plt.xlabel('value of epsilon')
plt.savefig(prefix + 'ratio-u' + suffix +'.png')