Ускорить программу, которая выполняет преобразование Фурье

procedure TForm1.Button2Click(Sender: TObject);
var
q,s,i,j           : integer;
v,u,x,y,f           : integer;
g                 : boolean;
scan              : pByteArray;
img               : Timage;
md                : real;
re,im             : real;
Pixel             : TColor;
begin
 
  M := bmp1.Height;
  N := bmp1.Width;
 
    begin
    for v := 0 to M - 1 do begin
       for u := 0 to N - 1 do begin
          md := 0;  Re := 0;  Im := 0;
          irmage[v,u].Re := 0;
          irmage[v,u].Im := 0;
          for y := 0 to M - 1 do begin
            scan := bmp1.ScanLine[y];
            for x := 0 to N - 1 do begin
              f:= 0;
              f := round(scan[x]);
              if odd(x+y) then f := f*(-1);
              Re := Re + cos((-1)*2*pi*(u*x/M+v*y/N))*f;
              Im := Im + sin((-1)*2*pi*(u*x/M+v*y/N))*f;
            end;
          end;
          irmage[v,u].Re := Re/(M*N);
          irmage[v,u].Im := Im/(M*N);
          irmage2[v,u].Re := Re/(M*N);
          irmage2[v,u].Im := Im/(M*N);
          md := sqrt( Re*Re/m*n + Im*Im/m*n );
          irmage[v,u].Modul := md;
          i := round (20 * ln(1+md));
          if i >= 255 then i := 255;
          bmp2.Canvas.Pixels[u,v] := rgb( i, i, i );
       end;
    end;
    image2.Picture.Bitmap.Assign(bmp2);
    bmp1.Free;
    bmp2.Free;
    end;
end;

procedure TForm1.Button3Click(Sender: TObject);
var
  i,j,y,x,u,v :integer;
  re,im,md    :real;
  f           :integer;
begin
     bmp1:=TBitmap.Create;
     bmp2:=TBitmap.Create;
     Jpg := TJPEGImage.Create;
     Jpg.LoadFromFile(opendialog1.FileName);
     Bmp2.Assign(Jpg);
    for y:=0 to M - 1 do  begin
      for x:=0 to N - 1 do  begin
         md := 0;  Re := 0;  Im := 0; f := 1;
         for V:=0 to M - 1 do begin
           for U:=0 to N - 1 do begin
             Re := Re + irmage2[v,u].Re * COS(2*pi*(u*x/M+v*y/N)) - irmage2[v,u].Im * SIN(2*pi*(u*x/M+v*y/N));
             Im := Im + irmage2[v,u].Im * COS(2*pi*(u*x/M+v*y/N)) + irmage2[v,u].Re * SIN(2*pi*(u*x/M+v*y/N));
              end;
          end;
          if odd(x+y) then f := f*(-1);
          md := sqrt(Re*Re + Im*Im);
          i := round(re*f);
          bmp2.Canvas.Pixels[x,y] := rgb(i,i,i);
      end;
    end;
    image3.Picture.Bitmap.Assign(bmp2);
end;
end.

  M := bmp2.Height;
  N := bmp2.Width;

     Jpg.LoadFromFile(opendialog1.FileName);
     Bmp2.Assign(Jpg);

          Re := Re + cos((-1)*2*pi*(u*x/M+v*y/N))*f;
          Im := Im + sin((-1)*2*pi*(u*x/M+v*y/N))*f;

uses ... Math;
...
procedure TForm1.Button2Click(Sender: TObject);
var
...
  s1, c1: extended;
begin
...
//          Re := Re + cos((-1)*2*pi*(u*x/M+v*y/N))*f; // Вместо этих строк
//          Im := Im + sin((-1)*2*pi*(u*x/M+v*y/N))*f;
          SinCos((-1)*2*pi*(u*x/M+v*y/N), s1, c1); // теперь вот эти строки
          Re := Re + c1*f;
          Im := Im + s1*f;
...