Improve failure rate when finding the metric and eigen-directions
Occasionally, simple-pe
fails to iteratively update the metric and eigen-directions because of the error: f(a) and f(b) must have different signs
. This stems from the function scipy.optimize.root_scalar
having a restrictive interval which brackets the root. This MR iteratively increases the bracket until the root can be found.
Edited by Charlie Hoy