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