using Plots
include("../common.jl")
include("../p_space.jl")

vertices = [0, 1 - 0.5im, 2, 6]
subdivisions = [64, 64, 64]
p, w = get_mesh(vertices, subdivisions)

scatter(real.(p), imag.(p))

# ComplexScaling-FV: Eq. (54)
V_pq(p, q) = -10 * g1(1, p, q)

H = get_H_matrix(V_pq, p, w)
evals = eigen(H).values

E_target = 0.258 - 0.164im
E = nearest(evals, E_target)

print("E = $E")
scatter(real.(evals), imag.(evals), xlim = (0,1))