triqs_cthyb.tail_fit.tail_fit

triqs_cthyb.tail_fit.tail_fit(Sigma_iw, fit_min_n=None, fit_max_n=None, fit_min_w=None, fit_max_w=None, fit_max_moment=None, fit_known_moments=None)[source]

Fit a high frequency 1/(iw)^n expansion of Sigma_iw and replace the high frequency part with the fitted high frequency expansion.

Either give frequency window to fit on in terms of matsubara frequencies index (fit_min_n/fit_max_n) or value (fit_min_w/fit_max_w).

Parameters:

Sigma_iw : Gf

Self-energy.

fit_min_n : int, optional, default=int(0.8*len(Sigma_iw.mesh))

Matsubara frequency index from which tail fitting should start.

fit_max_n : int, optional, default=int(len(Sigma_iw.mesh))

Matsubara frequency index at which tail fitting should end.

fit_min_w : float, optional

Matsubara frequency from which tail fitting should start.

fit_max_w : float, optional

Matsubara frequency at which tail fitting should end.

fit_max_moment : int, optional

Highest moment to fit in the tail of Sigma_iw.

fit_known_moments : ndarray.shape[order, Sigma_iw[0].target_shape], optional, default = None

Known moments of Sigma_iw, given as an numpy ndarray

Returns:

Sigma_iw : Gf

Self-energy.