for i in range(self.number_of_bands):
                energies[i] = vasp_data_energies.reshape(
                    number_of_kpoints * 2,  # factor of 2 for each kpoint block
                    self.number_of_bands,
                ).T[i][:number_of_kpoints]