Move TB tests to own object
More careful treatment of units. Test forces.
Replace pos -> x
Forces working at large M
Tweak test.
Refactor Tb code.
Fix lattice initialization. Choose default for
rcut in Si potential. Get TbMD code to compile.
Apply shift to GoodwinSi Hamiltonian.
Add calculation of density function and its
integral.