From 1e768f44960a69d3848b0eda5965d9ca0526b063 Mon Sep 17 00:00:00 2001 From: ysyapa Date: Fri, 18 Aug 2023 19:30:53 +0000 Subject: [PATCH] 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 672551f..030103e 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] *= s.multiplicity[i[1]] + Vs[i] *= sqrt(s.multiplicity[i[1]]) end return Vs end diff --git a/irrep.jl b/irrep.jl index 2df52aa..1027c3f 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 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] + 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]) end K_partial_x = K_partial_x[s.unique_i, s.unique_i]