From 3792647cbbb025b5e13181d04866bb8bc2532a54 Mon Sep 17 00:00:00 2001 From: ysyapa Date: Fri, 18 Aug 2023 19:38:55 +0000 Subject: [PATCH] Revert "2 bug fixes" --- common.jl | 2 +- irrep.jl | 14 +++++++------- 2 files changed, 8 insertions(+), 8 deletions(-) diff --git a/common.jl b/common.jl index 030103e..672551f 100644 --- a/common.jl +++ b/common.jl @@ -90,7 +90,7 @@ function calculate_Vs(s::system{T}, V_twobody::Function, ϕ::T, n_image::Int)::A end end end - Vs[i] *= sqrt(s.multiplicity[i[1]]) + Vs[i] *= s.multiplicity[i[1]] end return Vs end diff --git a/irrep.jl b/irrep.jl index 1027c3f..2df52aa 100644 --- a/irrep.jl +++ b/irrep.jl @@ -47,13 +47,13 @@ function sym_reduce(s, K_partial) K_partial_y = kron(kron(I, K_partial), I) K_partial_z = kron(kron(I, I), K_partial) - for (i, j) in enumerate(s.unique_i) - K_partial_x[j, :] *= sqrt(s.multiplicity[i]) - K_partial_x[:, j] *= sqrt(s.multiplicity[i]) - K_partial_y[j, :] *= sqrt(s.multiplicity[i]) - K_partial_y[:, j] *= sqrt(s.multiplicity[i]) - K_partial_z[j, :] *= sqrt(s.multiplicity[i]) - K_partial_z[:, j] *= sqrt(s.multiplicity[i]) + for i in s.unique_i + K_partial_x[i, :] *= s.multiplicity[i] + K_partial_x[:, i] *= s.multiplicity[i] + K_partial_y[i, :] *= s.multiplicity[i] + K_partial_y[:, i] *= s.multiplicity[i] + K_partial_z[i, :] *= s.multiplicity[i] + K_partial_z[:, i] *= s.multiplicity[i] end K_partial_x = K_partial_x[s.unique_i, s.unique_i]