20 lines
296 B
Julia
20 lines
296 B
Julia
include("../Hamiltonian.jl")
|
|
mode = cpu_tensor
|
|
T = Float32
|
|
|
|
V_gauss(r2) = -2 * exp(-r2 / 4)
|
|
|
|
d = 3
|
|
n = 3
|
|
N = 20
|
|
L = 15
|
|
n_imag = 1
|
|
ϕ = 0
|
|
|
|
s = system{T}(d, n, N, L)
|
|
H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
|
|
@time evals, _, info = eig(H, 5)
|
|
|
|
print(info.numops, " operations")
|
|
display(evals)
|