1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
TITLE DivU64x32.asm: 64-bit division function for IA-32
;------------------------------------------------------------------------------
;
; Copyright (c) 2004, Intel Corporation
; All rights reserved. This program and the accompanying materials
; are licensed and made available under the terms and conditions of the BSD License
; which accompanies this distribution. The full text of the license may be found at
; http://opensource.org/licenses/bsd-license.php
;
; THE PROGRAM IS DISTRIBUTED UNDER THE BSD LICENSE ON AN "AS IS" BASIS,
; WITHOUT WARRANTIES OR REPRESENTATIONS OF ANY KIND, EITHER EXPRESS OR IMPLIED.
;
; Module Name:
;
; DivU64x32.asm
;
; Abstract:
;
; 64-bit division function for IA-32
;
;------------------------------------------------------------------------------
.686P
.XMM
.MODEL SMALL
.CODE
DivU64x32 PROTO C Dividend: QWORD, Divisor: DWORD, Remainder: DWORD
DivU64x32 PROC C Dividend: QWORD, Divisor: DWORD, Remainder: DWORD
;------------------------------------------------------------------------------
; UINT64
; DivU64x32 (
; IN UINT64 Dividend,
; IN UINTN Divisor,
; OUT UINTN *Remainder OPTIONAL
; )
;
; Routine Description:
;
; This routine allows a 64 bit value to be divided with a 32 bit value returns
; 64bit result and the Remainder.
;
; Arguments:
;
; Dividend - dividend
; Divisor - divisor
; Remainder - buffer for remainder
;
; Returns:
;
; Dividend / Divisor
; Remainder = Dividend mod Divisor
;
; N.B. only works for 31bit divisors!!
;------------------------------------------------------------------------------
push ecx
;
; let edx contain the intermediate result of remainder
;
xor edx, edx
mov ecx, 64
_DivU64x32_Wend:
shl dword ptr Dividend[0], 1
rcl dword ptr Dividend[4], 1
rcl edx, 1
;
; If intermediate result of remainder is larger than
; or equal to divisor, then set the lowest bit of dividend,
; and subtract divisor from intermediate remainder
;
cmp edx, Divisor
jb _DivU64x32_Cont
bts dword ptr Dividend[0], 0
sub edx, Divisor
_DivU64x32_Cont:
loop _DivU64x32_Wend
cmp Remainder, 0
je _DivU64x32_Done
mov eax, Remainder
mov dword ptr [eax], edx
_DivU64x32_Done:
mov eax, dword ptr Dividend[0]
mov edx, dword ptr Dividend[4]
pop ecx
ret
DivU64x32 ENDP
END
|