Mathematically, they're equivalent, but what's the point of defining z at all? It's only used in that one sin. It looks like the author:
else if(x > -0.5/(GSL_DBL_EPSILON*M_PI)) is trying to avoid some kind of overflow/underflow. I'm guessing: double s = sin(M_PI*z); is correct ... but I'm no expert. Best wishes, Jim Ward On 4/1/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
