Refining some details:
* bisectY only if lb(y)+1e-8 < loup_est<ub(y)-0.1*diam(y))
* abrupt_decline_strict: caseB: alpha=0, step=0 (no more estimations)