yapanuwan
/
DVR-jl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Increase N

This commit is contained in:
ysyapa 2023-08-09 18:22:05 +00:00
parent e6f32295f5
commit 490e94262b
3 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

0
ComplexScaling-FV.dat Normal file
View File

2
ComplexScaling-FV.jl
View File

@ -7,7 +7,7 @@ V_gauss(r2) =
d = 3 d = 3
n = 2 n = 2
N = 16 N = 64
ϕ = pi/6 ϕ = pi/6
open("ComplexScaling-FV.dat", "w") do f open("ComplexScaling-FV.dat", "w") do f

2
Hamiltonian.jl
View File

@ -134,7 +134,7 @@ function eig(H::Hamiltonian{T}, levels::Int; resonances = !H.hermitian)::Tuple{V
x₀ = CUDA.rand(Complex{T}, vectorDims(H)...) x₀ = CUDA.rand(Complex{T}, vectorDims(H)...)
synchronize() synchronize()
end end
evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = H.hermitian, tol = tolerance) evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = H.hermitian, tol = tolerance, krylovdim = levels * 2)
resonances || info.converged < levels && throw(error("Not enough convergence")) # don't check convergence for resonances resonances || info.converged < levels && throw(error("Not enough convergence")) # don't check convergence for resonances
if H.hermitian evals = real.(evals) end if H.hermitian evals = real.(evals) end
if H.mode == gpu_cutensor # to avoid possible GPU memory leak if H.mode == gpu_cutensor # to avoid possible GPU memory leak