With submodules and procedure pointers it's much easier to than with objects to put all hamiltonians in one binary.
