Change default parameter

This commit is contained in:
ysyapa 2024-02-09 17:02:52 +00:00
parent 34d33629a2
commit 9696d27e2e
1 changed files with 1 additions and 1 deletions

View File

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