[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Re: Bug-gsl Digest, Vol 99, Issue 1
From: |
James Ward |
Subject: |
[Bug-gsl] Re: Bug-gsl Digest, Vol 99, Issue 1 |
Date: |
Tue, 5 Apr 2011 08:12:36 -0400 |
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, address@hidden <address@hidden> wrote:
> Send Bug-gsl mailing list submissions to
> address@hidden
>
> 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
> address@hidden
>
> You can reach the person managing the list at
> address@hidden
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Bug-gsl digest..."
>
- [Bug-gsl] Re: Bug-gsl Digest, Vol 99, Issue 1,
James Ward <=