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

Major bug fix

This commit is contained in:
ysyapa 2024-02-09 19:49:45 +00:00
parent 37e68b78cc
commit 01767fd713
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
EC_test.ipynb
View File

@ -97,7 +97,8 @@
"end\n", "end\n",
"\n", "\n",
"EC_basis = hcat(EC_basis, conj.(EC_basis)) # CA-EC\n", "EC_basis = hcat(EC_basis, conj.(EC_basis)) # CA-EC\n",
"N_EC = transpose(EC_basis) * EC_basis\n", "EC_basis_w = EC_basis .* w\n",
"N_EC = transpose(EC_basis_w) * EC_basis\n",
"\n", "\n",
"extrapolate_points = range(0.75, 0.40, 8)\n", "extrapolate_points = range(0.75, 0.40, 8)\n",
"\n", "\n",
@ -107,9 +108,8 @@
"for (j, c) in enumerate(extrapolate_points)\n", "for (j, c) in enumerate(extrapolate_points)\n",
" exact_E[j] = quick_pole_E(V_system(c))\n", " exact_E[j] = quick_pole_E(V_system(c))\n",
"\n", "\n",
" EC_basis_w = EC_basis .* w\n",
" H = get_H_matrix(V_system(c), p, w)\n", " H = get_H_matrix(V_system(c), p, w)\n",
" H_EC = transpose(EC_basis) * H * EC_basis\n", " H_EC = transpose(EC_basis_w) * H * EC_basis\n",
" evals = eigvals(H_EC, N_EC)\n", " evals = eigvals(H_EC, N_EC)\n",
" i = argmin(abs.(evals .- exact_E[j]))\n", " i = argmin(abs.(evals .- exact_E[j]))\n",
" extrapolate_E[j] = evals[i]\n", " extrapolate_E[j] = evals[i]\n",