Add support for NNN hamiltonians to dynamic expansion (uniform case).
Also allow the use to supply BB to take_step().