gem.utility.delta_fit

Functions

F_of_H(H, beta)

build_H(Lambda, Lambda_c, D, R)

dF_spectral(H, dH, beta[, tol])

dH_dD(dD, R)

dD is the perturbation to D.

dH_dLambda(dLambda, n)

dH_dLambdac(dLambdac, n)

dH_dR(dR, D)

dR is the perturbation to R.

fermi(eps, beta)

fermi_prime(eps, beta)

jacobian_LR(x, beta, Lambda_c, D, ...)

Jacobian for Lambda+R fit

jacobian_LR_movement(x, beta, Lambda_c, D, ...)

Stack original Jacobian with movement-penalty rows: factor * I.

jacobian_LcD(x, beta, Lambda, R, F11_target, ...)

Jacobian for Lambda_c+D

jacobian_LcD_movement(x, beta, Lambda, R, ...)

pack_params(Lambda, R)

residual_LR(x, beta, Lambda_c, D, ...)

Residual for Lambda+R fit.

residual_LR_movement(x, beta, Lambda_c, D, ...)

Original residuals + movement-penalty residuals sqrt(2*alpha)*(x-x0).

residual_LcD(x, beta, Lambda, R, F11_target, ...)

Residual for Lambda_c+D fit

residual_LcD_movement(x, beta, Lambda, R, ...)

solve_F_dF_LR_with_movement(beta, Lambda_c, ...)

Least-squares solve with quadratic movement penalty around starting x0.

solve_F_dF_LcD_with_movement(beta, Lambda, ...)

unpack_params(x, n, p)

update_hybridization_thermal_penalty(...[, ...])

Function to update the hybridization function parameters D and Lambda_c.

update_self_energy_thermal_penalty(Lambda0, ...)

Replacement of new_self_energy that includes movement penalty alpha * ||x-x0||^2.