Sqrt weights

This commit is contained in:
ysyapa 2023-08-18 20:57:54 +00:00
parent 443aaae391
commit 4063079f59
1 changed files with 6 additions and 6 deletions

View File

@ -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]