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