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

Orthogonality check

This commit is contained in:
Nuwan Yapa 2024-07-01 12:27:22 -04:00
parent 8431b4b7e3
commit ccb1b539c7
1 changed files with 8 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
test/3body_Berggren_orthogonality.jl Normal file
View File

@ -0,0 +1,8 @@
include("../berggren_3body_resonance.jl")
@time "Eigenvectors" evals, evecs = eigs(H, sigma=target, maxiter=5000, tol=1e-5, ritzvec=true, check=1)
weights_mat = spdiagm(repeat(kron(ws, ws), jmax + 1))
N = transpose(evecs) * weights_mat * evecs
display(abs.(N))