Версия для печати темы

Автор: TarasBer 4.01.2009 2:24

Итак у нас, например, такая задача. Вещественное число из диапазона от 0 до 65535 преобразовать в word взятием целой части.

Это программа для дельфы.

program TruncTest;

{$APPTYPE CONSOLE}

uses
  SysUtils, Windows;

var
  i: longint;
  d, d1: double;
  d_w: array [0 .. 3] of word absolute d;
  k: word;
  T: longint;

begin
  d := 1234.5678;

  T := GetTickCount;
  for i := 0 to $FFFFFF do
    k := trunc(d);
  T := GetTickCount - T;
  WriteLn('Standart trunc: ', T);

  T := GetTickCount;
  for i := 0 to $FFFFFF do
    k :=  (word(d_w[3] shl 11) or (d_w[2] shr 5) or $8000)
      shr (14 - (d_w[3] shr 4 and $03FF));
  T := GetTickCount - T;
  WriteLn('My trunc: ', T);

  T := GetTickCount;
  for i := 0 to $FFFFFF do
    d1 := 1/d;
  T := GetTickCount - T;
  WriteLn('Floating point division: ', T);

  ReadLn;
end.

var
  i: longint;
  d, d1: double;
  k: integer;
  Time: longint absolute $0040: $006C;
  T: longint;
  d_w: array [0 .. 3] of word absolute d;

begin
  d := 1234.5678;

  T := Time; while T = Time do; T := Time;
  for i := 0 to $FFFFFF do
    k := trunc(d);
  T := Time - T;
  WriteLn('Standart trunc: ', T);

  T := Time; while T = Time do; T := Time;
  for i := 0 to $FFFFFF do
    k := ((d_w[3] shl 11) or (d_w[2] shr 5) or $8000)
      shr (14 - (d_w[3] shr 4 and $03FF));
  T := Time - T;
  WriteLn('My trunc: ', T);

  T := Time; while T = Time do; T := Time;
  for i := 0 to $FFFFFF do
    d1 := 1/d;
  T := Time - T;
  WriteLn('Floating point division: ', T);
end.

Автор: Account 4.01.2009 3:21

Лечиться собственным кодом

Автор: volvo 4.01.2009 3:39

  randomize;
  T := GetTickCount;
  for i := 0 to $FFFFFF do
    k := trunc(1 + random * 1000);
  T := GetTickCount - T;
  WriteLn('Standart trunc: ', T);

  T := GetTickCount;
  for i := 0 to $FFFFFF do
    k := FastTrunc(1 + random * 1000);
  T := GetTickCount - T;
  WriteLn('FastTrunc: ', T);

  T := GetTickCount;
  for i := 0 to $FFFFFF do begin
    d := 1 + random * 1000;
    k :=  (word(d_w[3] shl 11) or (d_w[2] shr 5) or $8000)
      shr (14 - (d_w[3] shr 4 and $03FF));
  end;
  T := GetTickCount - T;
  WriteLn('My trunc: ', T);

  T := GetTickCount;
  for i := 0 to $FFFFFF do
    d1 := 1 / (1 + random * 1000);
  T := GetTickCount - T;
  WriteLn('Floating point division: ', T);

function myTrunc(X: single): integer;
asm
  CVTTSS2SI  eax, [ebp+offset X]
end;

Автор: TarasBer 4.01.2009 4:07

Дело в том, что я продолжаю гоняться за наивной детской мечтой реализовать полноценное 3д на Турбо Паскале. И в нём разница куда более заметна, чем в Дельфе.
И поэтому я не знаю, что такое FastTrunc, откуда он взялся, и зачем нужен тогда Trunc, если есть FastTrunc, а также, что за команда CVTTSS2SI.
Ещё вопрос - а стандартный trunc - это случайно не функция (всмысле что не макрос)? Может в этом всё дело?

И всё равно - непонятно, почему такая, казалось бы, простая команда не работает на порядок быстрее такой сложной команды, как деление?

Автор: volvo 4.01.2009 4:14

Автор: TarasBer 4.01.2009 4:23

Цитата(volvo @ 4.01.2009 0:14)

Автор: volvo 4.01.2009 5:39

var
  arr: array[0 .. $FFFFFF] of double;
...
  randomize;
  for I := 0 to $FFFFFF do begin
    arr[i] := random * 1000;
  end;

Автор: volvo 4.01.2009 16:09

Кстати, вот вариант от Merlyn-а, который перекрывает по скорости и "Standard Trunc" и "My Trunc" из первоначального кода:

const Half: double = 0.5 ;
function chomp(Flt : double) : integer;
var result: integer;
begin
  asm
    fld Flt
    fsub Half
    fistp result
  end;
  chomp := result;
end;

Автор: TarasBer 4.01.2009 19:50

const Half: double = 0.5 ;

asm
    fld d
    fsub Half
    fistp k
end