From c08e4690832f84ac2ea0e3fbc29bbcaaadff452f Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Mon, 4 Mar 2024 15:13:42 -0500 Subject: [PATCH] =?UTF-8?q?Complex-=CF=89=20idea=20testing?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- ho_basis_benchmark.jl | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/ho_basis_benchmark.jl b/ho_basis_benchmark.jl index c3d5c5e..38d5cab 100644 --- a/ho_basis_benchmark.jl +++ b/ho_basis_benchmark.jl @@ -2,7 +2,7 @@ using LinearAlgebra, Plots include("ho_basis.jl") include("p_space.jl") -ω = 0.5 +ω = -0.5im μ = 0.5 l = 0 V1 = -5 @@ -17,20 +17,20 @@ ls = fill(l, n_max + 1) T = sp_T_matrix(ns, ls; ω=ω, μ=μ) V = V1 .* V_Gaussian.(R1, l, ns, transpose(ns); ω=ω) + V2 .* V_Gaussian.(R2, l, ns, transpose(ns); ω=ω) -cs = range(1.35, 0.9, 5) +cs = range(1.25, 0.25, 10) -E = similar(cs) -bench_E = similar(cs) +E = similar(cs, ComplexF64) +bench_E = similar(cs, ComplexF64) for (j, c) in enumerate(cs) H = T + c .* V evals = eigvals(H) - bench_E[j] = quick_pole_E((p, q) -> c*(V1*g0(R1, p, q) + V2*g0(R2, p, q)), μ; cs_angle=0) + bench_E[j] = quick_pole_E((p, q) -> c*(V1*g0(R1, p, q) + V2*g0(R2, p, q)), μ; cs_angle=0.4, meshpoints=512) i = argmin(abs.(evals .- bench_E[j])) E[j] = evals[i] end scatter(real.(bench_E), imag.(bench_E), label="p-space") scatter!(real.(E), imag.(E), label="HO basis") -xlims!(-0.7,0) -ylims!(-0.1,0.1) +xlims!(-0.7,0.4) +ylims!(-0.25,0.25)