1
0
Fork 0
forked from suyu/suyu
suyu/src/common/uint128.cpp

42 lines
1 KiB
C++
Raw Normal View History

2019-02-16 01:04:11 +01:00
#ifdef _MSC_VER
#include <intrin.h>
#pragma intrinsic(_umul128)
#endif
2019-02-16 21:52:24 +01:00
#include <cstring>
2019-02-16 01:04:11 +01:00
#include "common/uint128.h"
namespace Common {
2019-02-16 21:52:24 +01:00
2019-02-16 01:04:11 +01:00
u128 Multiply64Into128(u64 a, u64 b) {
u128 result;
2019-02-16 21:52:24 +01:00
#ifdef _MSC_VER
2019-02-16 01:04:11 +01:00
result[0] = _umul128(a, b, &result[1]);
#else
unsigned __int128 tmp = a;
tmp *= b;
std::memcpy(&result, &tmp, sizeof(u128));
#endif
return result;
}
2019-02-16 21:52:24 +01:00
std::pair<u64, u64> Divide128On32(u128 dividend, u32 divisor) {
u64 remainder = dividend[0] % divisor;
u64 accum = dividend[0] / divisor;
if (dividend[1] == 0)
return {accum, remainder};
// We ignore dividend[1] / divisor as that overflows
2019-02-16 21:52:24 +01:00
const u64 first_segment = (dividend[1] % divisor) << 32;
accum += (first_segment / divisor) << 32;
2019-02-16 21:52:24 +01:00
const u64 second_segment = (first_segment % divisor) << 32;
accum += (second_segment / divisor);
remainder += second_segment % divisor;
2019-02-16 01:04:11 +01:00
if (remainder >= divisor) {
accum++;
remainder -= divisor;
}
return {accum, remainder};
}
} // namespace Common