Looking at the "kernel_sinl" function in "sincosl.c" file I found that the limit value for choosing the computation method is wrong with respect of the remaining code:

kernel_sinl (long double x, long double y, int iy)
{
 ...
 if (x < 0.148375L)        /* |x| < 0.1484375 */

The correct line should be:

 if (x < 0.1484375L)        /* |x| < 0.1484375 */

Probably a better approach would be defining this limit as a macro.

Regards,

Bob


Reply via email to