diff --git a/figures/sp-fd-2d-3d.png b/figures/sp-fd-2d-3d.png new file mode 100644 index 0000000..25235c7 Binary files /dev/null and b/figures/sp-fd-2d-3d.png differ diff --git a/figures/sp-fd-2d-3d.svg b/figures/sp-fd-2d-3d.svg new file mode 100644 index 0000000..ebc3ad0 --- /dev/null +++ b/figures/sp-fd-2d-3d.svg @@ -0,0 +1,1652 @@ + + + + + + + + 2023-03-19T10:35:10.109493 + image/svg+xml + + + Matplotlib v3.7.1, https://matplotlib.org/ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/graphs/sp-fd-2d-3d.py b/graphs/sp-fd-2d-3d.py new file mode 100644 index 0000000..a75bb77 --- /dev/null +++ b/graphs/sp-fd-2d-3d.py @@ -0,0 +1,50 @@ +from matplotlib import pyplot as plt + +from lib.lib import read_sp, read_xyz_alt, read_xy_alt, aggregate_sp_fd + +# %% +d2_raw = read_sp("../data/rust-sticking-probability", read_xy_alt) +d2 = aggregate_sp_fd(d2_raw) +d3_raw = read_sp("../data/rust-3d-offaxis-sp", read_xyz_alt) +d3 = aggregate_sp_fd(d3_raw) + +# %% +plt.plot( + d2.index, + d2.fd, + color='tab:blue', + label='D2' +) + +plt.fill_between( + d2.index, + d2.fd - d2.fd_std, + d2.fd + d2.fd_std, + alpha=0.2, + color='tab:blue', + label=f"Standard error band" +) + +plt.plot( + d3.index, + d3.fd, + color='tab:orange', + label='D3' +) + +plt.fill_between( + d3.index, + d3.fd - d3.fd_std, + d3.fd + d3.fd_std, + alpha=0.2, + color='tab:orange', + label=f"Standard error band" +) + +plt.xlabel("$p_{stick}$") +plt.ylabel("$fd$") +plt.legend() + +plt.savefig('../figures/sp-fd-2d-3d.svg') +plt.savefig('../figures/sp-fd-2d-3d.png') +plt.show()