diff --git a/p_space.jl b/p_space.jl index 562b616..8363fee 100644 --- a/p_space.jl +++ b/p_space.jl @@ -45,7 +45,7 @@ function identify_pole_i(p, evals, μ=0.5) return current_i 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) evals = eigvals(get_H_matrix(V_pq, p, w, μ)) return evals[identify_pole_i(p, evals, μ)]