arguments = []
for arg in sys.argv[1:]:
arguments.append(arg)
OUTFILE = 0
if np.size(arguments) >= 1:
OUTFILE = arguments[1]
plt.style.use('seaborn-dark-palette')
fig, ax = plt.subplots()
plt.grid(b=True, which='major', color='black', linewidth='0.4', alpha=0.3)
N = pf.get_N('e-q0s1_table1.dat')
N_array = np.arange(0, N, 2)
J_lm = np.zeros(len(N_array))
pf.plot_arrows(N_array, J_lm, ax, colour='tab:blue', llabel='local moment')
J_sc = 1.33742913 * np.ones(len(N_array))
pf.plot_arrows(N_array, J_sc, ax, colour='tab:red', llabel='strong coupling')
pf.plot_arrow_diagram(ax)
plt.legend(loc="upper center",
frameon=True,
framealpha=0.7,
facecolor='white',
edgecolor='white')
if OUTFILE != 0:
plt.savefig("{}".format(OUTFILE))
else:
nrg_eigenvalues = pf.read_h0_eigenvalues_lowest()
U_array = np.empty( int(N/2) )
GAMMA_array = np.empty( int(N/2) )
nrg_eigenvalue_step = np.empty( 6 )
for n in range(0,int(N/2)):
nrg_eigenvalues_step = nrg_eigenvalues[n,:]
res = opt.least_squares( F , [0.3,0.7] )
U_array[n] = res.x[0]
GAMMA_array[n] = res.x[1]
col = colours[colour]
pf.plot_arrows( U_array , GAMMA_array , ax , colour='black' , fixedpoints=False )
#ax.plot(GAMMA_array,U_array,color=col, linewidth=0 ,\
# label='${}={:3.0e}$, ${}={:3.0e}$'.format('U',pair[0],'\Gamma',pair[1]))
pf.include_fixedpoints_siam( ax , multi=True )
ax.set_yticks( np.arange(-9,8,1) )
legend = ax.legend(loc="upper right",frameon=True,framealpha=0.7,facecolor='white',edgecolor='white')
for handle in legend.legendHandles:
handle.set_linewidth(3.0)
if OUTFILE!=0:
plt.savefig("/graphs_and_data/arrowplots_multiple/{}".format(OUTFILE))
subprocess.call(['mv', 'new_param', 'param'])
subprocess.call('./1d_run_th_fl')
N = pf.get_N('e-q0s2_table1.dat')
nrg_eigenvalues = pf.read_h0_eigenvalues_lowest()
U_array = np.empty(int(N / 2))
GAMMA_array = np.empty(int(N / 2))
nrg_eigenvalue_step = np.empty(6)
for n in range(0, int(N / 2)):
nrg_eigenvalues_step = nrg_eigenvalues[n, :]
res = opt.least_squares(F, [0.3, 0.3], bounds=([-100, 0], [100, 100]))
U_array[n] = res.x[0]
GAMMA_array[n] = res.x[1]
pf.plot_arrows(U_array, GAMMA_array, ax, colour='black')
pf.include_fixedpoints(ax)
ax.set_xlim([-0.2, 1.2])
ax.set_xticks(np.arange(0, 1.3, 0.2))
ax.set_ylim([-3, 6])
ax.set_yticks(np.arange(-3, 5.5, 1))
ax.tick_params(labelsize=15)
if OUTFILE != 0:
plt.savefig(OUTFILE)
else:
plt.show()