#!/usr/bin/env python3
import sys
import py_functions as pf
import matplotlib.pyplot as plt
import numpy as np
arguments = []
for arg in sys.argv[1:]:
arguments.append(arg)
OUTNAME = 0
if len(arguments) >= 1:
OUTNAME = arguments[0]
plt.style.use('seaborn-paper')
fig = plt.figure()
ax1 = fig.add_axes([0.1, 0.51, 0.6, 0.39])
pf.plot_mag_susc_2(ax1, colour='black')
ax1.set_xticklabels([])
ax2 = fig.add_axes([0.1, 0.12, 0.6, 0.39])
pf.plot_entropy_2(ax2, colour='black')
if OUTNAME != 0:
plt.savefig(OUTNAME)
else:
plt.show()
OUTNAME = 0
if len(arguments) >= 1:
OUTNAME = arguments[0]
plt.style.use('seaborn-paper')
fig = plt.figure()
U, G = pf.get_u_gamma()
ax1 = fig.add_axes([0.2, 0.55, 0.55, 0.45])
pf.plot_mag_susc_2(ax1, colour='tab:purple', llabel=f"$U={U}$, $\Gamma={G}$")
ax1.set_xticklabels([])
ax1.set_yticks(np.arange(0, 1, 0.25))
ax2 = fig.add_axes([0.2, 0.15, 0.55, 0.40])
pf.plot_entropy_2(ax2, colour='tab:purple')
ax2.set_yticks(np.arange(0, 4, 1))
ax2.set_xticks(np.exp(np.arange(-20, 1, 4)))
ax2.semilogx()
ax1.legend(loc='upper left',
frameon=True,
facecolor='white',
edgecolor='white',
framealpha=0.8,
fontsize=15)
if OUTNAME != 0:
plt.savefig(OUTNAME)
else:
plt.show()
subprocess.call('./1e_run_th_fl')
N = pf.get_N('e-q0s1_table1.dat')
col = colours[colour]
pf.plot_mag_susc_2(ax_mag,
colour=col,
llabel='$U = {:3.0e}$ , $\Gamma = {:3.0e}$'.format(
pair[0], pair[1]),
llinestyle=lin,
labels=False)
ax_mag.tick_params(labelsize=15)
ax_mag.set_xticklabels([])
pf.plot_entropy_2(ax_entr, colour=col, llinestyle=lin, labels=False)
ax_entr.tick_params(labelsize=15)
ymax = 0.5
par_mag = ax_mag.twinx()
par_mag.set_ylim(ax_mag.get_ylim())
par_mag.set_yticks((0, 0.125, 0.25))
par_mag.set_yticklabels(('$0$ : SC', '$1/8$ : FO', '$1/4$ : LM'))
par_mag.tick_params(axis='both', which='major', labelsize=15)
ymax = mt.log(4) + 0.2
par_entr = ax_entr.twinx()
par_entr.set_ylim(ax_entr.get_ylim())
par_entr.set_yticks((0, mt.log(2), mt.log(4)))
par_entr.set_yticklabels(('$\ln(1)$ : SC', '$\ln(2)$ : LM', '$\ln(4)$ : FO'),
fontsize=15)
if line[0]=='J':
new_param.write( "Jkondo={}\n".format(j) )
else:
new_param.write(line)
param.close()
new_param.close()
subprocess.call([ 'mv' , 'new_param' , 'param' ])
subprocess.call( './1d_run_th_fl' )
N = pf.get_N( 'e-q0s1_table1.dat' )
col = colours[colour]
pf.plot_mag_susc_2( ax_mag , colour=col , mark=mark , llabel=f"$J={j}$" , twin=False )
ax_mag.set_xticklabels([])
pf.plot_entropy_2( ax_entr , colour=col , mark=mark , twin=False )
ymax=1
par = ax_entr.twinx()
par.set_ylim([ 0 , ymax ])
par.set_yticks(( 0 , mt.log(2) ))
par.set_yticklabels(( '$\ln(1)$ : SC' , '$\ln(2)$ : LM' ) )
par.tick_params( labelsize='large' )
ymax=0.5
par = ax_mag.twinx()
par.set_ylim([ 0 , ymax ])
par.set_yticks(( 0 , 0.25 ))
par.set_yticklabels(('$0$ : SC','$1/4$ : LM'))
par.tick_params( labelsize='large' )
for line in param:
if line=="" or line=="\n":
new_param.write("\n")
else:
if line[0]=='U':
new_param.write("U={}\n".format(pair[0]))
elif line[:5]=='Gamma':
new_param.write('Gamma={}\n'.format(pair[1]))
else:
new_param.write(line)
param.close()
new_param.close()
subprocess.call([ 'mv' , 'new_param' , 'param' ])
subprocess.call( './1d_run_th_fl' )
N = pf.get_N( 'e-q0s2_table1.dat' )
col = colours[colour]
pf.plot_mag_susc_2( ax_mag , colour=col , llabel='$U = {:3.0e}$ , $\Gamma = {:3.0e}$'.format(pair[0],pair[1]) )
pf.plot_entropy_2( ax_entr , colour=col )
ax_mag.legend( loc='upper left' , frameon=True , facecolor='white' , edgecolor='white' , framealpha=0.8 , ncol=2 )
if OUTFILE!=0:
plt.savefig("/graphs_and_data/arrowplots_multiple/{}".format(OUTFILE))
else:
plt.show()
arguments.append( arg )
OUTNAME = 0
if len(arguments)>=1:
OUTNAME = arguments[0]
plt.style.use( 'seaborn-paper' )
fig = plt.figure()
#colour = list( mcolors.TABLEAU_COLORS.values() )[0]
colour = 'black'
ax1 = fig.add_axes( [0.2, 0.55, 0.6, 0.43] )
pf.plot_mag_susc_2( ax1 , colour=colour , llabel=f'$U = {pf.get_u_gamma()[0]}$, $\Gamma = {pf.get_u_gamma()[1]}$' )
ax1.set_xticklabels([])
ax1.tick_params( axis='both' , which='major' , labelsize=15 )
ax2 = fig.add_axes( [0.2, 0.2, 0.6, 0.35] )
pf.plot_entropy_2( ax2 , colour=colour )
plt.xticks( fontsize=15 )
plt.yticks( fontsize=15 )
ax2.tick_params( axis='both' , which='major' , labelsize=15 )
ax1.legend( loc='upper left' , frameon=True , framealpha=0.7 , facecolor='white' , edgecolor='white' , fontsize=15 )
if OUTNAME!=0:
plt.savefig( OUTNAME )
else:
plt.show()