]>
Commit | Line | Data |
---|---|---|
69d35728 TS |
1 | /* |
2 | * Utility compute operations used by translated code. | |
3 | * | |
e494ead5 | 4 | * Copyright (c) 2003 Fabrice Bellard |
69d35728 TS |
5 | * Copyright (c) 2007 Aurelien Jarno |
6 | * | |
7 | * Permission is hereby granted, free of charge, to any person obtaining a copy | |
8 | * of this software and associated documentation files (the "Software"), to deal | |
9 | * in the Software without restriction, including without limitation the rights | |
10 | * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell | |
11 | * copies of the Software, and to permit persons to whom the Software is | |
12 | * furnished to do so, subject to the following conditions: | |
13 | * | |
14 | * The above copyright notice and this permission notice shall be included in | |
15 | * all copies or substantial portions of the Software. | |
16 | * | |
17 | * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR | |
18 | * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, | |
19 | * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL | |
20 | * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER | |
21 | * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, | |
22 | * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN | |
23 | * THE SOFTWARE. | |
24 | */ | |
25 | ||
facd2857 BS |
26 | #include <stdlib.h> |
27 | #include <stdint.h> | |
084ed5cc | 28 | #include "host-utils.h" |
69d35728 | 29 | |
5592a750 TS |
30 | //#define DEBUG_MULDIV |
31 | ||
e494ead5 | 32 | /* Long integer helpers */ |
7a51ad82 | 33 | #if !defined(__x86_64__) |
e494ead5 TS |
34 | static void add128 (uint64_t *plow, uint64_t *phigh, uint64_t a, uint64_t b) |
35 | { | |
36 | *plow += a; | |
37 | /* carry test */ | |
38 | if (*plow < a) | |
39 | (*phigh)++; | |
40 | *phigh += b; | |
41 | } | |
69d35728 | 42 | |
e494ead5 | 43 | static void neg128 (uint64_t *plow, uint64_t *phigh) |
69d35728 | 44 | { |
e494ead5 TS |
45 | *plow = ~*plow; |
46 | *phigh = ~*phigh; | |
47 | add128(plow, phigh, 1, 0); | |
48 | } | |
69d35728 | 49 | |
e494ead5 TS |
50 | static void mul64 (uint64_t *plow, uint64_t *phigh, uint64_t a, uint64_t b) |
51 | { | |
52 | uint32_t a0, a1, b0, b1; | |
53 | uint64_t v; | |
69d35728 | 54 | |
e494ead5 TS |
55 | a0 = a; |
56 | a1 = a >> 32; | |
69d35728 | 57 | |
e494ead5 TS |
58 | b0 = b; |
59 | b1 = b >> 32; | |
60 | ||
61 | v = (uint64_t)a0 * (uint64_t)b0; | |
62 | *plow = v; | |
63 | *phigh = 0; | |
64 | ||
65 | v = (uint64_t)a0 * (uint64_t)b1; | |
66 | add128(plow, phigh, v << 32, v >> 32); | |
67 | ||
68 | v = (uint64_t)a1 * (uint64_t)b0; | |
69 | add128(plow, phigh, v << 32, v >> 32); | |
70 | ||
71 | v = (uint64_t)a1 * (uint64_t)b1; | |
72 | *phigh += v; | |
69d35728 TS |
73 | } |
74 | ||
75 | /* Unsigned 64x64 -> 128 multiplication */ | |
e494ead5 | 76 | void mulu64 (uint64_t *plow, uint64_t *phigh, uint64_t a, uint64_t b) |
69d35728 | 77 | { |
e494ead5 | 78 | mul64(plow, phigh, a, b); |
e494ead5 TS |
79 | #if defined(DEBUG_MULDIV) |
80 | printf("mulu64: 0x%016llx * 0x%016llx = 0x%016llx%016llx\n", | |
81 | a, b, *phigh, *plow); | |
82 | #endif | |
83 | } | |
69d35728 | 84 | |
e494ead5 TS |
85 | /* Signed 64x64 -> 128 multiplication */ |
86 | void muls64 (uint64_t *plow, uint64_t *phigh, int64_t a, int64_t b) | |
87 | { | |
e494ead5 | 88 | int sa, sb; |
69d35728 | 89 | |
e494ead5 TS |
90 | sa = (a < 0); |
91 | if (sa) | |
92 | a = -a; | |
93 | sb = (b < 0); | |
94 | if (sb) | |
95 | b = -b; | |
96 | mul64(plow, phigh, a, b); | |
97 | if (sa ^ sb) { | |
98 | neg128(plow, phigh); | |
99 | } | |
e494ead5 TS |
100 | #if defined(DEBUG_MULDIV) |
101 | printf("muls64: 0x%016llx * 0x%016llx = 0x%016llx%016llx\n", | |
102 | a, b, *phigh, *plow); | |
69d35728 TS |
103 | #endif |
104 | } | |
7a51ad82 | 105 | #endif /* !defined(__x86_64__) */ |