def test_perform_inversion(self, job_file, old_eigenvectors,
                               expected_eigenvectors):
        # Setup
        Job = InitJob(job_file)

        # Action
        new_eigenvectors = Job.Electron.perform_inversion(Job, old_eigenvectors)

        # Result
        assert all(np.array_equal(new_eigenvectors[i], expected_eigenvectors[i])
                   for i in range(len(expected_eigenvectors)))