22 lines
431 B
Julia
22 lines
431 B
Julia
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))
|