1/*	$OpenBSD: conj_test.c,v 1.4 2021/12/13 18:04:28 deraadt Exp $	*/
2/*-
3 * Copyright (c) 2008 David Schultz <das@FreeBSD.org>
4 * All rights reserved.
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 * 1. Redistributions of source code must retain the above copyright
10 *    notice, this list of conditions and the following disclaimer.
11 * 2. Redistributions in binary form must reproduce the above copyright
12 *    notice, this list of conditions and the following disclaimer in the
13 *    documentation and/or other materials provided with the distribution.
14 *
15 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
16 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18 * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
19 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
20 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
21 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
22 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
23 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
24 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
25 * SUCH DAMAGE.
26 */
27
28/*
29 * Tests for conj{,f,l}()
30 */
31
32#include <complex.h>
33#include <fenv.h>
34#include <math.h>
35#include <stdio.h>
36
37#include "test-utils.h"
38
39#pragma	STDC CX_LIMITED_RANGE	OFF
40
41static float complex (*libconjf)(float complex) = conjf;
42static double complex (*libconj)(double complex) = conj;
43static long double complex (*libconjl)(long double complex) = conjl;
44static float (*libcrealf)(float complex) = crealf;
45static double (*libcreal)(double complex) = creal;
46static long double (*libcreall)(long double complex) = creall;
47static float (*libcimagf)(float complex) = cimagf;
48static double (*libcimag)(double complex) = cimag;
49static long double (*libcimagl)(long double complex) = cimagl;
50
51static const double tests[] = {
52	/* a +  bI */
53	0.0,	0.0,
54	0.0,	1.0,
55	1.0,	0.0,
56	-1.0,	0.0,
57	1.0,	-0.0,
58	0.0,	-1.0,
59	2.0,	4.0,
60	0.0,	INFINITY,
61	0.0,	-INFINITY,
62	INFINITY, 0.0,
63	NAN,	1.0,
64	1.0,	NAN,
65	NAN,	NAN,
66	-INFINITY, INFINITY,
67};
68
69ATF_TC_WITHOUT_HEAD(main);
70ATF_TC_BODY(main, tc)
71{
72	static const int ntests = sizeof(tests) / sizeof(tests[0]) / 2;
73	complex float in;
74	complex long double expected;
75	int i;
76
77	for (i = 0; i < ntests; i++) {
78		__real__ expected = __real__ in = tests[2 * i];
79		__imag__ in = tests[2 * i + 1];
80		__imag__ expected = -cimag(in);
81
82		ATF_REQUIRE(fpequal_cs(libcrealf(in), __real__ in, true));
83		ATF_REQUIRE(fpequal_cs(libcreal(in), __real__ in, true));
84		ATF_REQUIRE(fpequal_cs(libcreall(in), __real__ in, true));
85		ATF_REQUIRE(fpequal_cs(libcimagf(in), __imag__ in, true));
86		ATF_REQUIRE(fpequal_cs(libcimag(in), __imag__ in, true));
87		ATF_REQUIRE(fpequal_cs(libcimagl(in), __imag__ in, true));
88
89		ATF_REQUIRE_EQ(0, feclearexcept(FE_ALL_EXCEPT));
90		ATF_REQUIRE_MSG(
91		    cfpequal(libconjf(in), expected),
92		    "conjf(%#.2g + %#.2gI): wrong value", creal(in), cimag(in)
93		);
94		ATF_REQUIRE_EQ_MSG(0, fetestexcept(FE_ALL_EXCEPT),
95		    "conj(%#.2g + %#.2gI): threw an exception: %#x", creal(in),
96		    cimag(in), fetestexcept(FE_ALL_EXCEPT));
97
98		ATF_REQUIRE_EQ(0, feclearexcept(FE_ALL_EXCEPT));
99		ATF_REQUIRE_MSG(cfpequal(libconj(in), expected),
100		    "conj(%#.2g + %#.2gI): wrong value", creal(in), cimag(in));
101		ATF_REQUIRE_EQ_MSG(0, fetestexcept(FE_ALL_EXCEPT),
102		    "conj(%#.2g + %#.2gI): threw an exception: %#x", creal(in),
103		    cimag(in), fetestexcept(FE_ALL_EXCEPT));
104
105		ATF_REQUIRE_EQ(0, feclearexcept(FE_ALL_EXCEPT));
106		ATF_REQUIRE_MSG(cfpequal(libconjl(in), expected),
107		    "conjl(%#.2g + %#.2gI): wrong value", creal(in), cimag(in));
108		ATF_REQUIRE_EQ_MSG(0, fetestexcept(FE_ALL_EXCEPT),
109		    "conjl(%#.2g + %#.2gI): threw an exception: %#x", creal(in),
110		    cimag(in), fetestexcept(FE_ALL_EXCEPT));
111 	}
112}
113
114ATF_TP_ADD_TCS(tp)
115{
116	ATF_TP_ADD_TC(tp, main);
117
118	return (atf_no_error());
119}
120