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