1*3cab2bb3Spatrick // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
2*3cab2bb3Spatrick // See https://llvm.org/LICENSE.txt for license information.
3*3cab2bb3Spatrick // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
4*3cab2bb3Spatrick
5*3cab2bb3Spatrick #include "../int_math.h"
6*3cab2bb3Spatrick #include "DD.h"
7*3cab2bb3Spatrick
8*3cab2bb3Spatrick #define makeFinite(x) \
9*3cab2bb3Spatrick { \
10*3cab2bb3Spatrick (x).s.hi = crt_copysign(crt_isinf((x).s.hi) ? 1.0 : 0.0, (x).s.hi); \
11*3cab2bb3Spatrick (x).s.lo = 0.0; \
12*3cab2bb3Spatrick }
13*3cab2bb3Spatrick
14*3cab2bb3Spatrick #define zeroNaN(x) \
15*3cab2bb3Spatrick { \
16*3cab2bb3Spatrick if (crt_isnan((x).s.hi)) { \
17*3cab2bb3Spatrick (x).s.hi = crt_copysign(0.0, (x).s.hi); \
18*3cab2bb3Spatrick (x).s.lo = 0.0; \
19*3cab2bb3Spatrick } \
20*3cab2bb3Spatrick }
21*3cab2bb3Spatrick
__multc3(long double a,long double b,long double c,long double d)22*3cab2bb3Spatrick long double _Complex __multc3(long double a, long double b, long double c,
23*3cab2bb3Spatrick long double d) {
24*3cab2bb3Spatrick long double ac = __gcc_qmul(a, c);
25*3cab2bb3Spatrick long double bd = __gcc_qmul(b, d);
26*3cab2bb3Spatrick long double ad = __gcc_qmul(a, d);
27*3cab2bb3Spatrick long double bc = __gcc_qmul(b, c);
28*3cab2bb3Spatrick
29*3cab2bb3Spatrick DD real = {.ld = __gcc_qsub(ac, bd)};
30*3cab2bb3Spatrick DD imag = {.ld = __gcc_qadd(ad, bc)};
31*3cab2bb3Spatrick
32*3cab2bb3Spatrick if (crt_isnan(real.s.hi) && crt_isnan(imag.s.hi)) {
33*3cab2bb3Spatrick int recalc = 0;
34*3cab2bb3Spatrick
35*3cab2bb3Spatrick DD aDD = {.ld = a};
36*3cab2bb3Spatrick DD bDD = {.ld = b};
37*3cab2bb3Spatrick DD cDD = {.ld = c};
38*3cab2bb3Spatrick DD dDD = {.ld = d};
39*3cab2bb3Spatrick
40*3cab2bb3Spatrick if (crt_isinf(aDD.s.hi) || crt_isinf(bDD.s.hi)) {
41*3cab2bb3Spatrick makeFinite(aDD);
42*3cab2bb3Spatrick makeFinite(bDD);
43*3cab2bb3Spatrick zeroNaN(cDD);
44*3cab2bb3Spatrick zeroNaN(dDD);
45*3cab2bb3Spatrick recalc = 1;
46*3cab2bb3Spatrick }
47*3cab2bb3Spatrick
48*3cab2bb3Spatrick if (crt_isinf(cDD.s.hi) || crt_isinf(dDD.s.hi)) {
49*3cab2bb3Spatrick makeFinite(cDD);
50*3cab2bb3Spatrick makeFinite(dDD);
51*3cab2bb3Spatrick zeroNaN(aDD);
52*3cab2bb3Spatrick zeroNaN(bDD);
53*3cab2bb3Spatrick recalc = 1;
54*3cab2bb3Spatrick }
55*3cab2bb3Spatrick
56*3cab2bb3Spatrick if (!recalc) {
57*3cab2bb3Spatrick DD acDD = {.ld = ac};
58*3cab2bb3Spatrick DD bdDD = {.ld = bd};
59*3cab2bb3Spatrick DD adDD = {.ld = ad};
60*3cab2bb3Spatrick DD bcDD = {.ld = bc};
61*3cab2bb3Spatrick
62*3cab2bb3Spatrick if (crt_isinf(acDD.s.hi) || crt_isinf(bdDD.s.hi) ||
63*3cab2bb3Spatrick crt_isinf(adDD.s.hi) || crt_isinf(bcDD.s.hi)) {
64*3cab2bb3Spatrick zeroNaN(aDD);
65*3cab2bb3Spatrick zeroNaN(bDD);
66*3cab2bb3Spatrick zeroNaN(cDD);
67*3cab2bb3Spatrick zeroNaN(dDD);
68*3cab2bb3Spatrick recalc = 1;
69*3cab2bb3Spatrick }
70*3cab2bb3Spatrick }
71*3cab2bb3Spatrick
72*3cab2bb3Spatrick if (recalc) {
73*3cab2bb3Spatrick real.s.hi = CRT_INFINITY * (aDD.s.hi * cDD.s.hi - bDD.s.hi * dDD.s.hi);
74*3cab2bb3Spatrick real.s.lo = 0.0;
75*3cab2bb3Spatrick imag.s.hi = CRT_INFINITY * (aDD.s.hi * dDD.s.hi + bDD.s.hi * cDD.s.hi);
76*3cab2bb3Spatrick imag.s.lo = 0.0;
77*3cab2bb3Spatrick }
78*3cab2bb3Spatrick }
79*3cab2bb3Spatrick
80*3cab2bb3Spatrick long double _Complex z;
81*3cab2bb3Spatrick __real__ z = real.ld;
82*3cab2bb3Spatrick __imag__ z = imag.ld;
83*3cab2bb3Spatrick
84*3cab2bb3Spatrick return z;
85*3cab2bb3Spatrick }
86