10b57cec5SDimitry Andric // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
20b57cec5SDimitry Andric // See https://llvm.org/LICENSE.txt for license information.
30b57cec5SDimitry Andric // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
40b57cec5SDimitry Andric
50b57cec5SDimitry Andric #include "../int_math.h"
60b57cec5SDimitry Andric #include "DD.h"
70b57cec5SDimitry Andric
80b57cec5SDimitry Andric #define makeFinite(x) \
90b57cec5SDimitry Andric { \
100b57cec5SDimitry Andric (x).s.hi = crt_copysign(crt_isinf((x).s.hi) ? 1.0 : 0.0, (x).s.hi); \
110b57cec5SDimitry Andric (x).s.lo = 0.0; \
120b57cec5SDimitry Andric }
130b57cec5SDimitry Andric
140b57cec5SDimitry Andric #define zeroNaN(x) \
150b57cec5SDimitry Andric { \
160b57cec5SDimitry Andric if (crt_isnan((x).s.hi)) { \
170b57cec5SDimitry Andric (x).s.hi = crt_copysign(0.0, (x).s.hi); \
180b57cec5SDimitry Andric (x).s.lo = 0.0; \
190b57cec5SDimitry Andric } \
200b57cec5SDimitry Andric }
210b57cec5SDimitry Andric
__multc3(long double a,long double b,long double c,long double d)220b57cec5SDimitry Andric long double _Complex __multc3(long double a, long double b, long double c,
230b57cec5SDimitry Andric long double d) {
240b57cec5SDimitry Andric long double ac = __gcc_qmul(a, c);
250b57cec5SDimitry Andric long double bd = __gcc_qmul(b, d);
260b57cec5SDimitry Andric long double ad = __gcc_qmul(a, d);
270b57cec5SDimitry Andric long double bc = __gcc_qmul(b, c);
280b57cec5SDimitry Andric
290b57cec5SDimitry Andric DD real = {.ld = __gcc_qsub(ac, bd)};
300b57cec5SDimitry Andric DD imag = {.ld = __gcc_qadd(ad, bc)};
310b57cec5SDimitry Andric
320b57cec5SDimitry Andric if (crt_isnan(real.s.hi) && crt_isnan(imag.s.hi)) {
330b57cec5SDimitry Andric int recalc = 0;
340b57cec5SDimitry Andric
350b57cec5SDimitry Andric DD aDD = {.ld = a};
360b57cec5SDimitry Andric DD bDD = {.ld = b};
370b57cec5SDimitry Andric DD cDD = {.ld = c};
380b57cec5SDimitry Andric DD dDD = {.ld = d};
390b57cec5SDimitry Andric
400b57cec5SDimitry Andric if (crt_isinf(aDD.s.hi) || crt_isinf(bDD.s.hi)) {
410b57cec5SDimitry Andric makeFinite(aDD);
420b57cec5SDimitry Andric makeFinite(bDD);
430b57cec5SDimitry Andric zeroNaN(cDD);
440b57cec5SDimitry Andric zeroNaN(dDD);
450b57cec5SDimitry Andric recalc = 1;
460b57cec5SDimitry Andric }
470b57cec5SDimitry Andric
480b57cec5SDimitry Andric if (crt_isinf(cDD.s.hi) || crt_isinf(dDD.s.hi)) {
490b57cec5SDimitry Andric makeFinite(cDD);
500b57cec5SDimitry Andric makeFinite(dDD);
510b57cec5SDimitry Andric zeroNaN(aDD);
520b57cec5SDimitry Andric zeroNaN(bDD);
530b57cec5SDimitry Andric recalc = 1;
540b57cec5SDimitry Andric }
550b57cec5SDimitry Andric
560b57cec5SDimitry Andric if (!recalc) {
570b57cec5SDimitry Andric DD acDD = {.ld = ac};
580b57cec5SDimitry Andric DD bdDD = {.ld = bd};
590b57cec5SDimitry Andric DD adDD = {.ld = ad};
600b57cec5SDimitry Andric DD bcDD = {.ld = bc};
610b57cec5SDimitry Andric
620b57cec5SDimitry Andric if (crt_isinf(acDD.s.hi) || crt_isinf(bdDD.s.hi) ||
630b57cec5SDimitry Andric crt_isinf(adDD.s.hi) || crt_isinf(bcDD.s.hi)) {
640b57cec5SDimitry Andric zeroNaN(aDD);
650b57cec5SDimitry Andric zeroNaN(bDD);
660b57cec5SDimitry Andric zeroNaN(cDD);
670b57cec5SDimitry Andric zeroNaN(dDD);
680b57cec5SDimitry Andric recalc = 1;
690b57cec5SDimitry Andric }
700b57cec5SDimitry Andric }
710b57cec5SDimitry Andric
720b57cec5SDimitry Andric if (recalc) {
730b57cec5SDimitry Andric real.s.hi = CRT_INFINITY * (aDD.s.hi * cDD.s.hi - bDD.s.hi * dDD.s.hi);
740b57cec5SDimitry Andric real.s.lo = 0.0;
750b57cec5SDimitry Andric imag.s.hi = CRT_INFINITY * (aDD.s.hi * dDD.s.hi + bDD.s.hi * cDD.s.hi);
760b57cec5SDimitry Andric imag.s.lo = 0.0;
770b57cec5SDimitry Andric }
780b57cec5SDimitry Andric }
790b57cec5SDimitry Andric
800b57cec5SDimitry Andric long double _Complex z;
810b57cec5SDimitry Andric __real__ z = real.ld;
820b57cec5SDimitry Andric __imag__ z = imag.ld;
830b57cec5SDimitry Andric
840b57cec5SDimitry Andric return z;
850b57cec5SDimitry Andric }
86