ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 1 | /* |
| 2 | * Utility compute operations used by translated code. |
| 3 | * |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 4 | * Copyright (c) 2003 Fabrice Bellard |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 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 | |
Blue Swirl | facd285 | 2009-08-16 08:03:26 +0000 | [diff] [blame] | 26 | #include <stdlib.h> |
| 27 | #include <stdint.h> |
j_mayer | 084ed5c | 2007-11-05 13:01:41 +0000 | [diff] [blame] | 28 | #include "host-utils.h" |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 29 | |
ths | 5592a75 | 2007-10-26 22:35:02 +0000 | [diff] [blame] | 30 | //#define DEBUG_MULDIV |
| 31 | |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 32 | /* Long integer helpers */ |
j_mayer | 7a51ad8 | 2007-11-04 02:24:58 +0000 | [diff] [blame] | 33 | #if !defined(__x86_64__) |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 34 | static void add128 (uint64_t *plow, uint64_t *phigh, uint64_t a, uint64_t b) |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 35 | { |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 36 | *plow += a; |
| 37 | /* carry test */ |
| 38 | if (*plow < a) |
| 39 | (*phigh)++; |
| 40 | *phigh += b; |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 41 | } |
| 42 | |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 43 | static void neg128 (uint64_t *plow, uint64_t *phigh) |
| 44 | { |
| 45 | *plow = ~*plow; |
| 46 | *phigh = ~*phigh; |
| 47 | add128(plow, phigh, 1, 0); |
| 48 | } |
| 49 | |
| 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; |
| 54 | |
| 55 | a0 = a; |
| 56 | a1 = a >> 32; |
| 57 | |
| 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; |
| 73 | } |
| 74 | |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 75 | /* Unsigned 64x64 -> 128 multiplication */ |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 76 | void mulu64 (uint64_t *plow, uint64_t *phigh, uint64_t a, uint64_t b) |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 77 | { |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 78 | mul64(plow, phigh, a, b); |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 79 | #if defined(DEBUG_MULDIV) |
| 80 | printf("mulu64: 0x%016llx * 0x%016llx = 0x%016llx%016llx\n", |
| 81 | a, b, *phigh, *plow); |
| 82 | #endif |
| 83 | } |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 84 | |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 85 | /* Signed 64x64 -> 128 multiplication */ |
| 86 | void muls64 (uint64_t *plow, uint64_t *phigh, int64_t a, int64_t b) |
| 87 | { |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 88 | int sa, sb; |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 89 | |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 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 | } |
ths | e494ead | 2007-10-25 23:00:03 +0000 | [diff] [blame] | 100 | #if defined(DEBUG_MULDIV) |
| 101 | printf("muls64: 0x%016llx * 0x%016llx = 0x%016llx%016llx\n", |
| 102 | a, b, *phigh, *plow); |
ths | 69d3572 | 2007-05-16 11:59:40 +0000 | [diff] [blame] | 103 | #endif |
| 104 | } |
j_mayer | 7a51ad8 | 2007-11-04 02:24:58 +0000 | [diff] [blame] | 105 | #endif /* !defined(__x86_64__) */ |