Should paths for plane group 2 include parts of the BZ boundary? about brillouin.jl HOT 1 OPEN

function prep_bz_path(Rs)
    Gs = reciprocalbasis(Rs)
    cell = wignerseitz(Gs)
    kp = irrfbz_path(2, Rs)

    a, b = norm(Rs[1]), norm(Rs[2])
    cosβ = dot(Rs[1], Rs[2])/(a*b)
    sinβ = sin(acos(cosβ))
    Y = ((((1 + (b / a) * cosβ) / 2) / sinβ) / sinβ)
    N = (1 / 2 + (Y * a * cosβ) / b)
    pts = Dict(
        :Γ => [0, 0], 
        :B => [1/2, 0], 
        :Y => [0, 1/2,], 
        :A => [1/2, 1/2],
        # not included in path:
        :H => [1-N, -Y],
        :H₂ => [N, -1+Y],
        :H₄ => [-N, -Y],
        :H′ => [-1+N, Y],
        :H₂′ => [-N, 1-Y],
        :H₄′ => [N, Y],
        )

    merge!(kp.points, pts)
    return cell, kp
end

b = 1.0
β = π/2 + π/3
Rs = crystal(1.0, b, β)
cell, kp = prep_bz_path(Rs)

f = Figure();
ax = Axis(f[1,1]; aspect=1.0)

plot!(ax, cell)
plot!(ax, kp)
f

sl_β = Slider(f[1, 2], range = range(0.0+π/10, π-π/10, 100), startvalue = β, horizontal = false)
on(sl_β.value) do β
    Rs = crystal(1.0, b, β)
    cell, kp = prep_bz_path(Rs)
    empty!(ax)
    plot!(ax, cell)
    plot!(ax, kp)
end
f