yapanuwan
/
DVR-jl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

2 bug fixes

This commit is contained in:
ysyapa 2023-08-18 19:30:53 +00:00
parent f74f206d8b
commit 1e768f4496
2 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
common.jl
View File

@ -90,7 +90,7 @@ function calculate_Vs(s::system{T}, V_twobody::Function, ϕ::T, n_image::Int)::A
end end
end end
end end
Vs[i] *= s.multiplicity[i[1]] Vs[i] *= sqrt(s.multiplicity[i[1]])
end end
return Vs return Vs
end end

14
irrep.jl
View File

@ -47,13 +47,13 @@ function sym_reduce(s, K_partial)
K_partial_y = kron(kron(I, K_partial), I) K_partial_y = kron(kron(I, K_partial), I)
K_partial_z = kron(kron(I, I), K_partial) K_partial_z = kron(kron(I, I), K_partial)
for i in s.unique_i for (i, j) in enumerate(s.unique_i)
K_partial_x[i, :] *= s.multiplicity[i] K_partial_x[j, :] *= sqrt(s.multiplicity[i])
K_partial_x[:, i] *= s.multiplicity[i] K_partial_x[:, j] *= sqrt(s.multiplicity[i])
K_partial_y[i, :] *= s.multiplicity[i] K_partial_y[j, :] *= sqrt(s.multiplicity[i])
K_partial_y[:, i] *= s.multiplicity[i] K_partial_y[:, j] *= sqrt(s.multiplicity[i])
K_partial_z[i, :] *= s.multiplicity[i] K_partial_z[j, :] *= sqrt(s.multiplicity[i])
K_partial_z[:, i] *= s.multiplicity[i] K_partial_z[:, j] *= sqrt(s.multiplicity[i])
end end
K_partial_x = K_partial_x[s.unique_i, s.unique_i] K_partial_x = K_partial_x[s.unique_i, s.unique_i]