/* Specification. */
#include <math.h>
-/* s_atanl.c
+/* s_atanl.c
*
- * Inverse circular tangent for 128-bit long double precision
+ * Inverse circular tangent for 128-bit long double precision
* (arctangent)
*
*
*
*/
-#include "isnanl.h"
-
/* arctan(k/8), k = 0, ..., 82 */
static const long double atantbl[84] = {
0.0000000000000000000000000000000000000000E0L,
{
/* Infinity. */
if (sign)
- return -atantbl[83];
+ return -atantbl[83];
else
- return atantbl[83];
+ return atantbl[83];
}
if (sign)
else
{
/* Index of nearest table element.
- Roundoff to integer is asymmetrical to avoid cancellation when t < 0
+ Roundoff to integer is asymmetrical to avoid cancellation when t < 0
(cf. fdlibm). */
k = 8.0 * x + 0.25;
u = 0.125 * k;