оптимизировать для максимальной скорости выполнения программы, т.е. заменить ее кодом на ассемблере

[gusluk]

Есть процедура расчета некой функции. На делфи она выглядит так:

for i:=1 to N do
begin
     for j:=1 to M do
     begin
     if (Fx[i]>KrMin) and (Fx[i]<KrMax) and ((Fx[i]<-MinF)or (Fx[i]>MinF)) then    
     Fx[j]:=Fx[j]+ dx[i,j]*Fx[i];
     end;
Fx[i]:=0;
end;

Хочу ее оптимизировать для максимальной скорости выполнения программы, т.е. заменить ее кодом написанным на ассемблере, для этого начал потихоньку его изучать. Но есть вопрос: а есть ли в этом смысл, реально ли добиться серьезного увеличения скорости выполнения этой функции на ассемблере по сравнению с Delphi 7? И если реально что можете посоветовать по этому поводу?

[airyashov]

Для начала можно посмотреть асм код который делфи сделала.
В данном случае, проще разбить цикл на 4(2) потока, 4(2) ядра в помощь, никая оптимизация на ассемблере не даст выигрыш 3/4 времени из-за сложной адресации элементов массивов, хотя процентов 15-30% можно получить переписав все на fpu вручную.

[alexBlack]

Может сначала попробовать оптимизировать в Delphi. Например:

   for i:=1 to N do begin
      // условие, похоже зависит только от i (кроме одного случая)
      // если условие не выполняется нет смысла выполнять внутренний цикл
      B := (Fx[i]>KrMin) and (Fx[i]<KrMax) and ((Fx[i]<-MinF)or (Fx[i]>MinF));
      if B then begin
         for j:=1 to M do begin

            Fx[j]:=Fx[j]+ dx[i,j]*Fx[i];

            if j = i then begin
               // А здесь условие может поменяться
               // вычисляем его заново
               B := (Fx[i]>KrMin) and (Fx[i]<KrMax) and ((Fx[i]<-MinF)or (Fx[i]>MinF));
               // И опять если условие не выполняется нет смысла продолжать
               if not B then break;
            end;

         end;
      end;
      Fx[i]:=0;
   end;

Конечно, если KrMin, KrMax, MinF - это константы, а не функции с побочными эффектами.

[gusluk]

К сожалению асм код делфи для этой функции пока показать не могу, программа еще не дописана, но есть почти точно такая же функция которую я делал для проверки алгоритма, и позже его попробую выложить.
airyashov, что ты имел в виду когда сказал что проще разбить цикл на 4(2) потока? Я просто пихал этот цикл в два потока, т.к. машина двух ядерная, на никак его не разбивал.
alexBlack, спасибо, я тоже уже думал оптимизировать его на делфи, но сначала хотел узнать на сколько его можно оптимизировать на асм.
KrMin, KrMax, MinF задаются с исходными данными, и в процессе вычислений не меняются.

[gusluk]

Проанализировал асм код делфи, пришел к выводу что больше чем на 20% его не ускорить. В итоге решил его сначало оптимизировать как предложил
alexBlack. Посмотреть как его с компилирует делфи,и если будет возможность увеличить его производительность более чем на 30% то переписать на ассемблере, и пихнуть в два потока. Как вы считаете, это наиболее подходящий вариант, или можно что то получше придумать?

[gusluk]

Вот пример функции в делфи:

for i:=1 to 90 do
begin
     for j:=1 to 100 do
     begin
     if (Ln[i]>-5000) and (Ln[i]<5000) and ((Ln[i]<-1)or  (Ln[i]>1)) then Ln[j]:=Ln[j]+ Lsin[i,j]*Ln[i];

     end;
     Ln[i]:=0;
end;

и асм:

for i:=1 to 90 do
mov [ecx],$00000001
for j:=1 to 100 do
mov [ebx],$00000001
 if (Ln[i]>-5000) and (Ln[i]<5000) and ((Ln[i]<-1)or  (Ln[i]>1)) then Ln[j]:=Ln[j]+ Lsin[i,j]*Ln[i];
mov eax,[ecx]
fld qword ptr [$004523d0]
fstsw ax
sahf
jbe +$58
mov eax,[ecx]
fld qword ptr [edx+eax*8-$08]
fcomp dword ptr [$004523d4]
fstsw ax
sahf
jnb +$47
mov eax,[ecx]
fld qword ptr [edx+eax*8-$08]
fcomp dword ptr [$004523d8]
fstsw ax
sahf
jb +$11
mov eax,[ecx]
fld qword ptr [edx+eax*8-$08]
fcomp dword ptr [$004523dc]
fstsw ax
sahf
jbe +$25
imul eax,[ecx],$000003e8
lea eax,[esi+eax*8]
mov edi,[ebx]
fld qword ptr [eax+edi*8-$00001f48]
mov eax,[ecx]
fmul qword ptr [edx+eax*8-$08]
mov eax,[ebx]
fadd qword ptr [edx+eax*8-$08]
mov eax,[ebx]
fstp qword ptr [edx+eax*8-$08]
wait
end;
inc dword ptr [ebx]
for j:=1 to 100 do
cmp dword ptr [ebx],$65
jnz -$70
Ln[i]:=0;
mov eax,[ecx]
xor edi,edi
mov [edx+eax*8-$08],edi
mov [edx+eax*8-$04],edi
end;
inc dword ptr [ecx]
for i:=1 to 90 do
cmp dword ptr [ecx],$5b
jnz -$0000008d
end;
pop edi
pop esi
pop ebx
ret
add [eax],al
inc eax
pushfd
ids eax, [eax]
inc eax
pushfd
inc ebp
add [eax],al
cmp byte ptr [edi+$3f800000],$53
push esi
push edi
mov edx, $004590d0
mov ecx, $0139d4d8
mov ebx, $0139d4dc
mov esi,$00bfc210

[gusluk]

Странно, но вариант предложенный alexBlack дал выигрыш по времени около 20% , я ожидал большего. Попозже попробую выложить его асм код.