Small bug fix

This commit is contained in:
ysyapa 2024-02-09 16:40:14 +00:00
parent 11ce5304bb
commit 34d33629a2
1 changed files with 2 additions and 2 deletions

View File

@ -47,6 +47,6 @@ end
function quick_pole_E(V_pq, μ=0.5; cs_angle=0.5, cutoff=8.0, meshpoints=256)
p, w = get_mesh([0, cutoff * exp(-1im * cs_angle)], meshpoints)
evals = eigvals(get_H_matrix(V_pq, p, w))
return evals[identify_pole_i(p, evals)]
evals = eigvals(get_H_matrix(V_pq, p, w, μ))
return evals[identify_pole_i(p, evals, μ)]
end