lattice/pyerrors
1
0
Fork 0
mirror of https://github.com/fjosw/pyerrors.git synced 2025-11-29 20:46:54 +01:00
Code Issues Projects Releases Packages Wiki Activity

small switch of lines

Browse source
This commit is contained in:
Justus Kuhlmann 2025-11-03 16:04:59 +00:00
commit 1e00caf3c1
1 changed files with 4 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

7
pyerrors/input/misc.py
View file

@ -60,7 +60,10 @@ def fit_t0(t2E_dict: dict[float, Obs], fit_range: int, plot_fit: Optional[bool]=
fit_result = fit_lin(x, y) fit_result = fit_lin(x, y)
retval = (-fit_result[0] / fit_result[1])
if plot_fit is True: if plot_fit is True:
retval.gamma_method()
plt.figure() plt.figure()
gs = gridspec.GridSpec(2, 1, height_ratios=[3, 1], wspace=0.0, hspace=0.0) gs = gridspec.GridSpec(2, 1, height_ratios=[3, 1], wspace=0.0, hspace=0.0)
ax0 = plt.subplot(gs[0]) ax0 = plt.subplot(gs[0])
@ -72,8 +75,6 @@ def fit_t0(t2E_dict: dict[float, Obs], fit_range: int, plot_fit: Optional[bool]=
yplot = [fit_result[0] + fit_result[1] * xi for xi in xplot] yplot = [fit_result[0] + fit_result[1] * xi for xi in xplot]
[yi.gamma_method() for yi in yplot] [yi.gamma_method() for yi in yplot]
ax0.fill_between(xplot, y1=[yi.value - yi.dvalue for yi in yplot], y2=[yi.value + yi.dvalue for yi in yplot]) ax0.fill_between(xplot, y1=[yi.value - yi.dvalue for yi in yplot], y2=[yi.value + yi.dvalue for yi in yplot])
retval = (-fit_result[0] / fit_result[1])
retval.gamma_method()
ylim = ax0.get_ylim() ylim = ax0.get_ylim()
ax0.fill_betweenx(ylim, x1=retval.value - retval.dvalue, x2=retval.value + retval.dvalue, color='gray', alpha=0.4) ax0.fill_betweenx(ylim, x1=retval.value - retval.dvalue, x2=retval.value + retval.dvalue, color='gray', alpha=0.4)
ax0.set_ylim(ylim) ax0.set_ylim(ylim)
@ -96,7 +97,7 @@ def fit_t0(t2E_dict: dict[float, Obs], fit_range: int, plot_fit: Optional[bool]=
ax1.set_xlabel(r'$t/a^2$') ax1.set_xlabel(r'$t/a^2$')
plt.draw() plt.draw()
return -fit_result[0] / fit_result[1] return retval
def read_pbp(path: str, prefix: str, **kwargs): def read_pbp(path: str, prefix: str, **kwargs):