/* Test whether a basic csinl invocation works. */ /* This test is generated by misc/genbasic.c. */ #include #include #include #include #include "../basic.h" #pragma STDC FENV_ACCESS ON #define MF_UNSPEC1 (1 << 0) #define MF_UNSPEC2 (1 << 1) #define MF_MAYERR (1 << 2) #define MF_ANYSIGN1 (1 << 2) #define MF_ANYSIGN2 (1 << 3) // Soft fail on rounding errors and report only one. int imprecise; void test(int variant, long double complex input1, long double complex lower, long double complex expected, long double complex upper, int flags) { long double complex output = csinl(input1); if ( !(flags & MF_UNSPEC1) ) { long double real = creall(output); real = (flags & MF_ANYSIGN1) && real < 0.0 ? -real : real; long double lower_real = creall(lower); long double expected_real = creall(expected); long double upper_real = creall(upper); if ( !(isnan(expected_real) ? isnan(real) : isfinite(expected_real) && expected_real != 0.0 ? isfinite(real) && (real == expected_real || (lower_real < real && real < upper_real)) : real == expected_real) ) { if ( imprecise && isfinite(real) && isfinite(expected_real) ) return; warnx("(%d.) csinl(%.4Lf + i*%.4Lf).real = %.17La, not %.17La, diff %.17La, ratio %.20Lg", variant, creall(input1), cimagl(input1), real, expected_real, real - expected_real, real / expected_real); if ( !isfinite(real) || !isfinite(expected_real) ) exit(1); imprecise = 1; } long double imag = cimagl(output); imag = (flags & MF_ANYSIGN2) && imag < 0.0 ? -imag : imag; long double lower_imag = cimagl(lower); long double expected_imag = cimagl(expected); long double upper_imag = cimagl(upper); if ( !(isnan(expected_imag) ? isnan(imag) : isfinite(expected_imag) && expected_imag != 0.0 ? isfinite(imag) && (imag == expected_imag || (lower_imag < imag && imag < upper_imag)) : imag == expected_imag) ) { if ( imprecise && isfinite(imag) && isfinite(expected_imag) ) return; warnx("(%d.) csinl(%.4Lf + i*%.4Lf).imag = %.17La, not %.17La, diff %.17La, ratio %.20Lg", variant, creall(input1), cimagl(input1), imag, expected_imag, imag - expected_imag, imag / expected_imag); if ( !isfinite(imag) || !isfinite(expected_imag) ) exit(1); imprecise = 1; } } } int main(void) { test(1, CMPLXL(90.01, 13.37), CMPLXL(0x8.b167dccac6226c3p+15L, -0x8.eebaef0e510d094p+14L), CMPLXL(0x8.b167dccac6226c4p+15L, -0x8.eebaef0e510d093p+14L), CMPLXL(0x8.b167dccac6226c5p+15L, -0x8.eebaef0e510d092p+14L), 0); test(2, CMPLXL(-12.34, 13.37), CMPLXL(0x8.c62a6f4794b12ecp+13L, 0x9.861c6eb2b110b7p+15L), CMPLXL(0x8.c62a6f4794b12edp+13L, 0x9.861c6eb2b110b71p+15L), CMPLXL(0x8.c62a6f4794b12eep+13L, 0x9.861c6eb2b110b72p+15L), 0); test(3, CMPLXL(nanl(""), 13.37), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(4, CMPLXL(strtold("inf", NULL), 13.37), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(5, CMPLXL(strtold("-inf", NULL), 13.37), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(6, CMPLXL(0.0, 13.37), CMPLXL(-0x0.000000000000001p-16385L, 0x9.c5f12e7d5c95684p+15L), CMPLXL(0x0.0p+0L, 0x9.c5f12e7d5c95685p+15L), CMPLXL(0x0.000000000000001p-16385L, 0x9.c5f12e7d5c95686p+15L), 0); test(7, CMPLXL(90.01, -12.34), CMPLXL(0xc.69f0fa3b99ee0e6p+13L, 0xc.c183b03c4faef5dp+12L), CMPLXL(0xc.69f0fa3b99ee0e7p+13L, 0xc.c183b03c4faef5ep+12L), CMPLXL(0xc.69f0fa3b99ee0e8p+13L, 0xc.c183b03c4faef5fp+12L), 0); test(8, CMPLXL(-12.34, -12.34), CMPLXL(0xc.879655a638d024p+11L, -0xd.99b0c36ebd8895bp+13L), CMPLXL(0xc.879655a638d0241p+11L, -0xd.99b0c36ebd8895ap+13L), CMPLXL(0xc.879655a638d0242p+11L, -0xd.99b0c36ebd88959p+13L), 0); test(9, CMPLXL(nanl(""), -12.34), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(10, CMPLXL(strtold("inf", NULL), -12.34), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(11, CMPLXL(strtold("-inf", NULL), -12.34), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(12, CMPLXL(0.0, -12.34), CMPLXL(-0x0.000000000000001p-16385L, -0xd.f4d7cee6d57be0ep+13L), CMPLXL(0x0.0p+0L, -0xd.f4d7cee6d57be0dp+13L), CMPLXL(0x0.000000000000001p-16385L, -0xd.f4d7cee6d57be0cp+13L), 0); test(13, CMPLXL(90.01, nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(14, CMPLXL(-12.34, nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(15, CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(16, CMPLXL(strtold("inf", NULL), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(17, CMPLXL(strtold("-inf", NULL), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), CMPLXL(nanl(""), nanl("")), 0); test(18, CMPLXL(0.0, nanl("")), CMPLXL(-0x0.000000000000001p-16385L, nanl("")), CMPLXL(0x0.0p+0L, nanl("")), CMPLXL(0x0.000000000000001p-16385L, nanl("")), 0); test(19, CMPLXL(90.01, strtold("inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), 0); test(20, CMPLXL(-12.34, strtold("inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), 0); test(21, CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), 0 | MF_ANYSIGN2); test(22, CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), 0 | MF_ANYSIGN2); test(23, CMPLXL(strtold("-inf", NULL), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), 0 | MF_ANYSIGN2); test(24, CMPLXL(0.0, strtold("inf", NULL)), CMPLXL(-0x0.000000000000001p-16385L, strtold("inf", NULL)), CMPLXL(0x0.0p+0L, strtold("inf", NULL)), CMPLXL(0x0.000000000000001p-16385L, strtold("inf", NULL)), 0); test(25, CMPLXL(90.01, strtold("-inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("inf", NULL)), 0); test(26, CMPLXL(-12.34, strtold("-inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), 0); test(27, CMPLXL(nanl(""), strtold("-inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), 0 | MF_ANYSIGN2); test(28, CMPLXL(strtold("inf", NULL), strtold("-inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), 0 | MF_ANYSIGN2); test(29, CMPLXL(strtold("-inf", NULL), strtold("-inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), CMPLXL(nanl(""), strtold("inf", NULL)), 0 | MF_ANYSIGN2); test(30, CMPLXL(0.0, strtold("-inf", NULL)), CMPLXL(-0x0.000000000000001p-16385L, strtold("-inf", NULL)), CMPLXL(0x0.0p+0L, strtold("-inf", NULL)), CMPLXL(0x0.000000000000001p-16385L, strtold("-inf", NULL)), 0); test(31, CMPLXL(90.01, 0.0), CMPLXL(0xe.3b4642ba4986927p-4L, -0x0.000000000000001p-16385L), CMPLXL(0xe.3b4642ba4986928p-4L, -0x0.0p+0L), CMPLXL(0xe.3b4642ba4986929p-4L, 0x0.000000000000001p-16385L), 0); test(32, CMPLXL(-12.34, 0.0), CMPLXL(0xe.5d42e58f5a5c8abp-6L, -0x0.000000000000001p-16385L), CMPLXL(0xe.5d42e58f5a5c8acp-6L, 0x0.0p+0L), CMPLXL(0xe.5d42e58f5a5c8adp-6L, 0x0.000000000000001p-16385L), 0); test(33, CMPLXL(nanl(""), 0.0), CMPLXL(nanl(""), -0x0.000000000000001p-16385L), CMPLXL(nanl(""), 0x0.0p+0L), CMPLXL(nanl(""), 0x0.000000000000001p-16385L), 0); test(34, CMPLXL(strtold("inf", NULL), 0.0), CMPLXL(nanl(""), -0x0.000000000000001p-16385L), CMPLXL(nanl(""), 0x0.0p+0L), CMPLXL(nanl(""), 0x0.000000000000001p-16385L), 0); test(35, CMPLXL(strtold("-inf", NULL), 0.0), CMPLXL(nanl(""), -0x0.000000000000001p-16385L), CMPLXL(nanl(""), 0x0.0p+0L), CMPLXL(nanl(""), 0x0.000000000000001p-16385L), 0); test(36, CMPLXL(0.0, 0.0), CMPLXL(-0x0.000000000000001p-16385L, -0x0.000000000000001p-16385L), CMPLXL(0x0.0p+0L, 0x0.0p+0L), CMPLXL(0x0.000000000000001p-16385L, 0x0.000000000000001p-16385L), 0); return imprecise; }