To receive notifications about scheduled maintenance, please subscribe to the mailing-list gitlab-operations@sympa.ethz.ch. You can subscribe to the mailing-list at https://sympa.ethz.ch

Commit 7282acb3 authored by mgassner's avatar mgassner
Browse files

minor

parent 29a39373
......@@ -52,7 +52,6 @@ def MI_vs_Entropy(model_name, number_spins, param_range):
model_name = model_name
number_spins = number_spins
periodic = False
spin_inversion = None
params = np.linspace(param_range[0], param_range[1], int(param_range[2]))
sub_dim = 4
......@@ -69,7 +68,7 @@ def MI_vs_Entropy(model_name, number_spins, param_range):
else:
hamming_weight = number_spins // 2
model = spin_model(model_name=model_name, number_spins=number_spins, periodic=periodic,
param=param, hamming_weight=hamming_weight ,spin_inversion=spin_inversion)
param=param, hamming_weight=hamming_weight)
model.compute_ew_and_ev()
basis_states = model.basis.states
gs = model.eigenstates[:,0]
......
......@@ -19,7 +19,6 @@ def test_xxz():
model_name = 'xxz'
number_spins = 16
periodic = False
spin_inversion = None
params = np.linspace(-1.2, 3.0, 41)
entropies = []
for param in params:
......@@ -28,7 +27,7 @@ def test_xxz():
else:
hamming_weight = number_spins // 2
model = spin_model(model_name=model_name, number_spins=number_spins, periodic=periodic,
param=param, hamming_weight=hamming_weight, use_symmetries=False ,spin_inversion=spin_inversion)
param=param, hamming_weight=hamming_weight)
model.compute_ew_and_ev()
print('EIGENSTATE at h/J = ', param, 'is: ', model.eigenstates[:,0])
basis_states = model.basis.states
......@@ -36,7 +35,7 @@ def test_xxz():
print('Number Spins: ', model.basis.number_spins)
print('States', model.basis.states)
sub_dim = number_spins // 2
rhos = reduced_dm(sub_dim, number_spins, hamming_weight, gs, basis_states, spin_inversion)
rhos = reduced_dm(sub_dim, number_spins, hamming_weight, gs, basis_states)
entropy = compute_entropy(rhos)
entropies.append(entropy)
......@@ -48,7 +47,7 @@ def test_xxz():
plt.title('Model:' + model_name + ' Number spins:' + str(number_spins))#, ' Periodic:' + str(periodic))
plt.grid(True)
plt.savefig('output/' + 'entropy' + timestr + '.jpg')
plt.savefig('output/' + 'entropytest.jpg')
def test_area_law(number_spins, periodic, param_range):
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment