| 1 | // RUN: %clang_builtins %s %librt -o %t && %run %t |
|---|
| 2 | // REQUIRES: librt_has_divdf3 |
|---|
| 3 | |
|---|
| 4 | #include "int_lib.h" |
|---|
| 5 | #include <stdio.h> |
|---|
| 6 | |
|---|
| 7 | #include "fp_test.h" |
|---|
| 8 | |
|---|
| 9 | // Returns: a / b |
|---|
| 10 | COMPILER_RT_ABI double __divdf3(double a, double b); |
|---|
| 11 | |
|---|
| 12 | int test__divdf3(double a, double b, uint64_t expected) |
|---|
| 13 | { |
|---|
| 14 | double x = __divdf3(a, b); |
|---|
| 15 | int ret = compareResultD(result: x, expected); |
|---|
| 16 | |
|---|
| 17 | if (ret){ |
|---|
| 18 | printf(format: "error in test__divdf3(%.20e, %.20e) = %.20e, " |
|---|
| 19 | "expected %.20e\n" , a, b, x, |
|---|
| 20 | fromRep64(x: expected)); |
|---|
| 21 | } |
|---|
| 22 | return ret; |
|---|
| 23 | } |
|---|
| 24 | |
|---|
| 25 | int main() |
|---|
| 26 | { |
|---|
| 27 | // Returned NaNs are assumed to be qNaN by default |
|---|
| 28 | |
|---|
| 29 | // qNaN / any = qNaN |
|---|
| 30 | if (test__divdf3(a: makeQNaN64(), b: 3., UINT64_C(0x7ff8000000000000))) |
|---|
| 31 | return 1; |
|---|
| 32 | // NaN / any = NaN |
|---|
| 33 | if (test__divdf3(a: makeNaN64(UINT64_C(0x123)), b: 3., UINT64_C(0x7ff8000000000000))) |
|---|
| 34 | return 1; |
|---|
| 35 | // any / qNaN = qNaN |
|---|
| 36 | if (test__divdf3(a: 3., b: makeQNaN64(), UINT64_C(0x7ff8000000000000))) |
|---|
| 37 | return 1; |
|---|
| 38 | // any / NaN = NaN |
|---|
| 39 | if (test__divdf3(a: 3., b: makeNaN64(UINT64_C(0x123)), UINT64_C(0x7ff8000000000000))) |
|---|
| 40 | return 1; |
|---|
| 41 | |
|---|
| 42 | // +Inf / positive = +Inf |
|---|
| 43 | if (test__divdf3(a: makeInf64(), b: 3., UINT64_C(0x7ff0000000000000))) |
|---|
| 44 | return 1; |
|---|
| 45 | // +Inf / negative = -Inf |
|---|
| 46 | if (test__divdf3(a: makeInf64(), b: -3., UINT64_C(0xfff0000000000000))) |
|---|
| 47 | return 1; |
|---|
| 48 | // -Inf / positive = -Inf |
|---|
| 49 | if (test__divdf3(a: makeNegativeInf64(), b: 3., UINT64_C(0xfff0000000000000))) |
|---|
| 50 | return 1; |
|---|
| 51 | // -Inf / negative = +Inf |
|---|
| 52 | if (test__divdf3(a: makeNegativeInf64(), b: -3., UINT64_C(0x7ff0000000000000))) |
|---|
| 53 | return 1; |
|---|
| 54 | |
|---|
| 55 | // Inf / Inf = NaN |
|---|
| 56 | if (test__divdf3(a: makeInf64(), b: makeInf64(), UINT64_C(0x7ff8000000000000))) |
|---|
| 57 | return 1; |
|---|
| 58 | // 0.0 / 0.0 = NaN |
|---|
| 59 | if (test__divdf3(a: +0x0.0p+0, b: +0x0.0p+0, UINT64_C(0x7ff8000000000000))) |
|---|
| 60 | return 1; |
|---|
| 61 | // +0.0 / +Inf = +0.0 |
|---|
| 62 | if (test__divdf3(a: +0x0.0p+0, b: makeInf64(), UINT64_C(0x0))) |
|---|
| 63 | return 1; |
|---|
| 64 | // +Inf / +0.0 = +Inf |
|---|
| 65 | if (test__divdf3(a: makeInf64(), b: +0x0.0p+0, UINT64_C(0x7ff0000000000000))) |
|---|
| 66 | return 1; |
|---|
| 67 | |
|---|
| 68 | // positive / +0.0 = +Inf |
|---|
| 69 | if (test__divdf3(a: +1.0, b: +0x0.0p+0, UINT64_C(0x7ff0000000000000))) |
|---|
| 70 | return 1; |
|---|
| 71 | // positive / -0.0 = -Inf |
|---|
| 72 | if (test__divdf3(a: +1.0, b: -0x0.0p+0, UINT64_C(0xfff0000000000000))) |
|---|
| 73 | return 1; |
|---|
| 74 | // negative / +0.0 = -Inf |
|---|
| 75 | if (test__divdf3(a: -1.0, b: +0x0.0p+0, UINT64_C(0xfff0000000000000))) |
|---|
| 76 | return 1; |
|---|
| 77 | // negative / -0.0 = +Inf |
|---|
| 78 | if (test__divdf3(a: -1.0, b: -0x0.0p+0, UINT64_C(0x7ff0000000000000))) |
|---|
| 79 | return 1; |
|---|
| 80 | |
|---|
| 81 | // 1/3 |
|---|
| 82 | if (test__divdf3(a: 1., b: 3., UINT64_C(0x3fd5555555555555))) |
|---|
| 83 | return 1; |
|---|
| 84 | // smallest normal result |
|---|
| 85 | if (test__divdf3(a: 0x1.0p-1021, b: 2., UINT64_C(0x10000000000000))) |
|---|
| 86 | return 1; |
|---|
| 87 | |
|---|
| 88 | // divisor is exactly 1.0 |
|---|
| 89 | if (test__divdf3(a: 0x1.0p+0, b: 0x1.0p+0, UINT64_C(0x3ff0000000000000))) |
|---|
| 90 | return 1; |
|---|
| 91 | // divisor is truncated to exactly 1.0 in UQ1.31 |
|---|
| 92 | if (test__divdf3(a: 0x1.0p+0, b: 0x1.00000001p+0, UINT64_C(0x3fefffffffe00000))) |
|---|
| 93 | return 1; |
|---|
| 94 | |
|---|
| 95 | // smallest normal value divided by 2.0 |
|---|
| 96 | if (test__divdf3(a: 0x1.0p-1022, b: 2., UINT64_C(0x0008000000000000))) |
|---|
| 97 | return 1; |
|---|
| 98 | // smallest subnormal result |
|---|
| 99 | if (test__divdf3(a: 0x1.0p-1022, b: 0x1.0p+52, UINT64_C(0x0000000000000001))) |
|---|
| 100 | return 1; |
|---|
| 101 | |
|---|
| 102 | // some misc test cases obtained by fuzzing against h/w implementation |
|---|
| 103 | if (test__divdf3(a: 0x1.fdc239dd64735p-658, b: -0x1.fff9364c0843fp-948, UINT64_C(0xd20fdc8fc0ceffb1))) |
|---|
| 104 | return 1; |
|---|
| 105 | if (test__divdf3(a: -0x1.78abb261d47c8p+794, b: 0x1.fb01d537cc5aep+266, UINT64_C(0xe0e7c6148ffc23e3))) |
|---|
| 106 | return 1; |
|---|
| 107 | if (test__divdf3(a: -0x1.da7dfe6048b8bp-875, b: 0x1.ffc7ea3ff60a4p-610, UINT64_C(0xaf5dab1fe0269e2a))) |
|---|
| 108 | return 1; |
|---|
| 109 | if (test__divdf3(a: 0x1.0p-1022, b: 0x1.9p+5, UINT64_C(0x000051eb851eb852))) |
|---|
| 110 | return 1; |
|---|
| 111 | if (test__divdf3(a: 0x1.0p-1022, b: 0x1.0028p+41, UINT64_C(0x00000000000007ff))) |
|---|
| 112 | return 1; |
|---|
| 113 | if (test__divdf3(a: 0x1.0p-1022, b: 0x1.0028p+52, UINT64_C(0x1))) |
|---|
| 114 | return 1; |
|---|
| 115 | |
|---|
| 116 | // test 1 / (1 - eps(0.5)) = 1 + eps(1) |
|---|
| 117 | if (test__divdf3(a: 1.0, b: 0x1.fffffffffffffp-1, UINT64_C(0x3ff0000000000001))) |
|---|
| 118 | return 1; |
|---|
| 119 | |
|---|
| 120 | return 0; |
|---|
| 121 | } |
|---|
| 122 | |
|---|
