The following code:
while Assigned(p) do begin
np:= p^.next;
h:= leaf_hash(p^.data); <<-- inline routine
h:= h mod nhashprime;
p^.next:= nhashtab[h];
nhashtab[h]:= p;
p:= np;
end; { while }
Produces the following assembly:
hlife.pas.605: h:= leaf_hash(p^.data);
00000000005D4602 498B4018 mov rax,[r8+$18]
00000000005D4606 48C1E830 shr rax,$30
00000000005D460A 498B5018 mov rdx,[r8+$18]
00000000005D460E 48C1EA20 shr rdx,$20
00000000005D4612 81E2FFFF0000 and edx,$0000ffff
00000000005D4618 4D8B5818 mov r11,[r8+$18]
00000000005D461C 49C1EB10 shr r11,$10
00000000005D4620 4181E3FFFF0000 and r11d,$0000ffff
00000000005D4627 418B7018 mov esi,[r8+$18]
00000000005D462B 81E6FFFF0000 and esi,$0000ffff
00000000005D4631 488D34F6 lea rsi,[rsi+rsi*8]
00000000005D4635 4403DE add r11d,esi
00000000005D4638 4F8D1CDB lea r11,[r11+r11*8]
00000000005D463C 4103D3 add edx,r11d
00000000005D463F 488D14D2 lea rdx,[rdx+rdx*8]
00000000005D4643 03C2 add eax,edx
hlife.pas.606: h:= h mod nhashprime;
00000000005D4645 8BC0 mov eax,eax <<--- Why is there a NOP here?
00000000005D4647 4C63DB movsxd r11,rbx
00000000005D464A 4899 cwd
00000000005D464C 49F7FB idiv r11
00000000005D464F 488BC2 mov rax,rdx
hlife.pas.607: p^.next:= nhashtab[h];
00000000005D4652 488B5538 mov rdx,[rbp+$38]
Delphi inserts a NOP before the nhashtab[h]:= p;
line.
If the leaf_hash
function had been a normal function it would have made sense.
(no not really because the RET would still return to [5D4645] executing the nop)
But now this is not a jump target.
So I'm (just) curious, Why does it do this?
[EDIT] : SSCCE
OK I've got a SSCCE up {it's not very short but it will have to do.
Note the compiler settings (Debug + Win64)
unit Unit16;
interface
uses
Winapi.Windows, Winapi.Messages, System.SysUtils, System.Variants, System.Classes, Vcl.Graphics,
Vcl.Controls, Vcl.Forms, Vcl.Dialogs, Vcl.StdCtrls;
type
pnode = ^node;
tflavour = (tnode, tleaf, tleaf64);
node = record
case flavour: tflavour of
tnode: (next: pnode; (* hash link *)
nw, ne, sw, se: pnode; (* constant; nw not= 0 means nonleaf *)
res: pnode); (* cache *)
tleaf: (next1: pnode; (* hash link *)
isnode: pnode; (* must always be zero for leaves *)
nw1, ne1, sw1, se1: word; (* constant *)
res1, res2: word; (* constant *)
);
tleaf64: (next2: pnode; (* hash link *)
isnode1: pnode; (* must always be zero for leaves *)
data: Uint64; (* constant *)
res1a, res2a: word; (* constant *)
)
end;
ppnode = array of pnode;
THashBox = class(TPersistent)
strict private
leafhashpop: integer;
leafhashlimit: integer;
leafhashprime: integer;
leafhashtab: ppnode;
nodehashpop: integer;
nodehashlimit: integer;
nodehashprime: integer;
nodehashtab: ppnode;
private
TotalTime, Occurrences: Uint64;
StartTime, EndTime: Uint64;
procedure resize_leaves();
public
constructor Create;
end;
TForm16 = class(TForm)
Button1: TButton;
procedure Button1Click(Sender: TObject);
private
HashBox: THashBox;
public
end;
var
Form16: TForm16;
implementation
{$R *.dfm}
const
maxmem = 2000*1000*1000; {2GB}
var
alloced: cardinal;
function rdtsc: int64; assembler;
asm
{ xor eax,eax;
push rbx
cpuid
pop rbx }
rdtsc
end;
function node_hash(n: pnode): cardinal; { inline; } assembler; overload;
// var
// a: pnativeint;
// begin
// Result:= nativeint(n^.se) + 3 * (nativeint(n^.sw) + 3 * (nativeint(n^.ne) + 3 * nativeint(n^.nw) + 3));
asm
mov eax,[rcx+node.nw]
lea eax,[eax+eax*2+3]
add eax,[rcx+node.ne]
lea eax,[eax+eax*2]
add eax,[rcx+node.sw]
lea eax,[eax+eax*2]
add eax,[rcx+node.se]
end;
function leaf_hash(a, b, c, d: cardinal): cardinal; inline; overload;
begin
Result:= (d + 9 * (c + 9 * (b + 9 * a)))
end;
function leaf_hash(data: Uint64): cardinal; inline; overload;
begin
// Result:= d + 9 * (c + 9 * (b + 9 * a));
Result:= ((data shr 48) + 9 * (((data shr 32) and $FFFF) + 9 * (((data shr 16) and $FFFF) + 9 * (data and $FFFF))));
Inc(Result);
end;
procedure TForm16.Button1Click(Sender: TObject);
begin
HashBox:= THashBox.Create;
Hashbox.resize_leaves;
end;
function nextprime(old: integer): integer;
begin
Result:= 1009;
end;
constructor THashBox.Create;
begin
leafhashprime:= 7;
SetLength(leafhashtab, leafhashprime);
end;
procedure THashBox.resize_leaves();
var
i, i1, i2: integer;
nhashprime: Cardinal;
p: pnode;
nhashtab: ppnode;
np: pnode;
h: Integer;
n, n2: integer;
diff1, diff2: integer;
begin
nhashprime:= nextprime(4 * leafhashprime);
if (nhashprime * sizeof(pnode) > maxmem - alloced) then begin
leafhashlimit:= 2000 * 1000 * 1000;
exit;
end;
(*
* Don't let the hash table buckets take more than 4% of the
* memory. If we're starting to strain memory, let the buckets
* fill up a bit more.
*)
if (nhashprime > maxmem div 100) then begin
nhashprime:= nextprime(maxmem div 100);
if (nhashprime = leafhashprime) then begin
leafhashlimit:= 2000 * 1000 * 1000;
exit;
end;
end;
SetLength(nhashtab, nhashprime); //make a new table, do not resize the existing one.
alloced:= alloced + sizeof(pnode) * (nhashprime - leafhashprime);
diff1:= maxint;
for i1:= 0 to 100 do begin
n:= 0;
StartTime:= rdtsc;
for i:= 0 to leafhashprime - 1 do begin
p:= leafhashtab[i];
if Assigned(p) then begin
h:= node_hash(p);
h:= h mod nhashprime;
inc(n, h);
end;
end;
EndTime:= rdtsc;
if ((EndTime - StartTime) < diff1) then diff1:= (EndTime - StartTime);
end;
diff2:= maxint;
for i1:= 0 to 100 do begin
n2:= 0;
StartTime:= rdtsc;
for i:= 0 to leafhashprime - 1 do begin
p:= leafhashtab[i];
if Assigned(p) then begin
inc(n2);
end;
end;
EndTime:= rdtsc;
if (endtime - starttime) < diff2 then diff2:= endtime - starttime;
end;
TotalTime:= diff1 - diff2;
if n <> n2 then Occurrences:= nhashprime;
for i:= 0 to leafhashprime - 1 do begin
// for (p=hashtab[i]; p;) {
p:= leafhashtab[i];
while Assigned(p) do begin <<--- put a breakpoint here
np:= p^.next;
h:= leaf_hash(p^.data);
h:= h mod nhashprime;
p^.next:= nhashtab[h];
nhashtab[h]:= p;
p:= np;
end; { while }
end; { for i }
// free(hashtab);
leafhashtab:= nhashtab;
leafhashprime:= nhashprime;
leafhashlimit:= leafhashprime;
end;
end.
You will see this disassembly:
Unit16.pas.196: h:= h mod nhashprime;
000000000059CE4B 4863C0 movsxd rax,rax
000000000059CE4E 448B5528 mov r10d,[rbp+$28]
000000000059CE52 458BD2 mov r10d,r10d <<--- weird NOP here
000000000059CE55 4899 cwd
000000000059CE57 49F7FA idiv r10
000000000059CE5A 488BC2 mov rax,rdx
Unit16.pas.197: p^.next:= nhashtab[h];
000000000059CE5D 488B5538 mov rdx,[rbp+$38]
The answer is that
MOV EAX,EAX
Is not a no-op
On x64 manipulating the lower 32 bits of a 64 bit register will zero the upper 32 bits.
So the above instruction should really be read as:
MOVZX RAX,EAX
According to AMD
Results of 32-bit operations are implicitly zero extended to 64-bit values.
IMHO this is not a nop
for alignement, but it sounds to me just like un-optimized generated code, and wrong signing of your own variables.
h:= h mod nhashprime;
May be split into:
mov eax,eax new h = eax, old h = eax // which does not mean anything
movsxd r11,rbx convert with sign nhashprime stored in rbx into temp registry r11
cwd signed rax into rdx:rax
idiv r11 signed divide rdx:rax by r11 -> rax=quotient, rdx=remainder
mov rax,rdx store remainder rdx into rax
Did you try enabling the code generation optimization? I suppose it will fix the content of mov eax,eax
.
But your original code is also sub-optimized. You should use unsigned arithmetic in your case.
And, you should better use a power of two nhashprime
, the compute a simple and
binary operation instead of a slow division:
var h, nhashprimeMask: cardinal; // use UNSIGNED arithmetic here!
// here nhashprime is a POWER OF TWO (128,256,512,1024,2048...)
nhashprimeMask := nhashprime-1; // compute the mask
while Assigned(p) do begin
np:= p^.next;
h:= leaf_hash(p^.data) and nhashprimeMask;
p^.next:= nhashtab[h];
nhashtab[h]:= p;
p:= np;
end; { while }
This code will be much faster, and should much better compile.
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With