Nick Papior

See https://github.com/nschloe/quadpy/issues/250 where I posted a working solution (to be simplified and beautified :)

Could you make an example, I don't fully understand what you mean?

But this is just `MonkhorstPack(..., [x, y, z], trs=False)`, then you can do: ```python eig = bz.apply.array.eigh() plot(bz.k[:, 0], bz.k[:, 1], eig) ``` ? I don't think it deserves another...

No, that is not correct, the `size` argument define the size of the BZ, and `displacement` defines the center of the BZ. You can even do *better k-points for SCF*,...

As for > Planes that are orthogonal to a cartesian coordinate I don't exactly know what you mean, k-points are always in terms of the lattice vector. ?

So what you say is that you want to equi-sample inside the 3D triangle that spans Gamma -> M -> X -> R? Do you envision this as a surface...

> For 2D band structures the k sampling would need to be in a surface, because then you can represent them as: `xy` the coordinates of the surface and `z`...

I think it is a really good idea, however it will introduce some corner case problems: - `Geometry` this class should be trivially handled - `SparseGeometry` (Hamiltonians etc) how should...

Partially, since lots of code is relying on Geometry close calls. I am also worried about the documentation issue, i.e. which routines uses *atoms with orbitals* and which are using...

I agree, this is why I generally recommend one to create a `construct` function and pass that, it is a bit easier, but still *loong*. The problem of the `for...