diff --git a/irrep.jl b/irrep.jl index d1a4935..2df52aa 100644 --- a/irrep.jl +++ b/irrep.jl @@ -47,17 +47,6 @@ 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 s in 1:N^3 -# if labels[s] != s -# for mat in (K_partial_x, K_partial_y, K_partial_z) -# mat[labels[s], :] += mat[s, :] -# mat[s, :] = 0 -# mat[:, labels[s]] += mat[:, s] -# mat[:, s] = 0 -# end -# end -# end - for i in s.unique_i K_partial_x[i, :] *= s.multiplicity[i] K_partial_x[:, i] *= s.multiplicity[i]