n = 1
atoms = []
for i = -n:n
    for j = -n:n
        for k = -n:n
            push!(atoms, (a1 .+ [i j k]));
            push!(atoms, (a2 .+ [i j k]));
            push!(atoms, (a3 .+ [i j k]));
            push!(atoms, (a4 .+ [i j k]));
            push!(atoms, (a5 .+ [i j k]));
            push!(atoms, (a6 .+ [i j k]));
            push!(atoms, (a7 .+ [i j k]));
            push!(atoms, (a8 .+ [i j k]));
            push!(atoms, (a9 .+ [i j k]));
            push!(atoms, (a10 .+ [i j k]));
            push!(atoms, (a11 .+ [i j k]));
            push!(atoms, (a12 .+ [i j k]));
            push!(atoms, (a13 .+ [i j k]));
            push!(atoms, (a14 .+ [i j k]));
            push!(atoms, (a15 .+ [i j k]));
            push!(atoms, (a16 .+ [i j k]));
            push!(atoms, (a17 .+ [i j k]));
            push!(atoms, (a18 .+ [i j k]));
            push!(atoms, (a19 .+ [i j k]));
            push!(atoms, (a20 .+ [i j k]));
            push!(atoms, (a21 .+ [i j k]));
            push!(atoms, (a22 .+ [i j k]));
            push!(atoms, (a23 .+ [i j k]));
            push!(atoms, (a24 .+ [i j k]));
            push!(atoms, (a25 .+ [i j k]));
            push!(atoms, (a26 .+ [i j k]));
            push!(atoms, (a27 .+ [i j k]));
            push!(atoms, (a28 .+ [i j k]));
            push!(atoms, (a29 .+ [i j k]));
            push!(atoms, (a30 .+ [i j k]));
            push!(atoms, (a31 .+ [i j k]));
            push!(atoms, (a32 .+ [i j k]));
            push!(atoms, (a33 .+ [i j k]));
            push!(atoms, (a34 .+ [i j k]));
            push!(atoms, (a35 .+ [i j k]));
            push!(atoms, (a36 .+ [i j k]));
        end
    end
end
atoms = vcat(atoms...)
unit_cell = similar(atoms); mul!(unit_cell, atoms, M1);
unit_cell = round.(unit_cell, digits = 5)
unit_cell = unique(unit_cell, dims = 1)
unit_cell[.!((-a1 .< unit_cell[:, 1] .< a1) .| (-b1 .< unit_cell[:, 2] .< b1) .| (-c1 .< unit_cell[:, 3] .< c1)), :]

M = [a b*cosd(gamma) c*cosd(beta); 0 b*sind(gamma) (c/sind(gamma))*(cosd(alpha)-cosd(beta)*cosd(gamma)); 0 0 (c*((1-cosd(alpha)^2 - cosd(beta)^2 - cosd(gamma)^2 + (2)*cosd(alpha)* cosd(beta)* cosd(gamma))^0.5))/(sind(gamma))]
M1 = M .* abc
T1 = [0 0 0]
T2 = [0 0 0.5]
T3 = ([0 0 0] + [2/3 1/3 1/3]) - floor.([0 0 0] + [2/3 1/3 1/3])
T4 = ([0 0 0.5] + [2/3 1/3 1/3]) - floor.([0 0 0.5] + [2/3 1/3 1/3])
T5 = ([0 0 0] + [1/3 2/3 2/3]) - floor.([0 0 0] + [1/3 2/3 2/3])
T6 = ([0 0 0.5] + [1/3 2/3 2/3]) - floor.([0 0 0.5] + [1/3 2/3 2/3])
a1 = [x y z] + T1
a2 = [-y x - y z] + T1
a3 = [-x + y -x z] + T1
a4 = [y x -z] + T2
a5 = [x - y -y -z] + T2
a6 = [-x -x + y  -z] + T2
a7 = [-x  -y  -z] + T1
a8 = [y  -x + y  -z] + T1
a9 = [x - y  x  -z] + T1
a10 = [-y  -x  z] + T2
a11 = [-x + y  y  z] + T2
a12 = [x  x - y  z] + T2
a13 = [x  y  z] + T3
a14 = [-y  x - y  z] + T3
a15 = [-x + y  -x  z] + T3
a16 = [y  x  -z] + T4
a17 = [x - y  -y  -z] + T4
a18 = [-x  -x + y  -z] + T4
a19 = [-x  -y  -z] + T3
a20 = [y  -x + y  -z] + T3
a21 = [x - y  x  -z] + T3
a22 = [-y  -x  z] + T4
a23 = [-x + y  y  z] + T4
a24 = [x  x - y  z] + T4
a25 = [x  y  z] + T5
a26 = [-y  x - y  z] + T5
a27 = [-x + y  -x  z] + T5
a28 = [y  x  -z] + T6
a29 = [x - y  -y  -z] + T6
a30 = [-x  -x + y  -z] + T6
a31 = [-x  -y  -z] + T5
a32 = [y  -x + y  -z] + T5
a33 = [x - y  x  -z] + T5
a34 = [-y  -x  z] + T6
a35 = [-x + y  y  z] + T6
a36 = [x  x - y  z] + T6