# Load the data from the dsp file:
path = r'/Volumes/slab/Gerwin/Electron on helium/Maxwell/M018 Yggdrasil/M018V6/With L3 etch - just trap'
potential_file = r"ResonatorBiasSymmetricPotential_1100nm.dsp"
efield_file = r"2D Differential Mode Ex.dsp"
elements, nodes, elem_solution, bounding_box = import_data.load_dsp(os.path.join(path, potential_file))
xdata, ydata, Udata = interpolate_slow.prepare_for_interpolation(elements, nodes, elem_solution)
x0 = -2.0 # Starting point for y
k = 151 # This defines the sampling
xeval = anneal.construct_symmetric_y(x0, k)
xinterp, yinterp, Uinterp = interpolate_slow.evaluate_on_grid(xdata, ydata, Udata, xeval=xeval,
yeval=helium_height * 1E6, clim=(0.00, 1.00), plot_axes='xy',
linestyle='None',
cmap=plt.cm.viridis, plot_data=False,
**common.plot_opt("darkorange", msize=6))
# Mirror around the y-axis
xsize = len(Uinterp[0])
Uinterp_symmetric = np.zeros(2 * xsize)
Uinterp_symmetric[:xsize] = Uinterp[0]
Uinterp_symmetric[xsize:] = Uinterp[0][::-1]
x_symmetric = np.zeros(2 * xsize)
x_symmetric[:xsize] = xinterp[0]
x_symmetric[xsize:] = -xinterp[0][::-1]
elements, nodes, elem_solution, bounding_box = import_data.load_dsp(os.path.join(path, efield_file))
xdata, ydata, Udata = interpolate_slow.prepare_for_interpolation(elements, nodes, elem_solution)
k = 251 # This defines the sampling
xeval = anneal.construct_symmetric_y(x0, k)
fig0 = plt.figure(figsize=(5., 3.))
xinterp, yinterp, Uinterp = interpolate_slow.evaluate_on_grid(
xdata,
ydata,
Udata,
xeval=xeval,
yeval=np.linspace(-1, 1, 151),
clim=(0.00, 1.00),
plot_axes='xy',
linestyle='None',
cmap=plt.cm.viridis,
plot_data=True,
**common.plot_opt("darkorange", msize=6))
xinterp, yinterp, Uinterp = interpolate_slow.evaluate_on_grid(
xdata,
ydata,
Udata,
xeval=xeval,
yeval=h,
clim=(0.00, 1.00),
plot_axes='xy',
linestyle='None',
cmap=plt.cm.viridis,
plot_data=False,
**common.plot_opt("darkorange", msize=6))
# Mirror around the y-axis
x0 = -2.0 # Starting point for y
k = 151 # This defines the sampling
xeval = anneal.construct_symmetric_y(x0, k)
xinterp, yinterp, Uinterp = interpolate_slow.evaluate_on_grid(
xdata,
ydata,
Udata,
xeval=xeval,
yeval=h,
clim=(0.00, 1.00),
plot_axes='xy',
linestyle='None',
cmap=plt.cm.viridis,
plot_data=False,
**common.plot_opt("darkorange", msize=6))
# Mirror around the y-axis
xsize = len(Uinterp[0])
Uinterp_symmetric = np.zeros(2 * xsize)
Uinterp_symmetric[:xsize] = Uinterp[0]
Uinterp_symmetric[xsize:] = Uinterp[0][::-1]
x_symmetric = np.zeros(2 * xsize)
x_symmetric[:xsize] = xinterp[0]
x_symmetric[xsize:] = -xinterp[0][::-1]
elements, nodes, elem_solution, bounding_box = import_data.load_dsp(
os.path.join(path, efield_file))
xdata, ydata, Udata = interpolate_slow.prepare_for_interpolation(