Lots of work on computing Jacobians and cleaning up assignments.  Still to do is getting the Hamiltonians right and the init method, then on to boilerplate.