Almost done with surrogate Hamiltonian method. Added complex matrix class, chebyshev evaluation of time evolution operator where Hamiltonian is described in a matrix class. Other small additions and bug fixes too.