for i in range(self.number_of_bands):
                occupancy[i] = vasp_data_occupancies.reshape(
                    number_of_kpoints * 2, self.number_of_bands
                ).T[i][:number_of_kpoints]
                occupancy[self.number_of_bands + i] = vasp_data_occupancies.reshape(