1// RUN: %clang_builtins %s %librt -o %t && %run %t
2// REQUIRES: librt_has_divtf3
3
4#include "int_lib.h"
5#include <stdio.h>
6
7// The testcase currently assumes IEEE TF format, once that has been
8// fixed the defined(CRT_HAS_IEEE_TF) guard can be removed to enable it for
9// IBM 128 floats as well.
10#if defined(CRT_HAS_IEEE_TF)
11
12# include "fp_test.h"
13
14// Returns: a / b
15COMPILER_RT_ABI tf_float __divtf3(tf_float a, tf_float b);
16
17int test__divtf3(tf_float a, tf_float b, uint64_t expectedHi,
18 uint64_t expectedLo) {
19 tf_float x = __divtf3(a, b);
20 int ret = compareResultF128(result: x, expectedHi, expectedLo);
21
22 if (ret) {
23 printf(format: "error in test__divtf3(%.20Le, %.20Le) = %.20Le, "
24 "expected %.20Le\n",
25 a, b, x, fromRep128(hi: expectedHi, lo: expectedLo));
26 }
27 return ret;
28}
29
30char assumption_1[sizeof(tf_float) * CHAR_BIT == 128] = {0};
31
32#endif
33
34int main() {
35#if defined(CRT_HAS_IEEE_TF)
36 // Returned NaNs are assumed to be qNaN by default
37
38 // qNaN / any = qNaN
39 if (test__divtf3(a: makeQNaN128(), TF_C(0x1.23456789abcdefp+5),
40 UINT64_C(0x7fff800000000000), UINT64_C(0x0)))
41 return 1;
42 // NaN / any = NaN
43 if (test__divtf3(a: makeNaN128(UINT64_C(0x30000000)),
44 TF_C(0x1.23456789abcdefp+5), UINT64_C(0x7fff800000000000),
45 UINT64_C(0x0)))
46 return 1;
47 // any / qNaN = qNaN
48 if (test__divtf3(TF_C(0x1.23456789abcdefp+5), b: makeQNaN128(),
49 UINT64_C(0x7fff800000000000), UINT64_C(0x0)))
50 return 1;
51 // any / NaN = NaN
52 if (test__divtf3(TF_C(0x1.23456789abcdefp+5),
53 b: makeNaN128(UINT64_C(0x30000000)),
54 UINT64_C(0x7fff800000000000), UINT64_C(0x0)))
55 return 1;
56
57 // +Inf / positive = +Inf
58 if (test__divtf3(a: makeInf128(), TF_C(3.), UINT64_C(0x7fff000000000000),
59 UINT64_C(0x0)))
60 return 1;
61 // +Inf / negative = -Inf
62 if (test__divtf3(a: makeInf128(), b: -TF_C(3.), UINT64_C(0xffff000000000000),
63 UINT64_C(0x0)))
64 return 1;
65 // -Inf / positive = -Inf
66 if (test__divtf3(a: makeNegativeInf128(), TF_C(3.), UINT64_C(0xffff000000000000),
67 UINT64_C(0x0)))
68 return 1;
69 // -Inf / negative = +Inf
70 if (test__divtf3(a: makeNegativeInf128(), b: -TF_C(3.),
71 UINT64_C(0x7fff000000000000), UINT64_C(0x0)))
72 return 1;
73
74 // Inf / Inf = NaN
75 if (test__divtf3(a: makeInf128(), b: makeInf128(), UINT64_C(0x7fff800000000000),
76 UINT64_C(0x0)))
77 return 1;
78 // 0.0 / 0.0 = NaN
79 if (test__divtf3(a: +TF_C(0x0.0p+0), b: +TF_C(0x0.0p+0),
80 UINT64_C(0x7fff800000000000), UINT64_C(0x0)))
81 return 1;
82 // +0.0 / +Inf = +0.0
83 if (test__divtf3(a: +TF_C(0x0.0p+0), b: makeInf128(), UINT64_C(0x0), UINT64_C(0x0)))
84 return 1;
85 // +Inf / +0.0 = +Inf
86 if (test__divtf3(a: makeInf128(), b: +TF_C(0x0.0p+0), UINT64_C(0x7fff000000000000),
87 UINT64_C(0x0)))
88 return 1;
89
90 // positive / +0.0 = +Inf
91 if (test__divtf3(a: +TF_C(1.0), b: +TF_C(0x0.0p+0), UINT64_C(0x7fff000000000000),
92 UINT64_C(0x0)))
93 return 1;
94 // positive / -0.0 = -Inf
95 if (test__divtf3(a: +1.0L, b: -TF_C(0x0.0p+0), UINT64_C(0xffff000000000000),
96 UINT64_C(0x0)))
97 return 1;
98 // negative / +0.0 = -Inf
99 if (test__divtf3(a: -1.0L, b: +TF_C(0x0.0p+0), UINT64_C(0xffff000000000000),
100 UINT64_C(0x0)))
101 return 1;
102 // negative / -0.0 = +Inf
103 if (test__divtf3(TF_C(-1.0), b: -TF_C(0x0.0p+0), UINT64_C(0x7fff000000000000),
104 UINT64_C(0x0)))
105 return 1;
106
107 // 1/3
108 if (test__divtf3(TF_C(1.), TF_C(3.), UINT64_C(0x3ffd555555555555),
109 UINT64_C(0x5555555555555555)))
110 return 1;
111 // smallest normal result
112 if (test__divtf3(TF_C(0x1.0p-16381), TF_C(2.), UINT64_C(0x0001000000000000),
113 UINT64_C(0x0)))
114 return 1;
115
116 // divisor is exactly 1.0
117 if (test__divtf3(TF_C(0x1.0p+0), TF_C(0x1.0p+0), UINT64_C(0x3fff000000000000),
118 UINT64_C(0x0)))
119 return 1;
120 // divisor is truncated to exactly 1.0 in UQ1.63
121 if (test__divtf3(TF_C(0x1.0p+0), TF_C(0x1.0000000000000001p+0),
122 UINT64_C(0x3ffeffffffffffff), UINT64_C(0xfffe000000000000)))
123 return 1;
124
125 // smallest normal value divided by 2.0
126 if (test__divtf3(TF_C(0x1.0p-16382), b: 2.L, UINT64_C(0x0000800000000000),
127 UINT64_C(0x0)))
128 return 1;
129 // smallest subnormal result
130 if (test__divtf3(TF_C(0x1.0p-16382), TF_C(0x1p+112), UINT64_C(0x0),
131 UINT64_C(0x1)))
132 return 1;
133
134 // any / any
135 if (test__divtf3(TF_C(0x1.a23b45362464523375893ab4cdefp+5),
136 TF_C(0x1.eedcbaba3a94546558237654321fp-1),
137 UINT64_C(0x4004b0b72924d407), UINT64_C(0x0717e84356c6eba2)))
138 return 1;
139 if (test__divtf3(TF_C(0x1.a2b34c56d745382f9abf2c3dfeffp-50),
140 TF_C(0x1.ed2c3ba15935332532287654321fp-9),
141 UINT64_C(0x3fd5b2af3f828c9b), UINT64_C(0x40e51f64cde8b1f2)))
142 return 15;
143 if (test__divtf3(TF_C(0x1.2345f6aaaa786555f42432abcdefp+456),
144 TF_C(0x1.edacbba9874f765463544dd3621fp+6400),
145 UINT64_C(0x28c62e15dc464466), UINT64_C(0xb5a07586348557ac)))
146 return 1;
147 if (test__divtf3(TF_C(0x1.2d3456f789ba6322bc665544edefp-234),
148 TF_C(0x1.eddcdba39f3c8b7a36564354321fp-4455),
149 UINT64_C(0x507b38442b539266), UINT64_C(0x22ce0f1d024e1252)))
150 return 1;
151 if (test__divtf3(TF_C(0x1.2345f6b77b7a8953365433abcdefp+234),
152 TF_C(0x1.edcba987d6bb3aa467754354321fp-4055),
153 UINT64_C(0x50bf2e02f0798d36), UINT64_C(0x5e6fcb6b60044078)))
154 return 1;
155 if (test__divtf3(TF_C(6.72420628622418701252535563464350521E-4932), TF_C(2.),
156 UINT64_C(0x0001000000000000), UINT64_C(0)))
157 return 1;
158
159 // test 1 / (1 - eps(0.5)) = 1 + eps(1).
160 if (test__divtf3(a: 1.0L, TF_C(0x1.ffffffffffffffffffffffffffffp-1),
161 UINT64_C(0x3FFF000000000000), UINT64_C(1)))
162 return 1;
163
164#else
165 printf("skipped\n");
166
167#endif
168 return 0;
169}
170

source code of compiler-rt/test/builtins/Unit/divtf3_test.c