Вычисление корня "вручную"

[northener]

Сообщение от digitalis

Не думаю, что компилятор делает различие между двумя описаниями одного и того же действия.

В стародавние времена точно была разница. Но имхо уже достаточно давно машинный код будет одинаковый.

[Arigato]

Сообщение от digitalis

По логике, там даже не будет вызова процедуры, а in-line вставка.

А inc это и не процедура по сути. Есть набор инструкций, которые выглядят как процедуры/функции, вызываются как процедуры/функции, но функциями не являются. Скорее языковые конструкции.

[northener]

Сообщение от Arigato

А inc это и не процедура по сути. Есть набор инструкций, которые выглядят как процедуры/функции, вызываются как процедуры/функции, но функциями не являются. Скорее языковые конструкции.

Это часть той самой Compile Magic.

[digitalis]

Сообщение от Arigato

Есть набор инструкций, которые выглядят как процедуры/функции, вызываются как процедуры/функции, но функциями не являются.

Вызглядят и вызываются как... где, в исходнике - или сгенерированном коде? Если второе, то будучи функциями как таковыми, они будут вызываться стандартным образом: сохранение адреса возврата, PSW и используемых регистров, передача управления коду функции и восстановление сохранённого. В противном случае - если они не функции, а только похожи - при каждом "псевдовызове" будет генрироваться код, реализующий искомое. Тот самый in-line.

[Arigato]

Сообщение от digitalis

где, в исходнике - или сгенерированном коде?

В исходнике, конечно. Просто синтаксически похоже на процедуру/функцию, но ею не является. К примеру, в Лазарусе есть нажать Ctrl и навести на любую процедуру/функцию, она подсвечивается как ссылка. А вот такие как inc вообще никак не реагируют, как какой-нибудь for или if...

[Kronos913]

Сообщение от NetSpace

теперь запили отдельную библиотеку *.pas и скинь всем в тему. пусть используют, кто занимается расчётами.
и кстати, эту строки b:=b+k; тоже можно сделать через inc() - быстрее будет, там надо через запятую параметр k написать.

Две процедуры с interface

Function Sqrt64(var a:Int64):LongWord;
Function SqrtWord(const w:LongWord):LongWord;

Procedure SqrtEngine;
asm
  mov esi, 1073741824 {2^30}
  mov ebx, 32768      {2^15}

  @Begin1:
    cmp esi, edi
      jbe @End1

    shr esi, 2
    shr ebx, 1

    jmp @Begin1
  @End1:

  mov esi, 0

  @Begin2:
    mov ecx, esi
    add ecx, ebx
    mov eax, ecx
    mov edx, 0

    mul ecx

    cmp eax, edi
      ja @Big1

    add esi, ebx

    @Big1:
    shr ebx, 1

    cmp ebx, 0
      ja @Begin2

  mov eax, esi
end;

Function SqrtWord(const w:LongWord):LongWord;
asm
  push ebx
  push edi
  push esi

  mov edi, eax
  call SqrtEngine

  pop esi
  pop edi
  pop ebx
end;

Function Sqrt64(var a:Int64):LongWord;
asm
  push ebx
  push edi
  push esi

  mov edi, [eax]
  mov esi, [eax+4] {Çàãðóçèëè íà÷àëüíîå ÷èñëî}

  cmp esi, 0
    ja @Start1

  call SqrtEngine
  jmp @ExitFast

  @Start1:
  push ebp
  mov ebp, 1073741824
  mov ebx, 2147483648

  @Begin0:
    cmp ebp, esi
      jbe @End0

    shr ebp, 2
    shr ebx, 1

    jmp @Begin0
  @End0:

  mov ebp, 0
  @Begin1:
    mov ecx, ebp
    add ecx, ebx
    mov eax, ecx
    mov edx, 0

    mul ecx

    cmp edx, esi
      ja @Big1
      jb @Go1

    cmp eax, edi
      ja @Big1

    @Go1:
    add ebp, ebx

    @Big1:
    shr ebx, 1

    cmp ebx, 0
      ja @Begin1

  mov eax, ebp

  pop ebp

  @ExitFast:
  pop esi
  pop edi
  pop ebx
end;

[Vapaamies]

Сообщение от Kronos913

Есть переменная Int64, но я не нашел ее "положительного" аналога в виде некоего Word64 (Если говорить о 32 разрядном Delphi 7)

В Delphi 7 появился «недокументированный» тип UInt64. В чем именно его недоделанность, не изучал.

Сообщение от Arigato

А inc это и не процедура по сути. Есть набор инструкций, которые выглядят как процедуры/функции, вызываются как процедуры/функции, но функциями не являются. Скорее языковые конструкции.

В богомерзком Си оно называется заморским словом интринсик. В Delphi условно к ним также относятся Break и Continue, управляющие выполнением цикла. Это всё тлетворное влияние Turbo Pascal 7.0, с которого всё пошло.

Сообщение от Kronos913

Итог (на картинке): с числами до 2^32 мой код работает быстрее, а вот с числами после 2^32 - медленнее

Код с GetTickCount не вызывает доверия. Замеры скорости выполнения должны делаться через QueryPerformanceCounter. Вытащил код своего класса-обертки:

type
  TPerformanceCounter = class
  private
    FFrequency: UInt64;
    function GetValue: UInt64;
  public
    constructor Create;
    function MillisecondsFrom(StartValue: UInt64): Double;
    function MillisecondsBetween(Value1, Value2: UInt64): Double;

    property Frequency: UInt64 read FFrequency;
    property Value: UInt64 read GetValue;
  end;

constructor TPerformanceCounter.Create;
var
  Freq: Int64;
begin
  if QueryPerformanceFrequency(Freq) then // types of format and actual parameters must be equal
    FFrequency := Freq;
end;

function TPerformanceCounter.MillisecondsFrom(StartValue: UInt64): Double;
begin
  Result := (GetValue - StartValue) / FFrequency;
end;

function TPerformanceCounter.GetValue: UInt64;
var
  Rslt: Int64;
begin
  if QueryPerformanceCounter(Rslt) then
    Result := Rslt
  else
    Result := 0;
end;

function TPerformanceCounter.MillisecondsBetween(Value1, Value2: UInt64): Double;
begin
  Result := (Value2 - Value1) / FFrequency;
end;

[NetSpace]

так, c корнем квадратным разорались. теперь для вычислений надо бы ещё корень кубический и любой другой степени.
а есть модуль для быстрого вычисления числа в степени?
ведь корень - это по сути число в степени 1/2.

[digitalis]

Сообщение от NetSpace

c корнем квадратным разорались

Да, орали на всю округу

а есть модуль для быстрого вычисления числа в степени?

Через логарифм. Есть что-нибудь пошустрее? Не припомню.

[Arigato]

Сообщение от digitalis

Есть что-нибудь пошустрее? Не припомню.

Однозначно есть что-то шустрее стандартного Power из модуля Math. Натыкался на его разбор давным-давно, лет так 15 назад...