我编写了一些代码来测试 fsqrt 函数,但结果对我来说并不完全有意义。这是代码(delphi 中):
uses
mmsystem;
var
rand:longint=123456789;
function rng:longint;
asm
imul eax,[rand],$08088405
inc eax
mov [rand],eax
end;
function int_sqrt(adata:longint):longint;
asm
fnstcw word([esp-2])
// mov word([esp-4]),$1f3f // 80bit precision
mov word([esp-4]),$1c3f // 24bit precision
fldcw word([esp-4])
mov [esp-8],eax
fild longint([esp-8])
fsqrt
fistp longint([esp-8])
mov eax,[esp-8]
fldcw word([esp-2])
end;
procedure TForm1.FormCreate(Sender: TObject);
var
start,i,r,s1,s2:longint;
time0,time1:longint;
begin
timebeginperiod(1);
time0:=timegettime;
start:=1000000000;
for i:=(start+0) to (start+100000000) do begin
//r:=i;
r:=abs(rng);
// r:=2134567890;
// r:=$7fffffff;
s1:=int_sqrt(r);
s2:=trunc(sqrt(r));
if s1<>s2 then
showmessage('error: '+inttostr(r)+'/'+inttostr(s1)+'/'+inttostr(s2));
end;
time1:=timegettime;
timeendperiod(1);
showmessage('Milliseconds: '+inttostr(time1-time0));
end;
很简单,我正在寻找 int 的平方根。在 int_sqrt 中,其中一个精度行使 x87 使用 24 位精度作为 sqrt 精度,另一个使用 64 位精度。正如预期的那样,24 位版本的速度大大加快(10-20%,具体取决于输入)。
但这就是问题所在。我还没有找到一个 32 位(实际上是 31 位,最后一位是未使用的符号)int 在使用 24 位精度时返回错误结果!
到目前为止,我唯一的理论是,只有最终结果取决于精度,而不取决于源或任何中间缓冲区。这是有道理的,因为 31 位 int 的平方根的最大结果大小是 16 位。
是这样的吗?