-
Notifications
You must be signed in to change notification settings - Fork 0
/
fig_any-compare-pareto.py
50 lines (40 loc) · 1.58 KB
/
fig_any-compare-pareto.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
import matplotlib.pyplot as plt
import numpy as np
from matplotlib.backends.backend_pdf import PdfPages
plt.rc("text", usetex=True)
fNs = np.array(
[
"data/poiss10000000.pareto100_1000_10000.seed1186.psd.csv",
"data/const1000000.pareto100_1000_10000.seed27082.psd.csv",
"data/unif0_10000000.pareto100_1000_10000.seed15769.psd.csv",
"data/pareto100_1000_10000.pareto100_1000_10000.seed20358.psd.csv",
]
)
with PdfPages("figs/any-compare-pareto.pdf") as pdfFile:
fig = plt.figure(figsize=(3, 2))
plt.loglog()
plt.ylim([3e-8, 3e4])
plt.xlim([3e-6, 3e0])
ax = plt.gca()
ax.grid(False)
# ax.grid(visible=True, linestyle="--", dashes=(15, 15), c="#cccccc", lw=0.3)
ax.get_yaxis().set_tick_params(direction="out")
ax.get_xaxis().set_tick_params(direction="out")
ax.set_xlabel(r"$f$")
ax.set_ylabel(r"$S(f)$")
ax.set_xticks([1e-5, 1e-4, 1e-3, 1e-2, 1e-1, 1e0])
ax.set_yticks([1e-6, 1e-4, 1e-2, 1e0, 1e2, 1e4])
data = 10 ** np.loadtxt(fNs[0], delimiter=",")
data[:, 1] = data[:, 1]
plt.plot(data[:, 0], data[:, 1], "r-", lw=4)
data = 10 ** np.loadtxt(fNs[1], delimiter=",")
plt.plot(data[:, 0], data[:, 1], "g-", lw=4)
plt.plot(data[:, 0], 0.03 * data[:, 2], "k--")
data = 10 ** np.loadtxt(fNs[2], delimiter=",")
plt.plot(data[:, 0], data[:, 1], "b-", lw=4)
data = 10 ** np.loadtxt(fNs[3], delimiter=",")
plt.plot(data[:, 0], data[:, 1], "m-", lw=4)
plt.plot(data[:, 0], 10 * data[:, 2], "k--")
plt.tight_layout(w_pad=1, h_pad=1, pad=0.5)
pdfFile.savefig(fig)
plt.close()