Consider to upgrade to the toms748 algorithm

[ibell]

One example here:

sgtpy/sgtpy/vrmie_pure/density_solver.py

Line 104 in 5206013

rho_ext = brentq(dPsaft_fun, bracket[0], bracket[1],

.

The toms748 algorithm has all the guarantees of brentq but is always faster according to their derivations in the paper references in the scipy docs, hence no reason to not use it I think.

[gustavochm]

Hey Ian,

Thank you for raising this issue. You are right. I'll update the density solvers in the next version.

Gustavo