1
Fork 0

Compare commits

..

No commits in common. "a000c67e93edfffa9ca1004624fad1ef04012712" and "045b6ae439fc877ea0acadfa19808a45a030defe" have entirely different histories.

View file

@ -81,21 +81,18 @@ ax3d.quiver3D(
ax3d.set(xlabel="x (cm)", ylabel="y (cm)", zlabel="z (cm)") ax3d.set(xlabel="x (cm)", ylabel="y (cm)", zlabel="z (cm)")
theta = np.angle(raw_ts["x"] + 1j * raw_ts["y"]).mean() theta = np.angle(raw_ts["x"] + 1j * raw_ts["y"]).mean()
fig2dv, ax2dv = plt.subplots(figsize=(5/2.54, 2/3*10/2.54), dpi=200, constrained_layout=True) fig2dv, ax2dv = plt.subplots()
x0 = ts_flt["x"] * np.cos(theta) + ts_flt["y"] * np.sin(theta) x0 = ts_flt["x"] * np.cos(theta) + ts_flt["y"] * np.sin(theta)
#ax2dv.plot(x0, z0, c="#0066ff", lw=1) ax2dv.plot(x0, z0, c="#0066ff", lw=1)
ax2dv.quiver( ax2dv.quiver(
x0[:-1], x0[:-1],
z0[:-1], z0[:-1],
np.diff(x0)[:], np.diff(x0)[:],
np.diff(z0)[:], np.diff(z0)[:],
color="k", color="#0066ff",
scale_units="xy",
scale=1,
) )
ax2dv.grid(c="k", alpha=.2) ax2dv.grid()
ax2dv.set(aspect="equal", xlabel="x (cm)", ylabel="z (cm)") ax2dv.set(aspect="equal")
fig2dv.savefig("out_orbitals.pdf") fig2dv.savefig("out_orbitals.pdf")
fig2dv.savefig("out_orbitals.jpg")
plt.show() plt.show()