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

program Fractal;
uses windows;
type
BMP_HEADER=record
bfType:WORD; {BM}
FileSize:DWORD;
NoneUse:DWORD;
BytesToData:DWORD;
end;
BitMapInfo=record
StructSize:DWORD;
bmpWidth:DWORD;
bmpHeight:DWORD;
Planes:WORD;
BitPerPixel:WORD;
TypeCompression:DWORD;
DataSize:DWORD;
XPelsPerMeter:DWORD;
YPelsPerMeter:DWORD;
ColorUsed:DWORD;
ColorImportant:DWORD;
end;
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;
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
Zagolovok:BMP_HEADER;
Opis :BitMapInfo;
Pixels :^PixEkr;
Palette :array [1..256] of PixelType;
f1:file;
T,x,y:DWORD;
Odin,C,Z,A,B,D,E,F:Complex;
minx,miny,maxx,maxy,tx:single;
n:DWORD;
Begin
GetMem(Pixels,3*512*512);
Odin.Re:= 1.0;
Odin.Im:=-1.0;
B.Re:=1;
B.Im:=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;
for x:=1 to 512 do {}
begin
for y:=1 to 512 do
begin
C.Re:=minx+(maxx-minx)*x/512.0;
C.Im:=miny+(maxy-miny)*y/512.0;
Z:=C;
n:=0;
repeat
asm {C=A*B}
femms




movq mm0,Z {mm0=Z.Re,Z.Im}
pswapd mm1,Z {mm1=Z.Im,Z.Re}
pfmul mm1,mm0 {mm1=Z.Im*Z.Re,Z.Re*Z.Im}
pfmul mm0,Z {mm0=Z.Re*Z.Re,Z.Im*Z.Im}
pfpnacc mm0,mm1 {mm0=Z.Re*Z.Re-Z.Im*Z.Im,
Z.Im*Z.Re+Z.Re*Z.Im}
movq A,mm0 {збереження результату}


movq mm0,A {mm0=A.Re,A.Im}
pswapd mm1,Z {mm1=B.Im,B.Re}
pfmul mm1,mm0 {mm1=B.Im*A.Re,B.Re*A.Im}
pfmul mm0,Z {mm0=A.Re*B.Re,A.Im*B.Im}
pfpnacc mm0,mm1 {mm0=A.Re*B.Re-A.Im*B.Im,
B.Im*A.Re+B.Re*A.Im}
movq D,mm0 {збереження результату}


movq mm0,C {mm0={A.Re,A.Im}
pfsub mm0,D {mm0={A.Re-B.Re,A.Im-B.Im}
movq E,mm0 {збереження результату}


movq mm0,E {mm0={A.Re,A.Im}
pfadd mm0,B {mm0={A.Re+B.Re,A.Im+B.Im}
movq Z,mm0 {збереження результату}


pfmul mm0,mm0
pfacc mm0,mm0
movd tx,mm0




femms
end;
n:=n+1;
until( (tx>4.0) or (n=255) );{}
Pixels^[y,x]:=Palette[n];
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);
readln;
End.