Add user editable threshold for selecting states of the initial Hamiltonian

[Marius Retegan]

At the moment the threshold is set to the machine epsilon for double precision, which is very conservative.