• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 /* ===-- mulsc3.c - Implement __mulsc3 -------------------------------------===
2  *
3  *                     The LLVM Compiler Infrastructure
4  *
5  * This file is dual licensed under the MIT and the University of Illinois Open
6  * Source Licenses. See LICENSE.TXT for details.
7  *
8  * ===----------------------------------------------------------------------===
9  *
10  * This file implements __mulsc3 for the compiler_rt library.
11  *
12  * ===----------------------------------------------------------------------===
13  */
14 
15 #include "int_lib.h"
16 #include "int_math.h"
17 
18 /* Returns: the product of a + ib and c + id */
19 
20 COMPILER_RT_ABI Fcomplex
__mulsc3(float __a,float __b,float __c,float __d)21 __mulsc3(float __a, float __b, float __c, float __d)
22 {
23     float __ac = __a * __c;
24     float __bd = __b * __d;
25     float __ad = __a * __d;
26     float __bc = __b * __c;
27     Fcomplex z;
28     COMPLEX_REAL(z) = __ac - __bd;
29     COMPLEX_IMAGINARY(z) = __ad + __bc;
30     if (crt_isnan(COMPLEX_REAL(z)) && crt_isnan(COMPLEX_IMAGINARY(z)))
31     {
32         int __recalc = 0;
33         if (crt_isinf(__a) || crt_isinf(__b))
34         {
35             __a = crt_copysignf(crt_isinf(__a) ? 1 : 0, __a);
36             __b = crt_copysignf(crt_isinf(__b) ? 1 : 0, __b);
37             if (crt_isnan(__c))
38                 __c = crt_copysignf(0, __c);
39             if (crt_isnan(__d))
40                 __d = crt_copysignf(0, __d);
41             __recalc = 1;
42         }
43         if (crt_isinf(__c) || crt_isinf(__d))
44         {
45             __c = crt_copysignf(crt_isinf(__c) ? 1 : 0, __c);
46             __d = crt_copysignf(crt_isinf(__d) ? 1 : 0, __d);
47             if (crt_isnan(__a))
48                 __a = crt_copysignf(0, __a);
49             if (crt_isnan(__b))
50                 __b = crt_copysignf(0, __b);
51             __recalc = 1;
52         }
53         if (!__recalc && (crt_isinf(__ac) || crt_isinf(__bd) ||
54                           crt_isinf(__ad) || crt_isinf(__bc)))
55         {
56             if (crt_isnan(__a))
57                 __a = crt_copysignf(0, __a);
58             if (crt_isnan(__b))
59                 __b = crt_copysignf(0, __b);
60             if (crt_isnan(__c))
61                 __c = crt_copysignf(0, __c);
62             if (crt_isnan(__d))
63                 __d = crt_copysignf(0, __d);
64             __recalc = 1;
65         }
66         if (__recalc)
67         {
68             COMPLEX_REAL(z) = CRT_INFINITY * (__a * __c - __b * __d);
69             COMPLEX_IMAGINARY(z) = CRT_INFINITY * (__a * __d + __b * __c);
70         }
71     }
72     return z;
73 }
74