diff options
Diffstat (limited to 'buch/papers/laguerre/scripts/rel_error_shifted.py')
-rw-r--r-- | buch/papers/laguerre/scripts/rel_error_shifted.py | 13 |
1 files changed, 11 insertions, 2 deletions
diff --git a/buch/papers/laguerre/scripts/rel_error_shifted.py b/buch/papers/laguerre/scripts/rel_error_shifted.py index 1515c6e..dc9d177 100644 --- a/buch/papers/laguerre/scripts/rel_error_shifted.py +++ b/buch/papers/laguerre/scripts/rel_error_shifted.py @@ -1,10 +1,18 @@ if __name__ == "__main__": + import matplotlib as mpl import matplotlib.pyplot as plt import numpy as np import scipy.special import gamma_approx as ga + mpl.rcParams.update( + { + "mathtext.fontset": "stix", + "font.family": "serif", + "font.serif": "TeX Gyre Termes", + } + ) n = 8 # order of Laguerre polynomial N = 200 # number of points in interval @@ -19,7 +27,7 @@ if __name__ == "__main__": ax.semilogy(x, rel_error, label=f"$m={target}$", linewidth=3) gamma_lgo = ga.eval_laguerre_gamma(x, n=n, func="optimal_shifted") rel_error = np.abs(ga.calc_rel_error(gamma, gamma_lgo)) - ax.semilogy(x, rel_error, "m", linestyle="dotted", label="$m^*$", linewidth=3) + ax.semilogy(x, rel_error, "m", linestyle=":", label="$m^*$", linewidth=3) ax.set_xlim(x[0], x[-1]) ax.set_ylim(5e-9, 5e-8) ax.set_xlabel(r"$z$") @@ -27,5 +35,6 @@ if __name__ == "__main__": ax.set_xticks(np.linspace(0, 1, 11), minor=True) ax.grid(1, "both") ax.legend(ncol=1, fontsize=ga.fontsize) - fig.savefig(f"{ga.img_path}/rel_error_shifted.pgf") + # fig.savefig(f"{ga.img_path}/rel_error_shifted.pgf") + fig.savefig(f"{ga.img_path}/rel_error_shifted.pdf") # plt.show() |