Initial commit.

- Added support for siesta Hamiltonians
 * In particular the TSHS has gained a lot of attention as
   that has all information at once. I have tested the Hamiltonian
   creation vs. Inelastica with great success.
 * The HS and HSX formats are also read (but currently they
   need the xa and lasto arrays from the full simulation
   to be downscaled in requirements
 * We can also read in the density matrix.