Added Andrea's cpp code that computes derivaives of Hamiltonian in different coordinate and reference frames