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,
