# title += "<br>"
# title += "<br>w<sub>c</sub> = " + wc_str(config.wc)
# title += "<br>w<sub>a</sub> = " + wa_str(config.wa)
# title += "<br>g = " + g_str(config.g)
# title += "<br>γ = " + g_str(config.l)
# title += "</b>"
# print(config.path)
if config.T / 1e-6 >= 0.5:
PyPlot3D(
title=title,
z_csv=config.path + "/" + "z.csv",
x_csv=config.path + "/" + "x.csv",
y_csv=config.path + "/" + "t.csv",
# t_coeff=20000 / 1000 * (config.T / 1e-6),
# t_coeff=20000 / 1000 * (config.T / 1e-6),
t_coeff=1,
online=False,
path=config.path,
filename="Bipartite",
xaxis="states",
yaxis="time, " + T_str_mark(config.T),
y_scale=y_scale)
else:
PyPlot3D(
title=title,
z_csv=config.path + "/" + "z.csv",
x_csv=config.path + "/" + "x.csv",
y_csv=config.path + "/" + "t.csv",
# t_coeff=20000 / 1000 * (config.T / 1e-6),
t_coeff=1,
online=False,
run(ro_0=ro_0, H=H, dt=config.dt, nt=config.nt, config=config)
# -------------------------------------------------------------------------------------------------
from Common.STR import *
if not __debug__:
title = "<b>"
# title += "<br>wc = " + str(wc / GHz) + " GHz" # title += "<br>wa = "
# + str(np.array(wa) / GHz) + " GHz" # title += "<br>g = " +
# str(np.array(g) / MHz) + " MHz"
title = "<b>Bipartite oscillations</b>"
title += "<br>capacity = " + \
str(config.capacity) + ", n_1 = " + \
str(config.n_1) + ", n_2 = " + str(config.n_2)
title += "<br>t = " + T_str(config.T)
title += "</b>"
PyPlot3D(
title=title,
z_csv=config.path + "/" + "z.csv",
x_csv=config.path + "/" + "x.csv",
y_csv=config.path + "/" + "t.csv",
online=False,
path=config.path,
filename="Bipartite2",
xaxis="states",
yaxis="time, мкс",
y_scale=5,
)
# -------------------------------------------------------------------------------------------------
# if __debug__:
# w_0.print()
# -------------------------------------------------------------------------------------------------
ro_0 = DensityMatrix(w_0)
# if __debug__:
# ro_0.write_to_file(filename=ro_0_csv)
# -------------------------------------------------------------------------------------------------
run_ro(ro_0, H, dt=dt, nt=nt)
# -------------------------------------------------------------------------------------------------
if not __debug__:
title = "<b>Bipartite oscillations</b>" + "<br>capacity = " + str(
capacity) + ", n = " + str(n) + ", t = " + str(T / mks) + " мкс"
title += "<br>wc = " + str(wc / GHz) + " GHz"
title += "<br>wa = " + str(np.array(wa) / GHz) + " GHz"
title += "<br>g = " + str(np.array(g) / MHz) + " MHz"
PyPlot3D(title=title,
z_csv=path + "/" + "z.csv",
x_csv=path + "/" + "x.csv",
y_csv=path + "/" + "t.csv",
online=False,
path=path,
filename="Bipartite",
xaxis="states",
yaxis="time, мкс")
# -------------------------------------------------------------------------------------------------
elif config.T == 0.5 * config.mks:
y_scale = 0.01
elif config.T == 1 * config.mks:
y_scale = 7.5
# y_scale = 10
elif config.T == 5 * config.mks:
y_scale = 1
if not __debug__:
title = "<b>"
title += "capacity = " + str(config.capacity) + ", n = " + str(config.n)
title += "<br>w<sub>c</sub> = " + wc_str(config.wc)
title += "<br>w<sub>a</sub> = " + \
"[" + ", ".join([wa_str(i) for i in config.wa]) + "]"
title += "<br>g = " + "[" + ", ".join([g_str(i) for i in config.g]) + "]"
title += "<br>t = " + T_str(config.T)
title += "<br>l = " + wc_str(config.l)
title += "</b>"
PyPlot3D(title=title,
z_csv=config.path + "/" + "z.csv",
x_csv=config.path + "/" + "x.csv",
y_csv=config.path + "/" + "t.csv",
online=False,
path=config.path,
filename="BipartiteGeneralLindblad",
xaxis="states",
yaxis="time, " + T_str_mark(config.T),
y_scale=y_scale)
# -------------------------------------------------------------------------------------------------