Generalize loss interpolation
With #52 (closed), we obtain a better universality for user-provided loss functions. However right now we always split the loss of the parent interval proportionally into the children intervals. This prevents the user from guaranteeing that the intervals never become shorter than a certain size (e.g. machine precision) by means of merely redefining the loss.
I am not quite sure how we should address this though.