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

Minor refactoring

This commit is contained in:
Nuwan Yapa 2025-01-06 15:59:53 -05:00
parent 9462cd43bf
commit 660b0b4715
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
helper.jl
View File

@ -59,7 +59,7 @@ function gram_schmidt(vecs, ws, precision=1e-10, verbose=false)
for j in findall(mask[1 : (i - 1)]) for j in findall(mask[1 : (i - 1)])
vecs[i] -= proj(j, i) vecs[i] -= proj(j, i)
end end
norm = sum(vecs[i] .* ws .* vecs[i]) norm = c_product(i, i)
verbose && println("Mag. of basis vector $i = $(abs(norm))") verbose && println("Mag. of basis vector $i = $(abs(norm))")
if abs(norm) > precision if abs(norm) > precision
vecs[i] = vecs[i] ./ sqrt(norm) vecs[i] = vecs[i] ./ sqrt(norm)