The preconditioned residual vector is calculated for each band resulting in a 2*NBANDS basis-set
Within this subspace the NBANDS lowest eigenfunctions are calculated solving the eigenvalue problem
The NBANDS lowest eigenfunctions are used in the next step.