Minor optimization in bisection_iter

This is a tiny line that I had on cosma but didn't commit before cosma had its downtime.

It returns the minimum energy, when the lower bound is smaller or equal to the minimum energy and there is still a net cooling (LambdaNet_cgs < 0).

Edit: Sorry, merge request went into SWIFTsim instead of the colibre version.