(转)DELPHI求乱序数组最小值快速算法(Delphi MMX优化算法应用之三)

type

  PMyList = ^TMylist;

  TMylist = array of word;



function GetMinValue(List: PMyList):Word; stdcall;

var

  ForCount, ListCount,i: Dword;

  MAXVaule: int64;

  List0Address: Dword;



begin

  Result := $7FFF;

  ListCount := High(List^) + 1;

  if (ListCount < 16) or ((ListCount mod 16)<>0) then

  begin

   for I := 0 to ListCount - 1 do

      begin

         if List^[i] <  Result then  Result := List^[i];

      end;

  end else

  begin

    MAXVaule := $7FFF7FFF7FFF7FFF;

    List0Address := Dword(@List^[0]);

    ForCount := ListCount div 16; //4*MMX

    asm

      push eax

      push ecx

      push esi

      mov esi, List0Address  //开始地址

      mov ecx, ForCount      //长度

      MOVQ MM0,MAXVaule

      mov eax,0

    @CmpLoop:

      prefetchnta [esi + eax + 1024] // fetch ahead by 1024 bytes

      movq mm1, qword [esi+eax]

      movq mm2, qword [esi+eax+8]

      movq mm3, qword [esi+eax+16]

      movq mm4, qword [esi+eax+24]

      PMINSW MM0,MM1

      PMINSW MM0,MM2

      PMINSW MM0,MM3

      PMINSW MM0,MM4

      add eax,32

      sub ecx, 1

    jnz @CmpLoop

      MOVQ MM1,MM0

      PSRLQ MM1,$20  //右移

      PMINSW MM0,MM1

      MOVQ MM1,MM0

      PSRLQ MM1,$10  //右移

      PMINSW MM0,MM1

      Movd eax,mm0

      Mov  Result,AX

      Sfence

      Emms

      pop esi

      pop ecx

      pop eax

    end;

  end;



end;