2 * (c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
3 * See the copyright notice in the ACK home directory, in the file "Copyright".
5 * Author: Ceriel J.H. Jacobs
7 /* $Id: asin.c,v 1.4 1994/06/24 11:43:08 ceriel Exp $ */
11 #include "localmath.h"
14 asin_acos(double x, int cosfl)
20 -0.27368494524164255994e+2,
21 0.57208227877891731407e+2,
22 -0.39688862997540877339e+2,
23 0.10152522233806463645e+2,
24 -0.69674573447350646411e+0
27 -0.16421096714498560795e+3,
28 0.41714430248260412556e+3,
29 -0.38186303361750149284e+3,
30 0.15095270841030604719e+3,
31 -0.23823859153670238830e+2,
54 /* ??? avoid underflow ??? */
58 x += x * g * POLYNOM4(g, p) / POLYNOM5(g, q);
60 if (! negative) x = -x;
62 if ((cosfl == 0) == (i == 1)) {
63 x = (x + M_PI_4) + M_PI_4;
65 else if (cosfl && negative && i == 1) {
66 x = (x + M_PI_2) + M_PI_2;
68 if (! cosfl && negative) x = -x;
75 return asin_acos(x, 0);
81 return asin_acos(x, 1);