Example #1
0
#!/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()
Example #3
0
    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)
Example #4
0
            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' )
Example #5
0
    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()