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]