`quantile_bisect` fails to converge for intervals far from zero

[nworb-cire]

I have the following code:

d = MixtureModel(
    [
        Uniform{Float64}(-0.0001, 0.0001),
        LogNormal{Float64}(11.174347445936371, 1.6086247197750911)
    ],
    [
        0.8832,
        0.1168
    ]
)
quantile(d, 0.99)

When quantile_bisect is run, it eventually settles into lx::Float64 == 643851.9074493216 and rx::Float64 == 643851.9074493217. Since rx == nextfloat(lx) and rx - lx > tol, convergence criteria are never met and it loops infinitely.

[nworb-cire]

If lx isa AbstractFloat then it should suffice to add a check on this line whether rx != nextfloat(lx).

[devmotion]

Fixed by #2070