Vlákno názorů k článku Fraktály v počítačové grafice X od Ondrej Skutka - Takové optimalizace, jako zapamatování hodnoty výpočtu x*x, které...

Takové optimalizace, jako zapamatování hodnoty výpočtu x*x, které je potřeba na dalším řádku programu přece každý kompilátor musí umět.

Jestli to ve floatech a doublech odporuje IEEE opravdu nevím. Co jsem ale chtěl říct je, že dnešní kompilery dokážou jednoduché věci optimalizovat líp, než člověk. A tak místo přemýšlení nad takovými detaily je lepší spustit gcc s parametrem -O

Mimochodem — víte, že ta neoptimalizovaná Juliova množina je rychlejší, než ta upravená?

Průměrné hodnoty (ze tří pokusů) naměřené programem time.

cx=-0,1; cy=0,7; maxiter=100000:

verze	čas
1. verze	28,5s
2. verze	30,6s
1. verze, gcc -O	13,9s
2. verze, gcc -O	16,2s

Nedalo mi to a vypocet jsem vyzkousel na Pentiu M 1,6 GHz. Parametry vypoctu:
width=800
height=600
cx=0.0
cy=0.9
maxiter=1000

Neoptimalizovana verze:
Test 1 (puvodni kod)
in: 49896 (pocet pixelu uvnitr mnoziny - pro kontrolu vypoctu)
out: 430104 (pocet pixelu vne mnoziny)
total: 480000
time 1: 69

Test 2 (upravena iteracni smycka)
in: 49896
out: 430104
total: 480000
time 2: 58

V neoptimalizovane verzi (preklad pro i386) je upravena smycka o vic nez 10 procent rychlejsi.

Optimalizovana verze (-O9 -march=pentium)
Test 1 (puvodni kod)
in: 49896
out: 430104
total: 480000
time 1: 43

Test 2 (upravena iteracni smycka)
in: 49896
out: 430104
total: 480000
time 2: 44

Po optimalizaci se casy srovnaly, dokonce je ta upravena smycka pomalejsi :-o

Optimalizovana verze (-O9 -funroll-all-loops -march=pentium)
Test 1
in: 49896
out: 430104
total: 480000
time 1: 44

Test 2
in: 49896
out: 430104
total: 480000
time 2: 44

Tady je prekladac donucen k tomu, aby smycku rozbalil a preskladal kod. Casy se (cca) vyrovnaly, asi by to chtelo vice iteraci (ale rozdil je stejne pod jednim procentem, tim se nema cenu zabyvat).


Pro pobaveni jsem kod jeste prepsal do Javy a vysledky me DOST prekvapily (spis je nekde hacek, jinak si to neumim vysvetlit):

Beh jako klient (java -client FractalSpeed)
Test 1
in: 49896
out: 430104
total: 480000
Time1: 59

Test 2
in: 49896
out: 430104
total: 480000
Time2: 54

Ty casy celkem odpovidaji tem ceckovskym (a je videt, ze u vypoctu se Java cecka dost drzi, ty kecy o pomalosti patri minulosti).

Beh jake server (java -server FractalSpeed)
Test 1
in: 49896
out: 430104
total: 480000
Time1: 1

Test 2
in: 49896
out: 430104
total: 480000
Time2: 1

Ty casy si v tomto pripade opravdu nedokazu nijak vysvetlit (nejake kesovani predchoziho behu?). Opravdu to hned, tj. po cca jedne sekunde, po spusteni skoncilo, vysledky (pocet pixelu vne a uvnitr mnoziny) jsou pritom dobre.