| 1 | // RUN: %clang_builtins %s %librt -o %t && %run %t |
|---|
| 2 | |
|---|
| 3 | #define SINGLE_PRECISION |
|---|
| 4 | #include <fenv.h> |
|---|
| 5 | #include <math.h> |
|---|
| 6 | #include <stdio.h> |
|---|
| 7 | #include "fp_lib.h" |
|---|
| 8 | |
|---|
| 9 | int test__compiler_rt_fmaxf(fp_t x, fp_t y) { |
|---|
| 10 | fp_t crt_value = __compiler_rt_fmaxf(x, y); |
|---|
| 11 | fp_t libm_value = fmaxf(x: x, y: y); |
|---|
| 12 | // Consider +0 and -0 equal, and also disregard the sign/payload of two NaNs. |
|---|
| 13 | if (crt_value != libm_value && |
|---|
| 14 | !(crt_isnan(crt_value) && crt_isnan(libm_value))) { |
|---|
| 15 | printf(format: "error: in __compiler_rt_fmaxf(%a [%X], %a [%X]) = %a [%X] " |
|---|
| 16 | "!= %a [%X]\n" , |
|---|
| 17 | x, toRep(x), y, toRep(x: y), crt_value, toRep(x: crt_value), libm_value, |
|---|
| 18 | toRep(x: libm_value)); |
|---|
| 19 | return 1; |
|---|
| 20 | } |
|---|
| 21 | return 0; |
|---|
| 22 | } |
|---|
| 23 | |
|---|
| 24 | fp_t cases[] = { |
|---|
| 25 | -NAN, NAN, -INFINITY, INFINITY, -0.0, 0.0, -1, 1, -2, 2, |
|---|
| 26 | -0x1.0p-127, 0x1.0p-127, -0x1.0p-128, 0x1.0p-128, // subnormals |
|---|
| 27 | -1.001, 1.001, -1.002, 1.002, |
|---|
| 28 | }; |
|---|
| 29 | |
|---|
| 30 | int main() { |
|---|
| 31 | const unsigned N = sizeof(cases) / sizeof(cases[0]); |
|---|
| 32 | unsigned i, j; |
|---|
| 33 | for (i = 0; i < N; ++i) { |
|---|
| 34 | for (j = 0; j < N; ++j) { |
|---|
| 35 | if (test__compiler_rt_fmaxf(x: cases[i], y: cases[j])) return 1; |
|---|
| 36 | } |
|---|
| 37 | } |
|---|
| 38 | return 0; |
|---|
| 39 | } |
|---|
| 40 | |
|---|
