mirror of
https://github.com/llvm/llvm-project.git
synced 2025-04-26 18:26:05 +00:00
99c457dc2e
This re-lands cc0065a7d082f0bd322a538cf62cfaef1c8f89f8 in a way that keeps existing targets working. --------- Original commit message: #68132 ended up removing __multc3 & __divtc3 from compiler-rt library builds that have QUAD_PRECISION but not TF_MODE due to missing int128 support. I added support for QUAD_PRECISION to use the native hex float long double representation. --------- Co-authored-by: Sean Perry <perry@ca.ibm.com>
71 lines
2.0 KiB
C
71 lines
2.0 KiB
C
//===-- multc3.c - Implement __multc3 -------------------------------------===//
|
|
//
|
|
// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
|
|
// See https://llvm.org/LICENSE.txt for license information.
|
|
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
|
|
//
|
|
//===----------------------------------------------------------------------===//
|
|
//
|
|
// This file implements __multc3 for the compiler_rt library.
|
|
//
|
|
//===----------------------------------------------------------------------===//
|
|
|
|
#define QUAD_PRECISION
|
|
#include "fp_lib.h"
|
|
#include "int_lib.h"
|
|
#include "int_math.h"
|
|
|
|
#if defined(CRT_HAS_F128)
|
|
|
|
// Returns: the product of a + ib and c + id
|
|
|
|
COMPILER_RT_ABI Qcomplex __multc3(fp_t a, fp_t b, fp_t c, fp_t d) {
|
|
fp_t ac = a * c;
|
|
fp_t bd = b * d;
|
|
fp_t ad = a * d;
|
|
fp_t bc = b * c;
|
|
Qcomplex z;
|
|
COMPLEXTF_REAL(z) = ac - bd;
|
|
COMPLEXTF_IMAGINARY(z) = ad + bc;
|
|
if (crt_isnan(COMPLEXTF_REAL(z)) && crt_isnan(COMPLEXTF_IMAGINARY(z))) {
|
|
int recalc = 0;
|
|
if (crt_isinf(a) || crt_isinf(b)) {
|
|
a = crt_copysigntf(crt_isinf(a) ? 1 : 0, a);
|
|
b = crt_copysigntf(crt_isinf(b) ? 1 : 0, b);
|
|
if (crt_isnan(c))
|
|
c = crt_copysigntf(0, c);
|
|
if (crt_isnan(d))
|
|
d = crt_copysigntf(0, d);
|
|
recalc = 1;
|
|
}
|
|
if (crt_isinf(c) || crt_isinf(d)) {
|
|
c = crt_copysigntf(crt_isinf(c) ? 1 : 0, c);
|
|
d = crt_copysigntf(crt_isinf(d) ? 1 : 0, d);
|
|
if (crt_isnan(a))
|
|
a = crt_copysigntf(0, a);
|
|
if (crt_isnan(b))
|
|
b = crt_copysigntf(0, b);
|
|
recalc = 1;
|
|
}
|
|
if (!recalc &&
|
|
(crt_isinf(ac) || crt_isinf(bd) || crt_isinf(ad) || crt_isinf(bc))) {
|
|
if (crt_isnan(a))
|
|
a = crt_copysigntf(0, a);
|
|
if (crt_isnan(b))
|
|
b = crt_copysigntf(0, b);
|
|
if (crt_isnan(c))
|
|
c = crt_copysigntf(0, c);
|
|
if (crt_isnan(d))
|
|
d = crt_copysigntf(0, d);
|
|
recalc = 1;
|
|
}
|
|
if (recalc) {
|
|
COMPLEXTF_REAL(z) = CRT_INFINITY * (a * c - b * d);
|
|
COMPLEXTF_IMAGINARY(z) = CRT_INFINITY * (a * d + b * c);
|
|
}
|
|
}
|
|
return z;
|
|
}
|
|
|
|
#endif
|