roots-fortran icon indicating copy to clipboard operation
roots-fortran copied to clipboard

Published 20 hours ago • verified publisher icon jacobwilliams

Robust bisect

Open ivan-pi opened this issue 1 year ago • 4 comments

https://github.com/jacobwilliams/roots-fortran/blob/60524bc4f491a50ae1b2e538e1e486551812b9db/src/root_module.F90#L2524

The excerpt below is taken from Kahaner, Moler, & Nash, Numerical Methods and Software, 1989, Prentice-Hall, Inc.:

image

According to Herbie this rearrangement seems irrelevant in FP64, however I couldn't find an option to see what happens in lower precision: https://herbie.uwplse.org/demo/ebc6120a0224e43ec4e96ae4c79d8bf1537cfdf8.1.6/graph.html#history

ivan-pi avatar Aug 04 '22 12:08 ivan-pi