### args['theta'] = theta args['a1_x'] = a1_x args['a1_y'] = a1_y args['a2_x'] = a2_x args['a2_y'] = a2_y args['N1'] = N1 args['N2'] = N2 args['N_spins'] = N1 * N2 args['first_neighb'] = inputs['first_neighb'][0] ## ## ### start_time = timeit.default_timer() print 'building lattice :' neighbors_table = constructLattice(**args) args['neighbors_table'] = neighbors_table elapsed = timeit.default_timer() - start_time print 'time elapsed: %s ' % elapsed print '------' ############ MC specs args['max_cluster_size'] = int(inputs['max_cluster_size']) #1 ############ callibration: ### ##log_f = 1 ##f_factor= np.exp(log_f) args['f_factor_init'] = float(inputs['f_factor_init']) ## ##precision_f = 4 args['precision_f'] = int(inputs['precision_f']) #precision_f ##
a1_x = 1.0
a1_y = 0
theta = np.pi / 3
a2_x = np.cos(theta)
a2_y = np.sin(theta)
#
args['a1_x'] = a1_x
args['a1_y'] = a1_y
args['a2_x'] = a2_x
args['a2_y'] = a2_y
#########
#
args['N1'] = N1
args['N2'] = N2
args['N_spins'] = N1 * N2
args['first_neighb'] = True
neighbors_tables_list = constructLattice(**args)
neighbors_table = constructLattice(**args)
args['neighbors_table'] = neighbors_table
E_s = []
for i in range(2**(N1 * N2)):
config = getBinaryArray(N1 * N2, i)
config = [(-1)**config[i] for i in range(len(config))]
E = getEnergyOfConfig(config, **args)
E_s.append(E)
# print config, E
print np.min(E_s), np.max(E_s)
plt.hist(E_s, 100)
#### 4 * 4 ===> -21.1949051216 48.5294289008