Description

Tests have shown that this leads to better convergence than what we did in the code before. Moreover, now the behavior is uniform for the different eigensolvers.

News snippet

Do subspace diagonalization after eigensolver in SCF for better convergence

Checklist

I have checked that my code follows the Octopus coding standards
I have added tests for all the new features added in this request.