готов лаби 4 курс 1 сим / цеом / lipa44

program Fractal;
uses crt,windows;

  PixelType=record
    B,G,R:BYTE;
  end;
  Ekran =array [0..512*512-1] of PixelType;
  PixEkr=array [1..512,1..512] of PixelType;
  Complex=record
    Re:single;
    Im:single;
  end;
  DComplex=record
    C1:Complex;
    C2:Complex;
  end;
  ParaWord=record
    n1,n2:WORD;
  end;
const
  Shapka:array [1..54] of byte=($42,$4D,$36,0,$C,0,0,0,
                                0,0,$36,0,0,0,$28,0,
                                0,0,0,2,0,0,0,2,0,0,1,0,
                                $18,0,0,0,0,0,0,0,$C,0,
                                $12,$B,0,0,$12,$B,0,0,0,
                                0,0,0,0,0,0,0);
var
  Pixels   :^PixEkr;
  Palette  :array [1..256] of PixelType;

  f1:file;
  T,T1:int64;//QWORD;
  max,min,tx,Odin,C,Z:DComplex;
  minx,miny,maxx,maxy:single;
  n,x,y:DWORD;
  Kol:ParaWord;

Begin
  clrscr;
  GetMem(Pixels,3*512*512);
  Odin.C1.Re:= 1.0;
  Odin.C1.Im:= 1.0;
  Odin.C2.Re:=-1.0;
  Odin.C2.Im:=-1.0;
  max.C1.Re:= 4.0;
  max.C1.Im:= 4.0;
  max.C2.Re:= 4.0;
  max.C2.Im:= 4.0;
  min.C1.Re:=-4.0;
  min.C1.Im:=-4.0;
  min.C2.Re:=-4.0;
  min.C2.Im:=-4.0;

  minx:=-1.5;
  maxx:= 0.5;
  miny:=-1.0;
  maxy:= 1.0;
  for n:=1 to 256 do
  begin
    Palette[n].R:=255-abs(127-n);
    Palette[n].B:=256-n;
    Palette[n].G:=n-1;
  end;
  T:=GetTickCount;
      asm
        movaps  xmm7,Odin
        movaps  xmm6,max
        movaps  xmm5,min
      end;

  for x:=1 to 512 do 
  begin
    for y:=0 to 255 do
    begin
      C.C1.Re:=minx+(maxx-minx)*x/512.0;
      C.C1.Im:=miny+(maxy-miny)*(y*2.0)/512.0;
      C.C2.Re:=minx+(maxx-minx)*x/512.0;
      C.C2.Im:=miny+(maxy-miny)*(y*2.0+1.0)/512.0;
      Z:=C;
      KOL.N1:=1;
      KOL.N2:=1;
      ASM MOVAPS XMM0,Z END;
      REPEAT
        ASM
         
	movaps  xmm3,XMM0    {xmm0=y1,x1,y2,x2}		{Возводим в 4 степень - В}
	movaps  xmm1,XMM0    {xmm1=y3,x3,y4,x4}
	movaps  xmm2,xmm0 {xmm2=y1,x1,y2,x2}
	mulps   xmm2,xmm1 {xmm2=y1*y3,x1*x3,y2*y4,x2*x4}
	shufps  xmm1,xmm1,10110001b {xmm1=x3,y3,x4,y4}
	mulps   xmm1,xmm0 {xmm1=y1*x3,x1*y3,y2*x4,x2*y4}
	movaps  xmm0,xmm2 {xmm0=y1*y3,x1*x3,y2*y4,x2*x4}
	shufps  xmm0,xmm1,11011000b {xmm0=y1*y3,y2*y4,y1*x3,y2*x4}
	shufps  xmm2,xmm1,10001101b{xmm2= x1*x3, x2*x4,x1*y3,x2*y4}
	mulps   xmm0,Odin          {xmm0=-y1*y3,-y2*y4,y1*x3,y2*x4}
	addps   xmm0,xmm2 {xmm0=X1,X2,Y1,Y2}
	shufps  xmm0,xmm0,11011000b  {xmm0=Y1,X1,Y2,X2}
	movaps  A,xmm0


	movaps  xmm0,A    {xmm0=y1,x1,y2,x2}		{Возводим в 4 степень - В}
	movaps  xmm1,A    {xmm1=y3,x3,y4,x4}
	movaps  xmm2,xmm0 {xmm2=y1,x1,y2,x2}
	mulps   xmm2,xmm1 {xmm2=y1*y3,x1*x3,y2*y4,x2*x4}
	shufps  xmm1,xmm1,10110001b {xmm1=x3,y3,x4,y4}
	mulps   xmm1,xmm0 {xmm1=y1*x3,x1*y3,y2*x4,x2*y4}
	movaps  xmm0,xmm2 {xmm0=y1*y3,x1*x3,y2*y4,x2*x4}
	shufps  xmm0,xmm1,11011000b {xmm0=y1*y3,y2*y4,y1*x3,y2*x4}
	shufps  xmm2,xmm1,10001101b{xmm2= x1*x3, x2*x4,x1*y3,x2*y4}
	mulps   xmm0,Odin          {xmm0=-y1*y3,-y2*y4,y1*x3,y2*x4}
	addps   xmm0,xmm2 {xmm0=X1,X2,Y1,Y2}
	shufps  xmm0,xmm0,11011000b  {xmm0=Y1,X1,Y2,X2}
	movaps  B,xmm0


	movaps  xmm0,A    {xmm0=y1,x1,y2,x2}		{Возводим в 4 степень - В}
	movaps  xmm1,Z  {xmm1=y3,x3,y4,x4}
	movaps  xmm2,xmm0 {xmm2=y1,x1,y2,x2}
	mulps   xmm2,xmm1 {xmm2=y1*y3,x1*x3,y2*y4,x2*x4}
	shufps  xmm1,xmm1,10110001b {xmm1=x3,y3,x4,y4}
	mulps   xmm1,xmm0 {xmm1=y1*x3,x1*y3,y2*x4,x2*y4}
	movaps  xmm0,xmm2 {xmm0=y1*y3,x1*x3,y2*y4,x2*x4}
	shufps  xmm0,xmm1,11011000b {xmm0=y1*y3,y2*y4,y1*x3,y2*x4}
	shufps  xmm2,xmm1,10001101b{xmm2= x1*x3, x2*x4,x1*y3,x2*y4}
	mulps   xmm0,Odin          {xmm0=-y1*y3,-y2*y4,y1*x3,y2*x4}
	addps   xmm0,xmm2 {xmm0=X1,X2,Y1,Y2}
	shufps  xmm0,xmm0,11011000b  {xmm0=Y1,X1,Y2,X2}
	movaps  D,xmm0


		
	movaps xmm0,D
    	addps  xmm0,B
   	movaps Z,xmm0

	

	movaps xmm0,A
    	addps  xmm0,Z
   	movaps Z,xmm0


	movaps xmm0,Z
    	addps  xmm0,C
   	movaps Z,xmm0


	  MAXPS  XMM0,XMM5 
          MINPS  XMM0,XMM6 
          MOVAPS XMM3,XMM0
          MULPS  XMM3,XMM3 
          MOVAPS XMM1,XMM3
          SHUFPS XMM1,XMM1,10110001B
          ADDPS  XMM3,XMM1
          MOVAPS TX  ,XMM3


        END;
        IF (TX.C1.RE<4.0) THEN
          KOL.N1:=KOL.N1+1;
        IF (TX.C2.RE<4.0) THEN
          KOL.N2:=KOL.N2+1;
UNTIL( ((TX.C1.RE>4.0)AND(TX.C2.RE>4.0)) OR (KOL.N1=255) OR (KOL.N2=255));
      PIXELS^[Y*2+1,X]:=PALETTE[KOL.N1];
      PIXELS^[Y*2+2,X]:=PALETTE[KOL.N2];
    END;
  END;
  T:=GETTICKCOUNT-T;
  WRITELN('TIME CALCKULATE IS T=',T,' MS.');
  assign(f1,'Fractal.bmp');
  rewrite(f1,1);
  blockwrite(f1,Shapka,54);
  blockwrite(f1,Pixels^,3*512*512);
  FreeMem(Pixels,3*512*512);
  close(f1);
End.