Davidson diagonalization in Fortran

• Follows modern Fortran practices
• Allows both dense and matrix-free calculations.
• Implements different correction methods

The Davidson method is suitable for diagonal-dominant symmetric matrices, that are quite common
in certain scientific problems like electronic structure. The Davidson method could be not practical
for other kinds of symmetric matrices.