1*7f8c6829SSascha Wildner /*-
2*7f8c6829SSascha Wildner * Copyright (c) 2008-2013 David Schultz <das@FreeBSD.org>
3*7f8c6829SSascha Wildner * All rights reserved.
4*7f8c6829SSascha Wildner *
5*7f8c6829SSascha Wildner * Redistribution and use in source and binary forms, with or without
6*7f8c6829SSascha Wildner * modification, are permitted provided that the following conditions
7*7f8c6829SSascha Wildner * are met:
8*7f8c6829SSascha Wildner * 1. Redistributions of source code must retain the above copyright
9*7f8c6829SSascha Wildner * notice, this list of conditions and the following disclaimer.
10*7f8c6829SSascha Wildner * 2. Redistributions in binary form must reproduce the above copyright
11*7f8c6829SSascha Wildner * notice, this list of conditions and the following disclaimer in the
12*7f8c6829SSascha Wildner * documentation and/or other materials provided with the distribution.
13*7f8c6829SSascha Wildner *
14*7f8c6829SSascha Wildner * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
15*7f8c6829SSascha Wildner * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
16*7f8c6829SSascha Wildner * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
17*7f8c6829SSascha Wildner * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
18*7f8c6829SSascha Wildner * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
19*7f8c6829SSascha Wildner * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
20*7f8c6829SSascha Wildner * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
21*7f8c6829SSascha Wildner * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
22*7f8c6829SSascha Wildner * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
23*7f8c6829SSascha Wildner * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
24*7f8c6829SSascha Wildner * SUCH DAMAGE.
25*7f8c6829SSascha Wildner *
26*7f8c6829SSascha Wildner * $FreeBSD: src/tools/regression/lib/msun/test-invctrig.c,v 1.1 2013/05/31 00:27:55 svnexp Exp $
27*7f8c6829SSascha Wildner */
28*7f8c6829SSascha Wildner
29*7f8c6829SSascha Wildner /*
30*7f8c6829SSascha Wildner * Tests for casin[h](), cacos[h](), and catan[h]().
31*7f8c6829SSascha Wildner */
32*7f8c6829SSascha Wildner
33*7f8c6829SSascha Wildner #include <assert.h>
34*7f8c6829SSascha Wildner #include <complex.h>
35*7f8c6829SSascha Wildner #include <fenv.h>
36*7f8c6829SSascha Wildner #include <float.h>
37*7f8c6829SSascha Wildner #include <math.h>
38*7f8c6829SSascha Wildner #include <stdio.h>
39*7f8c6829SSascha Wildner
40*7f8c6829SSascha Wildner #define ALL_STD_EXCEPT (FE_DIVBYZERO | FE_INEXACT | FE_INVALID | \
41*7f8c6829SSascha Wildner FE_OVERFLOW | FE_UNDERFLOW)
42*7f8c6829SSascha Wildner #define OPT_INVALID (ALL_STD_EXCEPT & ~FE_INVALID)
43*7f8c6829SSascha Wildner #define OPT_INEXACT (ALL_STD_EXCEPT & ~FE_INEXACT)
44*7f8c6829SSascha Wildner #define FLT_ULP() ldexpl(1.0, 1 - FLT_MANT_DIG)
45*7f8c6829SSascha Wildner #define DBL_ULP() ldexpl(1.0, 1 - DBL_MANT_DIG)
46*7f8c6829SSascha Wildner #define LDBL_ULP() ldexpl(1.0, 1 - LDBL_MANT_DIG)
47*7f8c6829SSascha Wildner
48*7f8c6829SSascha Wildner #pragma STDC FENV_ACCESS ON
49*7f8c6829SSascha Wildner #pragma STDC CX_LIMITED_RANGE OFF
50*7f8c6829SSascha Wildner
51*7f8c6829SSascha Wildner /* Flags that determine whether to check the signs of the result. */
52*7f8c6829SSascha Wildner #define CS_REAL 1
53*7f8c6829SSascha Wildner #define CS_IMAG 2
54*7f8c6829SSascha Wildner #define CS_BOTH (CS_REAL | CS_IMAG)
55*7f8c6829SSascha Wildner
56*7f8c6829SSascha Wildner #ifdef DEBUG
57*7f8c6829SSascha Wildner #define debug(...) printf(__VA_ARGS__)
58*7f8c6829SSascha Wildner #else
59*7f8c6829SSascha Wildner #define debug(...) (void)0
60*7f8c6829SSascha Wildner #endif
61*7f8c6829SSascha Wildner
62*7f8c6829SSascha Wildner /*
63*7f8c6829SSascha Wildner * Test that a function returns the correct value and sets the
64*7f8c6829SSascha Wildner * exception flags correctly. The exceptmask specifies which
65*7f8c6829SSascha Wildner * exceptions we should check. We need to be lenient for several
66*7f8c6829SSascha Wildner * reasons, but mainly because on some architectures it's impossible
67*7f8c6829SSascha Wildner * to raise FE_OVERFLOW without raising FE_INEXACT.
68*7f8c6829SSascha Wildner *
69*7f8c6829SSascha Wildner * These are macros instead of functions so that assert provides more
70*7f8c6829SSascha Wildner * meaningful error messages.
71*7f8c6829SSascha Wildner *
72*7f8c6829SSascha Wildner * XXX The volatile here is to avoid gcc's bogus constant folding and work
73*7f8c6829SSascha Wildner * around the lack of support for the FENV_ACCESS pragma.
74*7f8c6829SSascha Wildner */
75*7f8c6829SSascha Wildner #define test_p(func, z, result, exceptmask, excepts, checksign) do { \
76*7f8c6829SSascha Wildner volatile long double complex _d = z; \
77*7f8c6829SSascha Wildner debug(" testing %s(%Lg + %Lg I) == %Lg + %Lg I\n", #func, \
78*7f8c6829SSascha Wildner creall(_d), cimagl(_d), creall(result), cimagl(result)); \
79*7f8c6829SSascha Wildner assert(feclearexcept(FE_ALL_EXCEPT) == 0); \
80*7f8c6829SSascha Wildner assert(cfpequal((func)(_d), (result), (checksign))); \
81*7f8c6829SSascha Wildner assert(((func), fetestexcept(exceptmask) == (excepts))); \
82*7f8c6829SSascha Wildner } while (0)
83*7f8c6829SSascha Wildner
84*7f8c6829SSascha Wildner /*
85*7f8c6829SSascha Wildner * Test within a given tolerance. The tolerance indicates relative error
86*7f8c6829SSascha Wildner * in ulps.
87*7f8c6829SSascha Wildner */
88*7f8c6829SSascha Wildner #define test_p_tol(func, z, result, tol) do { \
89*7f8c6829SSascha Wildner volatile long double complex _d = z; \
90*7f8c6829SSascha Wildner debug(" testing %s(%Lg + %Lg I) ~= %Lg + %Lg I\n", #func, \
91*7f8c6829SSascha Wildner creall(_d), cimagl(_d), creall(result), cimagl(result)); \
92*7f8c6829SSascha Wildner assert(cfpequal_tol((func)(_d), (result), (tol))); \
93*7f8c6829SSascha Wildner } while (0)
94*7f8c6829SSascha Wildner
95*7f8c6829SSascha Wildner /* These wrappers apply the identities f(conj(z)) = conj(f(z)). */
96*7f8c6829SSascha Wildner #define test(func, z, result, exceptmask, excepts, checksign) do { \
97*7f8c6829SSascha Wildner test_p(func, z, result, exceptmask, excepts, checksign); \
98*7f8c6829SSascha Wildner test_p(func, conjl(z), conjl(result), exceptmask, excepts, checksign); \
99*7f8c6829SSascha Wildner } while (0)
100*7f8c6829SSascha Wildner #define test_tol(func, z, result, tol) do { \
101*7f8c6829SSascha Wildner test_p_tol(func, z, result, tol); \
102*7f8c6829SSascha Wildner test_p_tol(func, conjl(z), conjl(result), tol); \
103*7f8c6829SSascha Wildner } while (0)
104*7f8c6829SSascha Wildner
105*7f8c6829SSascha Wildner /* Test the given function in all precisions. */
106*7f8c6829SSascha Wildner #define testall(func, x, result, exceptmask, excepts, checksign) do { \
107*7f8c6829SSascha Wildner test(func, x, result, exceptmask, excepts, checksign); \
108*7f8c6829SSascha Wildner test(func##f, x, result, exceptmask, excepts, checksign); \
109*7f8c6829SSascha Wildner } while (0)
110*7f8c6829SSascha Wildner #define testall_odd(func, x, result, exceptmask, excepts, checksign) do { \
111*7f8c6829SSascha Wildner testall(func, x, result, exceptmask, excepts, checksign); \
112*7f8c6829SSascha Wildner testall(func, -(x), -result, exceptmask, excepts, checksign); \
113*7f8c6829SSascha Wildner } while (0)
114*7f8c6829SSascha Wildner #define testall_even(func, x, result, exceptmask, excepts, checksign) do { \
115*7f8c6829SSascha Wildner testall(func, x, result, exceptmask, excepts, checksign); \
116*7f8c6829SSascha Wildner testall(func, -(x), result, exceptmask, excepts, checksign); \
117*7f8c6829SSascha Wildner } while (0)
118*7f8c6829SSascha Wildner
119*7f8c6829SSascha Wildner /*
120*7f8c6829SSascha Wildner * Test the given function in all precisions, within a given tolerance.
121*7f8c6829SSascha Wildner * The tolerance is specified in ulps.
122*7f8c6829SSascha Wildner */
123*7f8c6829SSascha Wildner #define testall_tol(func, x, result, tol) do { \
124*7f8c6829SSascha Wildner test_tol(func, x, result, (tol) * DBL_ULP()); \
125*7f8c6829SSascha Wildner test_tol(func##f, x, result, (tol) * FLT_ULP()); \
126*7f8c6829SSascha Wildner } while (0)
127*7f8c6829SSascha Wildner #define testall_odd_tol(func, x, result, tol) do { \
128*7f8c6829SSascha Wildner testall_tol(func, x, result, tol); \
129*7f8c6829SSascha Wildner testall_tol(func, -(x), -result, tol); \
130*7f8c6829SSascha Wildner } while (0)
131*7f8c6829SSascha Wildner #define testall_even_tol(func, x, result, tol) do { \
132*7f8c6829SSascha Wildner testall_tol(func, x, result, tol); \
133*7f8c6829SSascha Wildner testall_tol(func, -(x), result, tol); \
134*7f8c6829SSascha Wildner } while (0)
135*7f8c6829SSascha Wildner
136*7f8c6829SSascha Wildner /* XXX static const long double */
137*7f8c6829SSascha Wildner #define pi 3.14159265358979323846264338327950280L
138*7f8c6829SSascha Wildner #define c3pi 9.42477796076937971538793014983850839L
139*7f8c6829SSascha Wildner
140*7f8c6829SSascha Wildner /*
141*7f8c6829SSascha Wildner * Determine whether x and y are equal, with two special rules:
142*7f8c6829SSascha Wildner * +0.0 != -0.0
143*7f8c6829SSascha Wildner * NaN == NaN
144*7f8c6829SSascha Wildner * If checksign is 0, we compare the absolute values instead.
145*7f8c6829SSascha Wildner */
146*7f8c6829SSascha Wildner static int
fpequal(long double x,long double y,int checksign)147*7f8c6829SSascha Wildner fpequal(long double x, long double y, int checksign)
148*7f8c6829SSascha Wildner {
149*7f8c6829SSascha Wildner if (isnan(x) && isnan(y))
150*7f8c6829SSascha Wildner return (1);
151*7f8c6829SSascha Wildner if (checksign)
152*7f8c6829SSascha Wildner return (x == y && !signbit(x) == !signbit(y));
153*7f8c6829SSascha Wildner else
154*7f8c6829SSascha Wildner return (fabsl(x) == fabsl(y));
155*7f8c6829SSascha Wildner }
156*7f8c6829SSascha Wildner
157*7f8c6829SSascha Wildner static int
fpequal_tol(long double x,long double y,long double tol)158*7f8c6829SSascha Wildner fpequal_tol(long double x, long double y, long double tol)
159*7f8c6829SSascha Wildner {
160*7f8c6829SSascha Wildner fenv_t env;
161*7f8c6829SSascha Wildner int ret;
162*7f8c6829SSascha Wildner
163*7f8c6829SSascha Wildner if (isnan(x) && isnan(y))
164*7f8c6829SSascha Wildner return (1);
165*7f8c6829SSascha Wildner if (!signbit(x) != !signbit(y))
166*7f8c6829SSascha Wildner return (0);
167*7f8c6829SSascha Wildner if (x == y)
168*7f8c6829SSascha Wildner return (1);
169*7f8c6829SSascha Wildner if (tol == 0 || y == 0.0)
170*7f8c6829SSascha Wildner return (0);
171*7f8c6829SSascha Wildner
172*7f8c6829SSascha Wildner /* Hard case: need to check the tolerance. */
173*7f8c6829SSascha Wildner feholdexcept(&env);
174*7f8c6829SSascha Wildner ret = fabsl(x - y) <= fabsl(y * tol);
175*7f8c6829SSascha Wildner fesetenv(&env);
176*7f8c6829SSascha Wildner return (ret);
177*7f8c6829SSascha Wildner }
178*7f8c6829SSascha Wildner
179*7f8c6829SSascha Wildner static int
cfpequal(long double complex x,long double complex y,int checksign)180*7f8c6829SSascha Wildner cfpequal(long double complex x, long double complex y, int checksign)
181*7f8c6829SSascha Wildner {
182*7f8c6829SSascha Wildner return (fpequal(creal(x), creal(y), checksign & CS_REAL)
183*7f8c6829SSascha Wildner && fpequal(cimag(x), cimag(y), checksign & CS_IMAG));
184*7f8c6829SSascha Wildner }
185*7f8c6829SSascha Wildner
186*7f8c6829SSascha Wildner static int
cfpequal_tol(long double complex x,long double complex y,long double tol)187*7f8c6829SSascha Wildner cfpequal_tol(long double complex x, long double complex y, long double tol)
188*7f8c6829SSascha Wildner {
189*7f8c6829SSascha Wildner return (fpequal_tol(creal(x), creal(y), tol)
190*7f8c6829SSascha Wildner && fpequal_tol(cimag(x), cimag(y), tol));
191*7f8c6829SSascha Wildner }
192*7f8c6829SSascha Wildner
193*7f8c6829SSascha Wildner
194*7f8c6829SSascha Wildner /* Tests for 0 */
195*7f8c6829SSascha Wildner void
test_zero(void)196*7f8c6829SSascha Wildner test_zero(void)
197*7f8c6829SSascha Wildner {
198*7f8c6829SSascha Wildner long double complex zero = CMPLXL(0.0, 0.0);
199*7f8c6829SSascha Wildner
200*7f8c6829SSascha Wildner testall_tol(cacosh, zero, CMPLXL(0.0, pi / 2), 1);
201*7f8c6829SSascha Wildner testall_tol(cacosh, -zero, CMPLXL(0.0, -pi / 2), 1);
202*7f8c6829SSascha Wildner testall_tol(cacos, zero, CMPLXL(pi / 2, -0.0), 1);
203*7f8c6829SSascha Wildner testall_tol(cacos, -zero, CMPLXL(pi / 2, 0.0), 1);
204*7f8c6829SSascha Wildner
205*7f8c6829SSascha Wildner testall_odd(casinh, zero, zero, ALL_STD_EXCEPT, 0, CS_BOTH);
206*7f8c6829SSascha Wildner testall_odd(casin, zero, zero, ALL_STD_EXCEPT, 0, CS_BOTH);
207*7f8c6829SSascha Wildner
208*7f8c6829SSascha Wildner testall_odd(catanh, zero, zero, ALL_STD_EXCEPT, 0, CS_BOTH);
209*7f8c6829SSascha Wildner testall_odd(catan, zero, zero, ALL_STD_EXCEPT, 0, CS_BOTH);
210*7f8c6829SSascha Wildner }
211*7f8c6829SSascha Wildner
212*7f8c6829SSascha Wildner /*
213*7f8c6829SSascha Wildner * Tests for NaN inputs.
214*7f8c6829SSascha Wildner */
215*7f8c6829SSascha Wildner void
test_nan()216*7f8c6829SSascha Wildner test_nan()
217*7f8c6829SSascha Wildner {
218*7f8c6829SSascha Wildner long double complex nan_nan = CMPLXL(NAN, NAN);
219*7f8c6829SSascha Wildner long double complex z;
220*7f8c6829SSascha Wildner
221*7f8c6829SSascha Wildner /*
222*7f8c6829SSascha Wildner * IN CACOSH CACOS CASINH CATANH
223*7f8c6829SSascha Wildner * NaN,NaN NaN,NaN NaN,NaN NaN,NaN NaN,NaN
224*7f8c6829SSascha Wildner * finite,NaN NaN,NaN* NaN,NaN* NaN,NaN* NaN,NaN*
225*7f8c6829SSascha Wildner * NaN,finite NaN,NaN* NaN,NaN* NaN,NaN* NaN,NaN*
226*7f8c6829SSascha Wildner * NaN,Inf Inf,NaN NaN,-Inf ?Inf,NaN ?0,pi/2
227*7f8c6829SSascha Wildner * +-Inf,NaN Inf,NaN NaN,?Inf +-Inf,NaN +-0,NaN
228*7f8c6829SSascha Wildner * +-0,NaN NaN,NaN* pi/2,NaN NaN,NaN* +-0,NaN
229*7f8c6829SSascha Wildner * NaN,0 NaN,NaN* NaN,NaN* NaN,0 NaN,NaN*
230*7f8c6829SSascha Wildner *
231*7f8c6829SSascha Wildner * * = raise invalid
232*7f8c6829SSascha Wildner */
233*7f8c6829SSascha Wildner z = nan_nan;
234*7f8c6829SSascha Wildner testall(cacosh, z, nan_nan, ALL_STD_EXCEPT, 0, 0);
235*7f8c6829SSascha Wildner testall(cacos, z, nan_nan, ALL_STD_EXCEPT, 0, 0);
236*7f8c6829SSascha Wildner testall(casinh, z, nan_nan, ALL_STD_EXCEPT, 0, 0);
237*7f8c6829SSascha Wildner testall(casin, z, nan_nan, ALL_STD_EXCEPT, 0, 0);
238*7f8c6829SSascha Wildner testall(catanh, z, nan_nan, ALL_STD_EXCEPT, 0, 0);
239*7f8c6829SSascha Wildner testall(catan, z, nan_nan, ALL_STD_EXCEPT, 0, 0);
240*7f8c6829SSascha Wildner
241*7f8c6829SSascha Wildner z = CMPLXL(0.5, NAN);
242*7f8c6829SSascha Wildner testall(cacosh, z, nan_nan, OPT_INVALID, 0, 0);
243*7f8c6829SSascha Wildner testall(cacos, z, nan_nan, OPT_INVALID, 0, 0);
244*7f8c6829SSascha Wildner testall(casinh, z, nan_nan, OPT_INVALID, 0, 0);
245*7f8c6829SSascha Wildner testall(casin, z, nan_nan, OPT_INVALID, 0, 0);
246*7f8c6829SSascha Wildner testall(catanh, z, nan_nan, OPT_INVALID, 0, 0);
247*7f8c6829SSascha Wildner testall(catan, z, nan_nan, OPT_INVALID, 0, 0);
248*7f8c6829SSascha Wildner
249*7f8c6829SSascha Wildner z = CMPLXL(NAN, 0.5);
250*7f8c6829SSascha Wildner testall(cacosh, z, nan_nan, OPT_INVALID, 0, 0);
251*7f8c6829SSascha Wildner testall(cacos, z, nan_nan, OPT_INVALID, 0, 0);
252*7f8c6829SSascha Wildner testall(casinh, z, nan_nan, OPT_INVALID, 0, 0);
253*7f8c6829SSascha Wildner testall(casin, z, nan_nan, OPT_INVALID, 0, 0);
254*7f8c6829SSascha Wildner testall(catanh, z, nan_nan, OPT_INVALID, 0, 0);
255*7f8c6829SSascha Wildner testall(catan, z, nan_nan, OPT_INVALID, 0, 0);
256*7f8c6829SSascha Wildner
257*7f8c6829SSascha Wildner z = CMPLXL(NAN, INFINITY);
258*7f8c6829SSascha Wildner testall(cacosh, z, CMPLXL(INFINITY, NAN), ALL_STD_EXCEPT, 0, CS_REAL);
259*7f8c6829SSascha Wildner testall(cacosh, -z, CMPLXL(INFINITY, NAN), ALL_STD_EXCEPT, 0, CS_REAL);
260*7f8c6829SSascha Wildner testall(cacos, z, CMPLXL(NAN, -INFINITY), ALL_STD_EXCEPT, 0, CS_IMAG);
261*7f8c6829SSascha Wildner testall(casinh, z, CMPLXL(INFINITY, NAN), ALL_STD_EXCEPT, 0, 0);
262*7f8c6829SSascha Wildner testall(casin, z, CMPLXL(NAN, INFINITY), ALL_STD_EXCEPT, 0, CS_IMAG);
263*7f8c6829SSascha Wildner testall_tol(catanh, z, CMPLXL(0.0, pi / 2), 1);
264*7f8c6829SSascha Wildner testall(catan, z, CMPLXL(NAN, 0.0), ALL_STD_EXCEPT, 0, CS_IMAG);
265*7f8c6829SSascha Wildner
266*7f8c6829SSascha Wildner z = CMPLXL(INFINITY, NAN);
267*7f8c6829SSascha Wildner testall_even(cacosh, z, CMPLXL(INFINITY, NAN), ALL_STD_EXCEPT, 0,
268*7f8c6829SSascha Wildner CS_REAL);
269*7f8c6829SSascha Wildner testall_even(cacos, z, CMPLXL(NAN, INFINITY), ALL_STD_EXCEPT, 0, 0);
270*7f8c6829SSascha Wildner testall_odd(casinh, z, CMPLXL(INFINITY, NAN), ALL_STD_EXCEPT, 0,
271*7f8c6829SSascha Wildner CS_REAL);
272*7f8c6829SSascha Wildner testall_odd(casin, z, CMPLXL(NAN, INFINITY), ALL_STD_EXCEPT, 0, 0);
273*7f8c6829SSascha Wildner testall_odd(catanh, z, CMPLXL(0.0, NAN), ALL_STD_EXCEPT, 0, CS_REAL);
274*7f8c6829SSascha Wildner testall_odd_tol(catan, z, CMPLXL(pi / 2, 0.0), 1);
275*7f8c6829SSascha Wildner
276*7f8c6829SSascha Wildner z = CMPLXL(0.0, NAN);
277*7f8c6829SSascha Wildner /* XXX We allow a spurious inexact exception here. */
278*7f8c6829SSascha Wildner testall_even(cacosh, z, nan_nan, OPT_INVALID & ~FE_INEXACT, 0, 0);
279*7f8c6829SSascha Wildner testall_even_tol(cacos, z, CMPLXL(pi / 2, NAN), 1);
280*7f8c6829SSascha Wildner testall_odd(casinh, z, nan_nan, OPT_INVALID, 0, 0);
281*7f8c6829SSascha Wildner testall_odd(casin, z, CMPLXL(0.0, NAN), ALL_STD_EXCEPT, 0, CS_REAL);
282*7f8c6829SSascha Wildner testall_odd(catanh, z, CMPLXL(0.0, NAN), OPT_INVALID, 0, CS_REAL);
283*7f8c6829SSascha Wildner testall_odd(catan, z, nan_nan, OPT_INVALID, 0, 0);
284*7f8c6829SSascha Wildner
285*7f8c6829SSascha Wildner z = CMPLXL(NAN, 0.0);
286*7f8c6829SSascha Wildner testall(cacosh, z, nan_nan, OPT_INVALID, 0, 0);
287*7f8c6829SSascha Wildner testall(cacos, z, nan_nan, OPT_INVALID, 0, 0);
288*7f8c6829SSascha Wildner testall(casinh, z, CMPLXL(NAN, 0), ALL_STD_EXCEPT, 0, CS_IMAG);
289*7f8c6829SSascha Wildner testall(casin, z, nan_nan, OPT_INVALID, 0, 0);
290*7f8c6829SSascha Wildner testall(catanh, z, nan_nan, OPT_INVALID, 0, CS_IMAG);
291*7f8c6829SSascha Wildner testall(catan, z, CMPLXL(NAN, 0.0), ALL_STD_EXCEPT, 0, 0);
292*7f8c6829SSascha Wildner }
293*7f8c6829SSascha Wildner
294*7f8c6829SSascha Wildner void
test_inf(void)295*7f8c6829SSascha Wildner test_inf(void)
296*7f8c6829SSascha Wildner {
297*7f8c6829SSascha Wildner long double complex z;
298*7f8c6829SSascha Wildner
299*7f8c6829SSascha Wildner /*
300*7f8c6829SSascha Wildner * IN CACOSH CACOS CASINH CATANH
301*7f8c6829SSascha Wildner * Inf,Inf Inf,pi/4 pi/4,-Inf Inf,pi/4 0,pi/2
302*7f8c6829SSascha Wildner * -Inf,Inf Inf,3pi/4 3pi/4,-Inf --- ---
303*7f8c6829SSascha Wildner * Inf,finite Inf,0 0,-Inf Inf,0 0,pi/2
304*7f8c6829SSascha Wildner * -Inf,finite Inf,pi pi,-Inf --- ---
305*7f8c6829SSascha Wildner * finite,Inf Inf,pi/2 pi/2,-Inf Inf,pi/2 0,pi/2
306*7f8c6829SSascha Wildner */
307*7f8c6829SSascha Wildner z = CMPLXL(INFINITY, INFINITY);
308*7f8c6829SSascha Wildner testall_tol(cacosh, z, CMPLXL(INFINITY, pi / 4), 1);
309*7f8c6829SSascha Wildner testall_tol(cacosh, -z, CMPLXL(INFINITY, -c3pi / 4), 1);
310*7f8c6829SSascha Wildner testall_tol(cacos, z, CMPLXL(pi / 4, -INFINITY), 1);
311*7f8c6829SSascha Wildner testall_tol(cacos, -z, CMPLXL(c3pi / 4, INFINITY), 1);
312*7f8c6829SSascha Wildner testall_odd_tol(casinh, z, CMPLXL(INFINITY, pi / 4), 1);
313*7f8c6829SSascha Wildner testall_odd_tol(casin, z, CMPLXL(pi / 4, INFINITY), 1);
314*7f8c6829SSascha Wildner testall_odd_tol(catanh, z, CMPLXL(0, pi / 2), 1);
315*7f8c6829SSascha Wildner testall_odd_tol(catan, z, CMPLXL(pi / 2, 0), 1);
316*7f8c6829SSascha Wildner
317*7f8c6829SSascha Wildner z = CMPLXL(INFINITY, 0.5);
318*7f8c6829SSascha Wildner /* XXX We allow a spurious inexact exception here. */
319*7f8c6829SSascha Wildner testall(cacosh, z, CMPLXL(INFINITY, 0), OPT_INEXACT, 0, CS_BOTH);
320*7f8c6829SSascha Wildner testall_tol(cacosh, -z, CMPLXL(INFINITY, -pi), 1);
321*7f8c6829SSascha Wildner testall(cacos, z, CMPLXL(0, -INFINITY), OPT_INEXACT, 0, CS_BOTH);
322*7f8c6829SSascha Wildner testall_tol(cacos, -z, CMPLXL(pi, INFINITY), 1);
323*7f8c6829SSascha Wildner testall_odd(casinh, z, CMPLXL(INFINITY, 0), OPT_INEXACT, 0, CS_BOTH);
324*7f8c6829SSascha Wildner testall_odd_tol(casin, z, CMPLXL(pi / 2, INFINITY), 1);
325*7f8c6829SSascha Wildner testall_odd_tol(catanh, z, CMPLXL(0, pi / 2), 1);
326*7f8c6829SSascha Wildner testall_odd_tol(catan, z, CMPLXL(pi / 2, 0), 1);
327*7f8c6829SSascha Wildner
328*7f8c6829SSascha Wildner z = CMPLXL(0.5, INFINITY);
329*7f8c6829SSascha Wildner testall_tol(cacosh, z, CMPLXL(INFINITY, pi / 2), 1);
330*7f8c6829SSascha Wildner testall_tol(cacosh, -z, CMPLXL(INFINITY, -pi / 2), 1);
331*7f8c6829SSascha Wildner testall_tol(cacos, z, CMPLXL(pi / 2, -INFINITY), 1);
332*7f8c6829SSascha Wildner testall_tol(cacos, -z, CMPLXL(pi / 2, INFINITY), 1);
333*7f8c6829SSascha Wildner testall_odd_tol(casinh, z, CMPLXL(INFINITY, pi / 2), 1);
334*7f8c6829SSascha Wildner /* XXX We allow a spurious inexact exception here. */
335*7f8c6829SSascha Wildner testall_odd(casin, z, CMPLXL(0.0, INFINITY), OPT_INEXACT, 0, CS_BOTH);
336*7f8c6829SSascha Wildner testall_odd_tol(catanh, z, CMPLXL(0, pi / 2), 1);
337*7f8c6829SSascha Wildner testall_odd_tol(catan, z, CMPLXL(pi / 2, 0), 1);
338*7f8c6829SSascha Wildner }
339*7f8c6829SSascha Wildner
340*7f8c6829SSascha Wildner /* Tests along the real and imaginary axes. */
341*7f8c6829SSascha Wildner void
test_axes(void)342*7f8c6829SSascha Wildner test_axes(void)
343*7f8c6829SSascha Wildner {
344*7f8c6829SSascha Wildner static const long double nums[] = {
345*7f8c6829SSascha Wildner -2, -1, -0.5, 0.5, 1, 2
346*7f8c6829SSascha Wildner };
347*7f8c6829SSascha Wildner long double complex z;
348*7f8c6829SSascha Wildner int i;
349*7f8c6829SSascha Wildner
350*7f8c6829SSascha Wildner for (i = 0; i < sizeof(nums) / sizeof(nums[0]); i++) {
351*7f8c6829SSascha Wildner /* Real axis */
352*7f8c6829SSascha Wildner z = CMPLXL(nums[i], 0.0);
353*7f8c6829SSascha Wildner if (fabs(nums[i]) <= 1) {
354*7f8c6829SSascha Wildner testall_tol(cacosh, z, CMPLXL(0.0, acos(nums[i])), 1);
355*7f8c6829SSascha Wildner testall_tol(cacos, z, CMPLXL(acosl(nums[i]), -0.0), 1);
356*7f8c6829SSascha Wildner testall_tol(casin, z, CMPLXL(asinl(nums[i]), 0.0), 1);
357*7f8c6829SSascha Wildner testall_tol(catanh, z, CMPLXL(atanh(nums[i]), 0.0), 1);
358*7f8c6829SSascha Wildner } else {
359*7f8c6829SSascha Wildner testall_tol(cacosh, z,
360*7f8c6829SSascha Wildner CMPLXL(acosh(fabs(nums[i])),
361*7f8c6829SSascha Wildner (nums[i] < 0) ? pi : 0), 1);
362*7f8c6829SSascha Wildner testall_tol(cacos, z,
363*7f8c6829SSascha Wildner CMPLXL((nums[i] < 0) ? pi : 0,
364*7f8c6829SSascha Wildner -acosh(fabs(nums[i]))), 1);
365*7f8c6829SSascha Wildner testall_tol(casin, z,
366*7f8c6829SSascha Wildner CMPLXL(copysign(pi / 2, nums[i]),
367*7f8c6829SSascha Wildner acosh(fabs(nums[i]))), 1);
368*7f8c6829SSascha Wildner testall_tol(catanh, z,
369*7f8c6829SSascha Wildner CMPLXL(atanh(1 / nums[i]), pi / 2), 1);
370*7f8c6829SSascha Wildner }
371*7f8c6829SSascha Wildner testall_tol(casinh, z, CMPLXL(asinh(nums[i]), 0.0), 1);
372*7f8c6829SSascha Wildner testall_tol(catan, z, CMPLXL(atan(nums[i]), 0), 1);
373*7f8c6829SSascha Wildner
374*7f8c6829SSascha Wildner /* TODO: Test the imaginary axis. */
375*7f8c6829SSascha Wildner }
376*7f8c6829SSascha Wildner }
377*7f8c6829SSascha Wildner
378*7f8c6829SSascha Wildner void
test_small(void)379*7f8c6829SSascha Wildner test_small(void)
380*7f8c6829SSascha Wildner {
381*7f8c6829SSascha Wildner /*
382*7f8c6829SSascha Wildner * z = 0.75 + i 0.25
383*7f8c6829SSascha Wildner * acos(z) = Pi/4 - i ln(2)/2
384*7f8c6829SSascha Wildner * asin(z) = Pi/4 + i ln(2)/2
385*7f8c6829SSascha Wildner * atan(z) = atan(4)/2 + i ln(17/9)/4
386*7f8c6829SSascha Wildner */
387*7f8c6829SSascha Wildner static const struct {
388*7f8c6829SSascha Wildner complex long double z;
389*7f8c6829SSascha Wildner complex long double acos_z;
390*7f8c6829SSascha Wildner complex long double asin_z;
391*7f8c6829SSascha Wildner complex long double atan_z;
392*7f8c6829SSascha Wildner } tests[] = {
393*7f8c6829SSascha Wildner { CMPLXL(0.75L, 0.25L),
394*7f8c6829SSascha Wildner CMPLXL(pi / 4, -0.34657359027997265470861606072908828L),
395*7f8c6829SSascha Wildner CMPLXL(pi / 4, 0.34657359027997265470861606072908828L),
396*7f8c6829SSascha Wildner CMPLXL(0.66290883183401623252961960521423782L,
397*7f8c6829SSascha Wildner 0.15899719167999917436476103600701878L) },
398*7f8c6829SSascha Wildner };
399*7f8c6829SSascha Wildner int i;
400*7f8c6829SSascha Wildner
401*7f8c6829SSascha Wildner for (i = 0; i < sizeof(tests) / sizeof(tests[0]); i++) {
402*7f8c6829SSascha Wildner testall_tol(cacos, tests[i].z, tests[i].acos_z, 2);
403*7f8c6829SSascha Wildner testall_odd_tol(casin, tests[i].z, tests[i].asin_z, 2);
404*7f8c6829SSascha Wildner testall_odd_tol(catan, tests[i].z, tests[i].atan_z, 2);
405*7f8c6829SSascha Wildner }
406*7f8c6829SSascha Wildner }
407*7f8c6829SSascha Wildner
408*7f8c6829SSascha Wildner /* Test inputs that might cause overflow in a sloppy implementation. */
409*7f8c6829SSascha Wildner void
test_large(void)410*7f8c6829SSascha Wildner test_large(void)
411*7f8c6829SSascha Wildner {
412*7f8c6829SSascha Wildner
413*7f8c6829SSascha Wildner /* TODO: Write these tests */
414*7f8c6829SSascha Wildner }
415*7f8c6829SSascha Wildner
416*7f8c6829SSascha Wildner int
main(int argc,char * argv[])417*7f8c6829SSascha Wildner main(int argc, char *argv[])
418*7f8c6829SSascha Wildner {
419*7f8c6829SSascha Wildner
420*7f8c6829SSascha Wildner printf("1..6\n");
421*7f8c6829SSascha Wildner
422*7f8c6829SSascha Wildner test_zero();
423*7f8c6829SSascha Wildner printf("ok 1 - invctrig zero\n");
424*7f8c6829SSascha Wildner
425*7f8c6829SSascha Wildner test_nan();
426*7f8c6829SSascha Wildner printf("ok 2 - invctrig nan\n");
427*7f8c6829SSascha Wildner
428*7f8c6829SSascha Wildner test_inf();
429*7f8c6829SSascha Wildner printf("ok 3 - invctrig inf\n");
430*7f8c6829SSascha Wildner
431*7f8c6829SSascha Wildner test_axes();
432*7f8c6829SSascha Wildner printf("ok 4 - invctrig axes\n");
433*7f8c6829SSascha Wildner
434*7f8c6829SSascha Wildner test_small();
435*7f8c6829SSascha Wildner printf("ok 5 - invctrig small\n");
436*7f8c6829SSascha Wildner
437*7f8c6829SSascha Wildner test_large();
438*7f8c6829SSascha Wildner printf("ok 6 - invctrig large\n");
439*7f8c6829SSascha Wildner
440*7f8c6829SSascha Wildner return (0);
441*7f8c6829SSascha Wildner }
442