root / host-utils.c @ 67fc07d3
History | View | Annotate | Download (2.5 kB)
1 |
/*
|
---|---|
2 |
* Utility compute operations used by translated code.
|
3 |
*
|
4 |
* Copyright (c) 2007 Aurelien Jarno
|
5 |
*
|
6 |
* Permission is hereby granted, free of charge, to any person obtaining a copy
|
7 |
* of this software and associated documentation files (the "Software"), to deal
|
8 |
* in the Software without restriction, including without limitation the rights
|
9 |
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
10 |
* copies of the Software, and to permit persons to whom the Software is
|
11 |
* furnished to do so, subject to the following conditions:
|
12 |
*
|
13 |
* The above copyright notice and this permission notice shall be included in
|
14 |
* all copies or substantial portions of the Software.
|
15 |
*
|
16 |
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
17 |
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
18 |
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
|
19 |
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
20 |
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
21 |
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
|
22 |
* THE SOFTWARE.
|
23 |
*/
|
24 |
|
25 |
#include "vl.h" |
26 |
|
27 |
/* Signed 64x64 -> 128 multiplication */
|
28 |
|
29 |
void muls64(int64_t *phigh, int64_t *plow, int64_t a, int64_t b)
|
30 |
{ |
31 |
#if defined(__x86_64__)
|
32 |
__asm__ ("imul %0\n\t"
|
33 |
: "=d" (*phigh), "=a" (*plow) |
34 |
: "a" (a), "0" (b) |
35 |
); |
36 |
#else
|
37 |
int64_t ph; |
38 |
uint64_t pm1, pm2, pl; |
39 |
|
40 |
pl = (uint64_t)((uint32_t)a) * (uint64_t)((uint32_t)b); |
41 |
pm1 = (a >> 32) * (uint32_t)b;
|
42 |
pm2 = (uint32_t)a * (b >> 32);
|
43 |
ph = (a >> 32) * (b >> 32); |
44 |
|
45 |
ph += (int64_t)pm1 >> 32;
|
46 |
ph += (int64_t)pm2 >> 32;
|
47 |
pm1 = (uint64_t)((uint32_t)pm1) + (uint64_t)((uint32_t)pm2) + (pl >> 32);
|
48 |
|
49 |
*phigh = ph + ((int64_t)pm1 >> 32);
|
50 |
*plow = (pm1 << 32) + (uint32_t)pl;
|
51 |
#endif
|
52 |
} |
53 |
|
54 |
/* Unsigned 64x64 -> 128 multiplication */
|
55 |
void mulu64(uint64_t *phigh, uint64_t *plow, uint64_t a, uint64_t b)
|
56 |
{ |
57 |
#if defined(__x86_64__)
|
58 |
__asm__ ("mul %0\n\t"
|
59 |
: "=d" (*phigh), "=a" (*plow) |
60 |
: "a" (a), "0" (b) |
61 |
); |
62 |
#else
|
63 |
uint64_t ph, pm1, pm2, pl; |
64 |
|
65 |
pl = (uint64_t)((uint32_t)a) * (uint64_t)((uint32_t)b); |
66 |
pm1 = (a >> 32) * (uint32_t)b;
|
67 |
pm2 = (uint32_t)a * (b >> 32);
|
68 |
ph = (a >> 32) * (b >> 32); |
69 |
|
70 |
ph += pm1 >> 32;
|
71 |
ph += pm2 >> 32;
|
72 |
pm1 = (uint64_t)((uint32_t)pm1) + (uint64_t)((uint32_t)pm2) + (pl >> 32);
|
73 |
|
74 |
*phigh = ph + (pm1 >> 32);
|
75 |
*plow = (pm1 << 32) + (uint32_t)pl;
|
76 |
#endif
|
77 |
} |