From 4063079f59841f92f6d5da329ed341f1d8e224a5 Mon Sep 17 00:00:00 2001 From: ysyapa Date: Fri, 18 Aug 2023 20:57:54 +0000 Subject: [PATCH] Sqrt weights --- irrep.jl | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/irrep.jl b/irrep.jl index 2df52aa..cb2bd8c 100644 --- a/irrep.jl +++ b/irrep.jl @@ -48,12 +48,12 @@ function sym_reduce(s, K_partial) 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] + K_partial_x[i, :] *= sqrt(s.multiplicity[i]) + K_partial_x[:, i] *= sqrt(s.multiplicity[i]) + K_partial_y[i, :] *= sqrt(s.multiplicity[i]) + K_partial_y[:, i] *= sqrt(s.multiplicity[i]) + K_partial_z[i, :] *= sqrt(s.multiplicity[i]) + K_partial_z[:, i] *= sqrt(s.multiplicity[i]) end K_partial_x = K_partial_x[s.unique_i, s.unique_i]