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 ec58be1a authored by mgassner's avatar mgassner
Browse files

minor changes

parent e4eecabc
......@@ -15,18 +15,35 @@ parser.add_argument("--boundary_condition", "-b", dest='periodic' ,type=bool, re
args = parser.parse_args()
def test_xxz(number_spins, param_range, periodic=False):
def test_xxz(number_spins, param_range, periodic=True):
model_name = 'xxz'
params = np.linspace(param_range[0], param_range[1], int(param_range[2]))
entropies = []
hamming_weight = number_spins // 2
hamming_weight = None #number_spins // 2
for param in params:
if param > -0.95:
hamming_weight = number_spins // 2
model = spin_model(model_name=model_name, number_spins=number_spins, periodic=periodic,
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
gs = model.eigenstates[:,0]
# just for testing
if number_spins == 16 and param <= -1.1 and param > -1.4:
print(gs[0])
printi = []
for i, x in enumerate(gs):
if abs(x) > 1e-1:
printi.append([x, basis_states[i]])
print(printi)
if number_spins == 10 and param < -2.0:
print(gs[0])
printi = []
for i, x in enumerate(gs):
if abs(x) > 1e-2:
printi.append([x, basis_states[i]])
print(printi)
#print('Number Spins: ', model.basis.number_spins)
#print('States', model.basis.states)
sub_dim = number_spins // 2
......
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