Coupling intervals are specified at the beginning of the simulation along with total number of simulations. bar_main partitions lambdas so that they would be equidistant with respect to the ddG function, however the new lambda matrix needs to fit in the coupling interval. This might not work as expected.