1 |
|
---|
2 | /*============================================================================
|
---|
3 |
|
---|
4 | This C header file is part of the SoftFloat IEEE Floating-Point Arithmetic
|
---|
5 | Package, Release 3e, by John R. Hauser.
|
---|
6 |
|
---|
7 | Copyright 2017 The Regents of the University of California. All rights
|
---|
8 | reserved.
|
---|
9 |
|
---|
10 | Redistribution and use in source and binary forms, with or without
|
---|
11 | modification, are permitted provided that the following conditions are met:
|
---|
12 |
|
---|
13 | 1. Redistributions of source code must retain the above copyright notice,
|
---|
14 | this list of conditions, and the following disclaimer.
|
---|
15 |
|
---|
16 | 2. Redistributions in binary form must reproduce the above copyright notice,
|
---|
17 | this list of conditions, and the following disclaimer in the documentation
|
---|
18 | and/or other materials provided with the distribution.
|
---|
19 |
|
---|
20 | 3. Neither the name of the University nor the names of its contributors may
|
---|
21 | be used to endorse or promote products derived from this software without
|
---|
22 | specific prior written permission.
|
---|
23 |
|
---|
24 | THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS "AS IS", AND ANY
|
---|
25 | EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
|
---|
26 | WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE, ARE
|
---|
27 | DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE FOR ANY
|
---|
28 | DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
|
---|
29 | (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
|
---|
30 | LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
|
---|
31 | ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
|
---|
32 | (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
|
---|
33 | SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
---|
34 |
|
---|
35 | =============================================================================*/
|
---|
36 |
|
---|
37 | #ifndef opts_GCC_h
|
---|
38 | #define opts_GCC_h 1
|
---|
39 |
|
---|
40 | #ifdef INLINE
|
---|
41 |
|
---|
42 | #include <stdint.h>
|
---|
43 | #include "primitiveTypes.h"
|
---|
44 |
|
---|
45 | #ifdef SOFTFLOAT_BUILTIN_CLZ
|
---|
46 |
|
---|
47 | INLINE uint_fast8_t softfloat_countLeadingZeros16( uint16_t a )
|
---|
48 | { return a ? __builtin_clz( a ) - 16 : 16; }
|
---|
49 | #define softfloat_countLeadingZeros16 softfloat_countLeadingZeros16
|
---|
50 |
|
---|
51 | INLINE uint_fast8_t softfloat_countLeadingZeros32( uint32_t a )
|
---|
52 | { return a ? __builtin_clz( a ) : 32; }
|
---|
53 | #define softfloat_countLeadingZeros32 softfloat_countLeadingZeros32
|
---|
54 |
|
---|
55 | INLINE uint_fast8_t softfloat_countLeadingZeros64( uint64_t a )
|
---|
56 | { return a ? __builtin_clzll( a ) : 64; }
|
---|
57 | #define softfloat_countLeadingZeros64 softfloat_countLeadingZeros64
|
---|
58 |
|
---|
59 | #endif
|
---|
60 |
|
---|
61 | #ifdef SOFTFLOAT_INTRINSIC_INT128
|
---|
62 |
|
---|
63 | INLINE struct uint128 softfloat_mul64ByShifted32To128( uint64_t a, uint32_t b )
|
---|
64 | {
|
---|
65 | union { unsigned __int128 ui; struct uint128 s; } uZ;
|
---|
66 | uZ.ui = (unsigned __int128) a * ((uint_fast64_t) b<<32);
|
---|
67 | return uZ.s;
|
---|
68 | }
|
---|
69 | #define softfloat_mul64ByShifted32To128 softfloat_mul64ByShifted32To128
|
---|
70 |
|
---|
71 | INLINE struct uint128 softfloat_mul64To128( uint64_t a, uint64_t b )
|
---|
72 | {
|
---|
73 | union { unsigned __int128 ui; struct uint128 s; } uZ;
|
---|
74 | uZ.ui = (unsigned __int128) a * b;
|
---|
75 | return uZ.s;
|
---|
76 | }
|
---|
77 | #define softfloat_mul64To128 softfloat_mul64To128
|
---|
78 |
|
---|
79 | INLINE
|
---|
80 | struct uint128 softfloat_mul128By32( uint64_t a64, uint64_t a0, uint32_t b )
|
---|
81 | {
|
---|
82 | union { unsigned __int128 ui; struct uint128 s; } uZ;
|
---|
83 | uZ.ui = ((unsigned __int128) a64<<64 | a0) * b;
|
---|
84 | return uZ.s;
|
---|
85 | }
|
---|
86 | #define softfloat_mul128By32 softfloat_mul128By32
|
---|
87 |
|
---|
88 | INLINE
|
---|
89 | void
|
---|
90 | softfloat_mul128To256M(
|
---|
91 | uint64_t a64, uint64_t a0, uint64_t b64, uint64_t b0, uint64_t *zPtr )
|
---|
92 | {
|
---|
93 | unsigned __int128 z0, mid1, mid, z128;
|
---|
94 | z0 = (unsigned __int128) a0 * b0;
|
---|
95 | mid1 = (unsigned __int128) a64 * b0;
|
---|
96 | mid = mid1 + (unsigned __int128) a0 * b64;
|
---|
97 | z128 = (unsigned __int128) a64 * b64;
|
---|
98 | z128 += (unsigned __int128) (mid < mid1)<<64 | mid>>64;
|
---|
99 | mid <<= 64;
|
---|
100 | z0 += mid;
|
---|
101 | z128 += (z0 < mid);
|
---|
102 | zPtr[indexWord( 4, 0 )] = z0;
|
---|
103 | zPtr[indexWord( 4, 1 )] = z0>>64;
|
---|
104 | zPtr[indexWord( 4, 2 )] = z128;
|
---|
105 | zPtr[indexWord( 4, 3 )] = z128>>64;
|
---|
106 | }
|
---|
107 | #define softfloat_mul128To256M softfloat_mul128To256M
|
---|
108 |
|
---|
109 | #endif
|
---|
110 |
|
---|
111 | #endif
|
---|
112 |
|
---|
113 | #endif
|
---|
114 |
|
---|