I did some timing tests, and it looks like Moler & Morrision is about twice as slow as gsl_hypot(), up to the same accuracy (1e-15).
http://steve.hollasch.net/cgindex/math/pythag-root.txt suggests this is due to M&M being cubically convergent, while gsl_hypot() is quadratically convergent. So ... never mind! Jim Ward On 4/12/11, [email protected] <[email protected]> wrote: > Send Bug-gsl mailing list submissions to > [email protected] > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.gnu.org/mailman/listinfo/bug-gsl > or, via email, send a message with subject or body 'help' to > [email protected] > > You can reach the person managing the list at > [email protected] > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Bug-gsl digest..." > _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
