Sweden-Number/dlls/ntdll/large_int.c

980 lines
29 KiB
C

/*
* Large integer functions
*
* Copyright 2000 Alexandre Julliard
* Copyright 2003 Thomas Mertes
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this library; if not, write to the Free Software
* Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
*/
#include <stdarg.h>
#include "ntstatus.h"
#define WIN32_NO_STATUS
#include "windef.h"
#include "winternl.h"
#include "wine/asm.h"
#ifndef _WIN64
/*
* Note: we use LONGLONG instead of LARGE_INTEGER, because
* the latter is a structure and the calling convention for
* returning a structure would not be binary-compatible.
*
* FIXME: for platforms that don't have a native LONGLONG type,
* we should define LONGLONG as a structure similar to LARGE_INTEGER
* and do everything by hand. You are welcome to do it...
*/
/******************************************************************************
* RtlLargeIntegerAdd (NTDLL.@)
*
* Add two 64 bit integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to add to a.
*
* RETURNS
* The sum of a and b.
*/
LONGLONG WINAPI RtlLargeIntegerAdd( LONGLONG a, LONGLONG b )
{
return a + b;
}
/******************************************************************************
* RtlLargeIntegerSubtract (NTDLL.@)
*
* Subtract two 64 bit integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to subtract from a.
*
* RETURNS
* The difference of a and b.
*/
LONGLONG WINAPI RtlLargeIntegerSubtract( LONGLONG a, LONGLONG b )
{
return a - b;
}
/******************************************************************************
* RtlLargeIntegerNegate (NTDLL.@)
*
* Negate a 64 bit integer.
*
* PARAMS
* a [I] Initial number.
*
* RETURNS
* The value of a negated.
*/
LONGLONG WINAPI RtlLargeIntegerNegate( LONGLONG a )
{
return -a;
}
/******************************************************************************
* RtlLargeIntegerShiftLeft (NTDLL.@)
*
* Perform a shift left on a 64 bit integer.
*
* PARAMS
* a [I] Initial number.
* count [I] Number of bits to shift by
*
* RETURNS
* The value of a following the shift.
*/
LONGLONG WINAPI RtlLargeIntegerShiftLeft( LONGLONG a, INT count )
{
return a << count;
}
/******************************************************************************
* RtlLargeIntegerShiftRight (NTDLL.@)
*
* Perform a shift right on a 64 bit integer.
*
* PARAMS
* a [I] Initial number.
* count [I] Number of bits to shift by
*
* RETURNS
* The value of a following the shift.
*/
LONGLONG WINAPI RtlLargeIntegerShiftRight( LONGLONG a, INT count )
{
return (ULONGLONG)a >> count;
}
/******************************************************************************
* RtlLargeIntegerArithmeticShift (NTDLL.@)
*
* Perform an arithmetic shift right on a 64 bit integer.
*
* PARAMS
* a [I] Initial number.
* count [I] Number of bits to shift by
*
* RETURNS
* The value of a following the shift.
*/
LONGLONG WINAPI RtlLargeIntegerArithmeticShift( LONGLONG a, INT count )
{
/* FIXME: gcc does arithmetic shift here, but it may not be true on all platforms */
return a >> count;
}
/******************************************************************************
* RtlLargeIntegerDivide (NTDLL.@)
*
* Divide one 64 bit unsigned integer by another, with remainder.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by
* rem [O] Destination for remainder
*
* RETURNS
* The dividend of a and b. If rem is non-NULL it is set to the remainder.
*
* FIXME
* Should it be signed division instead?
*/
ULONGLONG WINAPI RtlLargeIntegerDivide( ULONGLONG a, ULONGLONG b, ULONGLONG *rem )
{
ULONGLONG ret = a / b;
if (rem) *rem = a - ret * b;
return ret;
}
/******************************************************************************
* RtlConvertLongToLargeInteger (NTDLL.@)
*
* Convert a 32 bit integer into 64 bits.
*
* PARAMS
* a [I] Number to convert
*
* RETURNS
* a.
*/
LONGLONG WINAPI RtlConvertLongToLargeInteger( LONG a )
{
return a;
}
/******************************************************************************
* RtlConvertUlongToLargeInteger (NTDLL.@)
*
* Convert a 32 bit unsigned integer into 64 bits.
*
* PARAMS
* a [I] Number to convert
*
* RETURNS
* a.
*/
ULONGLONG WINAPI RtlConvertUlongToLargeInteger( ULONG a )
{
return a;
}
/******************************************************************************
* RtlEnlargedIntegerMultiply (NTDLL.@)
*
* Multiply two integers giving a 64 bit integer result.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to multiply a by.
*
* RETURNS
* The product of a and b.
*/
LONGLONG WINAPI RtlEnlargedIntegerMultiply( INT a, INT b )
{
return (LONGLONG)a * b;
}
/******************************************************************************
* RtlEnlargedUnsignedMultiply (NTDLL.@)
*
* Multiply two unsigned integers giving a 64 bit unsigned integer result.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to multiply a by.
*
* RETURNS
* The product of a and b.
*/
ULONGLONG WINAPI RtlEnlargedUnsignedMultiply( UINT a, UINT b )
{
return (ULONGLONG)a * b;
}
/******************************************************************************
* RtlEnlargedUnsignedDivide (NTDLL.@)
*
* Divide one 64 bit unsigned integer by a 32 bit unsigned integer, with remainder.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by
* remptr [O] Destination for remainder
*
* RETURNS
* The dividend of a and b. If remptr is non-NULL it is set to the remainder.
*/
UINT WINAPI RtlEnlargedUnsignedDivide( ULONGLONG a, UINT b, UINT *remptr )
{
#if defined(__i386__) && defined(__GNUC__)
UINT ret, rem;
__asm__("divl %4"
: "=a" (ret), "=d" (rem)
: "0" ((UINT)a), "1" ((UINT)(a >> 32)), "g" (b) );
if (remptr) *remptr = rem;
return ret;
#else
UINT ret = a / b;
if (remptr) *remptr = a % b;
return ret;
#endif
}
/******************************************************************************
* RtlExtendedLargeIntegerDivide (NTDLL.@)
*
* Divide one 64 bit integer by a 32 bit integer, with remainder.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by
* rem [O] Destination for remainder
*
* RETURNS
* The dividend of a and b. If rem is non-NULL it is set to the remainder.
*/
LONGLONG WINAPI RtlExtendedLargeIntegerDivide( LONGLONG a, INT b, INT *rem )
{
LONGLONG ret = a / b;
if (rem) *rem = a - b * ret;
return ret;
}
/******************************************************************************
* RtlExtendedIntegerMultiply (NTDLL.@)
*
* Multiply one 64 bit integer by another 32 bit integer.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to multiply a by.
*
* RETURNS
* The product of a and b.
*/
LONGLONG WINAPI RtlExtendedIntegerMultiply( LONGLONG a, INT b )
{
return a * b;
}
/******************************************************************************
* RtlExtendedMagicDivide (NTDLL.@)
*
* Allows replacing a division by a longlong constant with a multiplication by
* the inverse constant.
*
* RETURNS
* (dividend * inverse_divisor) >> (64 + shift)
*
* NOTES
* If the divisor of a division is constant, the constants inverse_divisor and
* shift must be chosen such that inverse_divisor = 2^(64 + shift) / divisor.
* Then we have RtlExtendedMagicDivide(dividend,inverse_divisor,shift) ==
* dividend * inverse_divisor / 2^(64 + shift) == dividend / divisor.
*
* The Parameter inverse_divisor although defined as LONGLONG is used as
* ULONGLONG.
*/
#define LOWER_32(A) ((A) & 0xffffffff)
#define UPPER_32(A) ((A) >> 32)
LONGLONG WINAPI RtlExtendedMagicDivide(
LONGLONG dividend, /* [I] Dividend to be divided by the constant divisor */
LONGLONG inverse_divisor, /* [I] Constant computed manually as 2^(64+shift) / divisor */
INT shift) /* [I] Constant shift chosen to make inverse_divisor as big as possible for 64 bits */
{
ULONGLONG dividend_high;
ULONGLONG dividend_low;
ULONGLONG inverse_divisor_high;
ULONGLONG inverse_divisor_low;
ULONGLONG ah_bl;
ULONGLONG al_bh;
LONGLONG result;
int positive;
if (dividend < 0) {
dividend_high = UPPER_32((ULONGLONG) -dividend);
dividend_low = LOWER_32((ULONGLONG) -dividend);
positive = 0;
} else {
dividend_high = UPPER_32((ULONGLONG) dividend);
dividend_low = LOWER_32((ULONGLONG) dividend);
positive = 1;
} /* if */
inverse_divisor_high = UPPER_32((ULONGLONG) inverse_divisor);
inverse_divisor_low = LOWER_32((ULONGLONG) inverse_divisor);
ah_bl = dividend_high * inverse_divisor_low;
al_bh = dividend_low * inverse_divisor_high;
result = (LONGLONG) ((dividend_high * inverse_divisor_high +
UPPER_32(ah_bl) +
UPPER_32(al_bh) +
UPPER_32(LOWER_32(ah_bl) + LOWER_32(al_bh) +
UPPER_32(dividend_low * inverse_divisor_low))) >> shift);
if (positive) {
return result;
} else {
return -result;
} /* if */
}
/*************************************************************************
* RtlInterlockedCompareExchange64 (NTDLL.@)
*/
#ifdef __GCC_HAVE_SYNC_COMPARE_AND_SWAP_8
LONGLONG WINAPI RtlInterlockedCompareExchange64( LONGLONG *dest, LONGLONG xchg, LONGLONG compare )
{
return __sync_val_compare_and_swap( dest, compare, xchg );
}
#else
__ASM_STDCALL_FUNC(RtlInterlockedCompareExchange64, 20,
"push %ebx\n\t"
__ASM_CFI(".cfi_adjust_cfa_offset 4\n\t")
__ASM_CFI(".cfi_rel_offset %ebx,0\n\t")
"push %esi\n\t"
__ASM_CFI(".cfi_adjust_cfa_offset 4\n\t")
__ASM_CFI(".cfi_rel_offset %esi,0\n\t")
"movl 12(%esp),%esi\n\t"
"movl 16(%esp),%ebx\n\t"
"movl 20(%esp),%ecx\n\t"
"movl 24(%esp),%eax\n\t"
"movl 28(%esp),%edx\n\t"
"lock; cmpxchg8b (%esi)\n\t"
"pop %esi\n\t"
__ASM_CFI(".cfi_same_value %esi\n\t")
__ASM_CFI(".cfi_adjust_cfa_offset -4\n\t")
"pop %ebx\n\t"
__ASM_CFI(".cfi_same_value %ebx\n\t")
__ASM_CFI(".cfi_adjust_cfa_offset -4\n\t")
"ret $20")
#endif
#endif /* _WIN64 */
/******************************************************************************
* RtlLargeIntegerToChar [NTDLL.@]
*
* Convert an unsigned large integer to a character string.
*
* RETURNS
* Success: STATUS_SUCCESS. str contains the converted number
* Failure: STATUS_INVALID_PARAMETER, if base is not 0, 2, 8, 10 or 16.
* STATUS_BUFFER_OVERFLOW, if str would be larger than length.
* STATUS_ACCESS_VIOLATION, if str is NULL.
*
* NOTES
* Instead of base 0 it uses 10 as base.
* Writes at most length characters to the string str.
* Str is '\0' terminated when length allows it.
* When str fits exactly in length characters the '\0' is omitted.
* If value_ptr is NULL it crashes, as the native function does.
*
* DIFFERENCES
* - Accept base 0 as 10 instead of crashing as native function does.
* - The native function does produce garbage or STATUS_BUFFER_OVERFLOW for
* base 2, 8 and 16 when the value is larger than 0xFFFFFFFF.
*/
NTSTATUS WINAPI RtlLargeIntegerToChar(
const ULONGLONG *value_ptr, /* [I] Pointer to the value to be converted */
ULONG base, /* [I] Number base for conversion (allowed 0, 2, 8, 10 or 16) */
ULONG length, /* [I] Length of the str buffer in bytes */
PCHAR str) /* [O] Destination for the converted value */
{
ULONGLONG value = *value_ptr;
CHAR buffer[65];
PCHAR pos;
CHAR digit;
ULONG len;
if (base == 0) {
base = 10;
} else if (base != 2 && base != 8 && base != 10 && base != 16) {
return STATUS_INVALID_PARAMETER;
} /* if */
pos = &buffer[64];
*pos = '\0';
do {
pos--;
digit = value % base;
value = value / base;
if (digit < 10) {
*pos = '0' + digit;
} else {
*pos = 'A' + digit - 10;
} /* if */
} while (value != 0L);
len = &buffer[64] - pos;
if (len > length) {
return STATUS_BUFFER_OVERFLOW;
} else if (str == NULL) {
return STATUS_ACCESS_VIOLATION;
} else if (len == length) {
memcpy(str, pos, len);
} else {
memcpy(str, pos, len + 1);
} /* if */
return STATUS_SUCCESS;
}
/**************************************************************************
* RtlInt64ToUnicodeString (NTDLL.@)
*
* Convert a large unsigned integer to a '\0' terminated unicode string.
*
* RETURNS
* Success: STATUS_SUCCESS. str contains the converted number
* Failure: STATUS_INVALID_PARAMETER, if base is not 0, 2, 8, 10 or 16.
* STATUS_BUFFER_OVERFLOW, if str is too small to hold the string
* (with the '\0' termination). In this case str->Length
* is set to the length, the string would have (which can
* be larger than the MaximumLength).
*
* NOTES
* Instead of base 0 it uses 10 as base.
* If str is NULL it crashes, as the native function does.
*
* DIFFERENCES
* - Accept base 0 as 10 instead of crashing as native function does.
* - Do not return STATUS_BUFFER_OVERFLOW when the string is long enough.
* The native function does this when the string would be longer than 31
* characters even when the string parameter is long enough.
* - The native function does produce garbage or STATUS_BUFFER_OVERFLOW for
* base 2, 8 and 16 when the value is larger than 0xFFFFFFFF.
*/
NTSTATUS WINAPI RtlInt64ToUnicodeString(
ULONGLONG value, /* [I] Value to be converted */
ULONG base, /* [I] Number base for conversion (allowed 0, 2, 8, 10 or 16) */
UNICODE_STRING *str) /* [O] Destination for the converted value */
{
WCHAR buffer[65];
PWCHAR pos;
WCHAR digit;
if (base == 0) {
base = 10;
} else if (base != 2 && base != 8 && base != 10 && base != 16) {
return STATUS_INVALID_PARAMETER;
} /* if */
pos = &buffer[64];
*pos = '\0';
do {
pos--;
digit = value % base;
value = value / base;
if (digit < 10) {
*pos = '0' + digit;
} else {
*pos = 'A' + digit - 10;
} /* if */
} while (value != 0L);
str->Length = (&buffer[64] - pos) * sizeof(WCHAR);
if (str->Length >= str->MaximumLength) {
return STATUS_BUFFER_OVERFLOW;
} else {
memcpy(str->Buffer, pos, str->Length + sizeof(WCHAR));
} /* if */
return STATUS_SUCCESS;
}
#ifdef __i386__
/* those builtin functions use stdcall calling convention, but compilers reference them without stdcall declarations */
#if defined(__MINGW32__) || defined(_MSC_VER)
LONGLONG WINAPI _alldiv( LONGLONG a, LONGLONG b ) asm(__ASM_NAME("_alldiv"));
LONGLONG WINAPI _allmul( LONGLONG a, LONGLONG b ) asm(__ASM_NAME("_allmul"));
LONGLONG WINAPI _allrem( LONGLONG a, LONGLONG b ) asm(__ASM_NAME("_allrem"));
ULONGLONG WINAPI _aulldiv( ULONGLONG a, ULONGLONG b ) asm(__ASM_NAME("_aulldiv"));
ULONGLONG WINAPI _aullrem( ULONGLONG a, ULONGLONG b ) asm(__ASM_NAME("_aullrem"));
#endif
static ULONGLONG udivmod(ULONGLONG a, ULONGLONG b, ULONGLONG *rem)
{
const ULARGE_INTEGER n = { .QuadPart = a };
const ULARGE_INTEGER d = { .QuadPart = b };
DWORD sr, carry, index;
ULARGE_INTEGER q, r;
const unsigned n_uword_bits = 32;
const unsigned n_udword_bits = 64;
/* special cases, X is unknown, K != 0 */
if (n.u.HighPart == 0)
{
if (d.u.HighPart == 0)
{
/* 0 X / 0 X */
if (rem) *rem = n.u.LowPart % d.u.LowPart;
return n.u.LowPart / d.u.LowPart;
}
/* 0 X / K X */
if (rem) *rem = n.u.LowPart;
return 0;
}
/* n.u.HighPart != 0 */
if (d.u.LowPart == 0)
{
if (d.u.HighPart == 0)
{
/* K X / 0 0 */
if (rem) *rem = n.u.HighPart % d.u.LowPart;
return n.u.HighPart / d.u.LowPart;
}
/* d.u.HighPart != 0 */
if (n.u.LowPart == 0) {
/* K 0 / K 0 */
if (rem)
{
r.u.HighPart = n.u.HighPart % d.u.HighPart;
r.u.LowPart = 0;
*rem = r.QuadPart;
}
return n.u.HighPart / d.u.HighPart;
}
/* K K / K 0 */
if ((d.u.HighPart & (d.u.HighPart - 1)) == 0) /* if d is a power of 2 */
{
if (rem)
{
r.u.LowPart = n.u.LowPart;
r.u.HighPart = n.u.HighPart & (d.u.HighPart - 1);
*rem = r.QuadPart;
}
BitScanForward(&index, d.u.HighPart);
return n.u.HighPart >> index;
}
/* K K / K 0 */
BitScanReverse(&index, d.u.HighPart);
BitScanReverse(&sr, n.u.HighPart);
sr -= index;
/* 0 <= sr <= n_uword_bits - 2 or sr large */
if (sr > n_uword_bits - 2)
{
if (rem) *rem = n.QuadPart;
return 0;
}
++sr;
/* 1 <= sr <= n_uword_bits - 1 */
/* q.QuadPart = n.QuadPart << (n_udword_bits - sr); */
q.u.LowPart = 0;
q.u.HighPart = n.u.LowPart << (n_uword_bits - sr);
/* r.QuadPart = n.QuadPart >> sr; */
r.u.HighPart = n.u.HighPart >> sr;
r.u.LowPart = (n.u.HighPart << (n_uword_bits - sr)) | (n.u.LowPart >> sr);
}
else /* d.u.LowPart != 0 */
{
if (d.u.HighPart == 0)
{
/* K X / 0 K */
if ((d.u.LowPart & (d.u.LowPart - 1)) == 0) /* if d is a power of 2 */
{
if (rem) *rem = n.u.LowPart & (d.u.LowPart - 1);
if (d.u.LowPart == 1) return n.QuadPart;
BitScanForward(&sr, d.u.LowPart);
q.u.HighPart = n.u.HighPart >> sr;
q.u.LowPart = (n.u.HighPart << (n_uword_bits - sr)) | (n.u.LowPart >> sr);
return q.QuadPart;
}
BitScanReverse(&index, d.u.LowPart);
BitScanReverse(&sr, n.u.HighPart);
sr = 1 + n_uword_bits + sr - index;
/* 2 <= sr <= n_udword_bits - 1
* q.QuadPart = n.QuadPart << (n_udword_bits - sr);
* r.QuadPart = n.QuadPart >> sr; */
if (sr == n_uword_bits)
{
q.u.LowPart = 0;
q.u.HighPart = n.u.LowPart;
r.u.HighPart = 0;
r.u.LowPart = n.u.HighPart;
}
else if (sr < n_uword_bits) /* 2 <= sr <= n_uword_bits - 1 */
{
q.u.LowPart = 0;
q.u.HighPart = n.u.LowPart << (n_uword_bits - sr);
r.u.HighPart = n.u.HighPart >> sr;
r.u.LowPart = (n.u.HighPart << (n_uword_bits - sr)) | (n.u.LowPart >> sr);
}
else /* n_uword_bits + 1 <= sr <= n_udword_bits - 1 */
{
q.u.LowPart = n.u.LowPart << (n_udword_bits - sr);
q.u.HighPart = (n.u.HighPart << (n_udword_bits - sr)) |
(n.u.LowPart >> (sr - n_uword_bits));
r.u.HighPart = 0;
r.u.LowPart = n.u.HighPart >> (sr - n_uword_bits);
}
}
else
{
/* K X / K K */
BitScanReverse(&index, d.u.HighPart);
BitScanReverse(&sr, n.u.HighPart);
sr -= index;
/* 0 <= sr <= n_uword_bits - 1 or sr large */
if (sr > n_uword_bits - 1)
{
if (rem) *rem = n.QuadPart;
return 0;
}
++sr;
/* 1 <= sr <= n_uword_bits
* q.QuadPart = n.QuadPart << (n_udword_bits - sr); */
q.u.LowPart = 0;
if (sr == n_uword_bits)
{
q.u.HighPart = n.u.LowPart;
r.u.HighPart = 0;
r.u.LowPart = n.u.HighPart;
}
else
{
q.u.HighPart = n.u.LowPart << (n_uword_bits - sr);
r.u.HighPart = n.u.HighPart >> sr;
r.u.LowPart = (n.u.HighPart << (n_uword_bits - sr)) | (n.u.LowPart >> sr);
}
}
}
/* Not a special case
* q and r are initialized with:
* q.QuadPart = n.QuadPart << (n_udword_bits - sr);
* r.QuadPart = n.QuadPart >> sr;
* 1 <= sr <= n_udword_bits - 1 */
carry = 0;
for (; sr > 0; --sr)
{
LONGLONG s;
/* r:q = ((r:q) << 1) | carry */
r.u.HighPart = (r.u.HighPart << 1) | (r.u.LowPart >> (n_uword_bits - 1));
r.u.LowPart = (r.u.LowPart << 1) | (q.u.HighPart >> (n_uword_bits - 1));
q.u.HighPart = (q.u.HighPart << 1) | (q.u.LowPart >> (n_uword_bits - 1));
q.u.LowPart = (q.u.LowPart << 1) | carry;
/* if (r.QuadPart >= d.QuadPart)
* {
* r.QuadPart -= d.QuadPart;
* carry = 1;
* }
*/
s = (LONGLONG)(d.QuadPart - r.QuadPart - 1) >> (n_udword_bits - 1);
carry = s & 1;
r.QuadPart -= d.QuadPart & s;
}
q.QuadPart = (q.QuadPart << 1) | carry;
if (rem) *rem = r.QuadPart;
return q.QuadPart;
}
/******************************************************************************
* _alldiv (NTDLL.@)
*
* Divide two 64 bit unsigned integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by.
*
* RETURNS
* The dividend of a and b.
*/
LONGLONG WINAPI _alldiv( LONGLONG a, LONGLONG b )
{
LONGLONG s_a = a >> 63; /* s_a = a < 0 ? -1 : 0 */
LONGLONG s_b = b >> 63; /* s_b = b < 0 ? -1 : 0 */
a = (a ^ s_a) - s_a; /* negate if s_a == -1 */
b = (b ^ s_b) - s_b; /* negate if s_b == -1 */
s_a ^= s_b; /* sign of quotient */
return (udivmod(a, b, NULL) ^ s_a) - s_a; /* negate if s_a == -1 */
}
/******************************************************************************
* _allmul (NTDLL.@)
*
* Multiply two 64 bit integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to multiply a by.
*
* RETURNS
* The product of a and b.
*/
LONGLONG WINAPI _allmul( LONGLONG a, LONGLONG b )
{
LARGE_INTEGER x = { .QuadPart = a };
LARGE_INTEGER y = { .QuadPart = b };
LARGE_INTEGER r;
unsigned int t;
const int bits_in_word_2 = 16;
const unsigned int lower_mask = ~0u >> bits_in_word_2;
r.u.LowPart = (x.u.LowPart & lower_mask) * (y.u.LowPart & lower_mask);
t = r.u.LowPart >> bits_in_word_2;
r.u.LowPart &= lower_mask;
t += (x.u.LowPart >> bits_in_word_2) * (y.u.LowPart & lower_mask);
r.u.LowPart += (t & lower_mask) << bits_in_word_2;
r.u.HighPart = t >> bits_in_word_2;
t = r.u.LowPart >> bits_in_word_2;
r.u.LowPart &= lower_mask;
t += (y.u.LowPart >> bits_in_word_2) * (x.u.LowPart & lower_mask);
r.u.LowPart += (t & lower_mask) << bits_in_word_2;
r.u.HighPart += t >> bits_in_word_2;
r.u.HighPart += (x.u.LowPart >> bits_in_word_2) * (y.u.LowPart >> bits_in_word_2);
r.u.HighPart += x.u.HighPart * y.u.LowPart + x.u.LowPart * y.u.HighPart;
return r.QuadPart;
}
/******************************************************************************
* _allrem (NTDLL.@)
*
* Calculate the remainder after dividing two 64 bit integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by.
*
* RETURNS
* The remainder of a divided by b.
*/
LONGLONG WINAPI _allrem( LONGLONG a, LONGLONG b )
{
LONGLONG s = b >> 63; /* s = b < 0 ? -1 : 0 */
ULONGLONG r;
b = (b ^ s) - s; /* negate if s == -1 */
s = a >> 63; /* s = a < 0 ? -1 : 0 */
a = (a ^ s) - s; /* negate if s == -1 */
udivmod(a, b, &r);
return ((LONGLONG)r ^ s) - s; /* negate if s == -1 */
}
/******************************************************************************
* _aulldiv (NTDLL.@)
*
* Divide two 64 bit unsigned integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by.
*
* RETURNS
* The dividend of a and b.
*/
ULONGLONG WINAPI _aulldiv( ULONGLONG a, ULONGLONG b )
{
return udivmod(a, b, NULL);
}
/******************************************************************************
* _allshl (NTDLL.@)
*
* Shift a 64 bit integer to the left.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to shift a by to the left.
*
* RETURNS
* The left-shifted value.
*/
LONGLONG WINAPI _allshl( LONGLONG a, LONG b )
{
return a << b;
}
/******************************************************************************
* _allshr (NTDLL.@)
*
* Shift a 64 bit integer to the right.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to shift a by to the right.
*
* RETURNS
* The right-shifted value.
*/
LONGLONG WINAPI _allshr( LONGLONG a, LONG b )
{
return a >> b;
}
/******************************************************************************
* _alldvrm (NTDLL.@)
*
* Divide two 64 bit integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by.
*
* RETURNS
* Returns the quotient of a and b in edx:eax.
* Returns the remainder of a and b in ebx:ecx.
*/
__ASM_GLOBAL_FUNC( _alldvrm,
"pushl %ebp\n\t"
__ASM_CFI(".cfi_adjust_cfa_offset 4\n\t")
__ASM_CFI(".cfi_rel_offset %ebp,0\n\t")
"movl %esp,%ebp\n\t"
__ASM_CFI(".cfi_def_cfa_register %ebp\n\t")
"pushl 20(%ebp)\n\t"
"pushl 16(%ebp)\n\t"
"pushl 12(%ebp)\n\t"
"pushl 8(%ebp)\n\t"
"call " __ASM_NAME("_allrem") "\n\t"
"movl %edx,%ebx\n\t"
"pushl %eax\n\t"
"pushl 20(%ebp)\n\t"
"pushl 16(%ebp)\n\t"
"pushl 12(%ebp)\n\t"
"pushl 8(%ebp)\n\t"
"call " __ASM_NAME("_alldiv") "\n\t"
"popl %ecx\n\t"
"leave\n\t"
__ASM_CFI(".cfi_def_cfa %esp,4\n\t")
__ASM_CFI(".cfi_same_value %ebp\n\t")
"ret $16" )
/******************************************************************************
* _aullrem (NTDLL.@)
*
* Calculate the remainder after dividing two 64 bit unsigned integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by.
*
* RETURNS
* The remainder of a divided by b.
*/
ULONGLONG WINAPI _aullrem( ULONGLONG a, ULONGLONG b )
{
ULONGLONG r;
udivmod(a, b, &r);
return r;
}
/******************************************************************************
* _aullshr (NTDLL.@)
*
* Shift a 64 bit unsigned integer to the right.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to shift a by to the right.
*
* RETURNS
* The right-shifted value.
*/
ULONGLONG WINAPI _aullshr( ULONGLONG a, LONG b )
{
return a >> b;
}
/******************************************************************************
* _aulldvrm (NTDLL.@)
*
* Divide two 64 bit unsigned integers.
*
* PARAMS
* a [I] Initial number.
* b [I] Number to divide a by.
*
* RETURNS
* Returns the quotient of a and b in edx:eax.
* Returns the remainder of a and b in ebx:ecx.
*/
__ASM_GLOBAL_FUNC( _aulldvrm,
"pushl %ebp\n\t"
__ASM_CFI(".cfi_adjust_cfa_offset 4\n\t")
__ASM_CFI(".cfi_rel_offset %ebp,0\n\t")
"movl %esp,%ebp\n\t"
__ASM_CFI(".cfi_def_cfa_register %ebp\n\t")
"pushl 20(%ebp)\n\t"
"pushl 16(%ebp)\n\t"
"pushl 12(%ebp)\n\t"
"pushl 8(%ebp)\n\t"
"call " __ASM_NAME("_aullrem") "\n\t"
"movl %edx,%ebx\n\t"
"pushl %eax\n\t"
"pushl 20(%ebp)\n\t"
"pushl 16(%ebp)\n\t"
"pushl 12(%ebp)\n\t"
"pushl 8(%ebp)\n\t"
"call " __ASM_NAME("_aulldiv") "\n\t"
"popl %ecx\n\t"
"leave\n\t"
__ASM_CFI(".cfi_def_cfa %esp,4\n\t")
__ASM_CFI(".cfi_same_value %ebp\n\t")
"ret $16" )
#endif /* __i386__ */