From 56a88089382c0856d8f21279d9229d4bb80175fb Mon Sep 17 00:00:00 2001 From: ysyapa Date: Wed, 16 Aug 2023 22:06:19 +0000 Subject: [PATCH] Small bug fixed --- irrep.jl | 40 ++++++++++++++++++++++++---------------- 1 file changed, 24 insertions(+), 16 deletions(-) diff --git a/irrep.jl b/irrep.jl index b4ddfae..249d19f 100644 --- a/irrep.jl +++ b/irrep.jl @@ -2,24 +2,32 @@ using DelimitedFiles rotations = readdlm("rotations.mat", ',', Int, '\n') rotations = reshape(rotations, (24, 3, 3)) -N = 24 -ks = -N÷2:N÷2-1 -lattice = hcat((collect.(Iterators.product(ks, ks, ks)))...) -labels = reshape(collect(1:N^3), (N, N, N)) - -for r in 1:24 - rotated_lattice = Matrix(rotations[r, :, :]) * lattice - for index in 1:N^3 - rotated_lattice_point = rotated_lattice[:, index] - (i, j, k) = mod1.(rotated_lattice_point .+ (N÷2 + 1), N) - old_label = max(labels[index], labels[i, j, k]) - new_label = min(labels[index], labels[i, j, k]) - if old_label != new_label - for o in findall(isequal(old_label), labels) - labels[o] = new_label +function get_A1_labels(N::Int) + rotations = readdlm("rotations.csv", ',', Int, '\n') + rotations = reshape(rotations, (24, 3, 3)) + + ks = -N÷2:N÷2-1 + lattice = hcat((collect.(Iterators.product(ks, ks, ks)))...) + labels = reshape(collect(1:N^3), (N, N, N)) + + for r in 1:24 + rotated_lattice = Matrix(rotations[r, :, :]) * lattice + for index in 1:N^3 + rotated_lattice_point = rotated_lattice[:, index] + (i, j, k) = mod1.(rotated_lattice_point .+ (N÷2 + 1), N) + old_label = max(labels[index], labels[i, j, k]) + new_label = min(labels[index], labels[i, j, k]) + if old_label != new_label + for o in findall(isequal(old_label), labels) + labels[o] = new_label + end end end end + + return labels end -unique(labels) \ No newline at end of file +function sym_reduce(N::Int, K_full) + +end \ No newline at end of file