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

Minor optimization

This commit is contained in:
ysyapa 2024-02-09 17:21:04 +00:00
parent 9696d27e2e
commit 37e68b78cc
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
EC_test.ipynb
View File

@ -97,6 +97,7 @@
"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",
"\n", "\n",
"extrapolate_points = range(0.75, 0.40, 8)\n", "extrapolate_points = range(0.75, 0.40, 8)\n",
"\n", "\n",
@ -109,7 +110,6 @@
" EC_basis_w = EC_basis .* w\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) * H * EC_basis\n",
" N_EC = transpose(EC_basis) * 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",