t0=0.52, U=1, pbc=True, field=32.9, F0=10, a=4, a_scale=1, J_scale=1, tracking=1) """generating our class of scaled parameters""" """this is used for most of the calculations""" lat = hhg(field=param.field, nup=param.N_up, ndown=param.N_down, nx=param.L, ny=0, U=param.U, t=param.t0, F0=param.F0, a=param.a, pbc=param.pbc) """setup our evolution time parameters""" t_p = time_evolution_params(perimeter_params=lat, cycles=2, nsteps=int(2e3), plotting=1) """prepare to load our data to be plotted""" outfile = './Data/expectations:{}sites-{}up-{}down-{}t0-{}U-{}cycles-{}steps-{}pbc:a_scale={:.2f}-J_scale={:.2f}.npz'.\ format(param.L, param.N_up, param.N_down, param.t0, param.U, t_p.cycles, t_p.n_steps, param.pbc, param.a_scale , param.J_scale) expectations = np.load(outfile) """plot out expectations"""
param = unscaledparam(L=6, t0=0.52, U=0.5, pbc=True, field=32.9, F0=9, a=4, a_scale=1, J_scale=1, tracking=1) lat = hhg(field=param.field, nup=param.N_up, ndown=param.N_down, nx=param.L, ny=0, U=param.U, t=param.t0, F0=param.F0, a=param.a, pbc=param.pbc) acceptable_current = UnivariateSpline(np.linspace(0, 2 / lat.freq, num=len(a_d['current'])), a_d['current'], s=0) param_track = unscaledparam(L=6, t0=0.52, U=0.5, pbc=True,