Commit 35d5ec85 authored by Antonino D'Anna's avatar Antonino D'Anna

bug fix

parent c90e092c
...@@ -19,7 +19,7 @@ function __Q(nu,chi,nmc) ...@@ -19,7 +19,7 @@ function __Q(nu,chi,nmc)
z[j] = randn() z[j] = randn()
z2[j] = z[j]^2 z2[j] = z[j]^2
end end
aux = sum(z2[i]*eig[i] for i in in eachindex(eig)) aux = sum(z2[i]*eig[i] for i in eachindex(eig))
if aux < chi if aux < chi
Q+=1.0 Q+=1.0
end end
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment