Dotaz: efektivita procházení matice

[Dadam]

Zdravím,
mám 32bitový Arch linux, procesor Core2Duo, gcc 4.5.0, optimalizace vypnuté, měřím pomocí gprof. Je nějaký důvod, proč by měl být následující kód

        for (i = 0; i < in->h; i++) {
                for (j = 0; j < in->w; j++) {
                        b = 0;
                        for (k = -r; k <= r; k++) for (l = -r; l <= r; l++) { 
                                if (in->arr[i+l][j+k] <= p) b++;
                        }
                        if (b > (max-1)/2) out->arr[i][j] = 0;
                        else out->arr[i][j] = 255;
                }
        }

pomalejší než tento:

        for (i = 0; i < in->h; i++) {
                for (j = 0; j < in->w; j++) {
                        b = 0;
                        for (k = 0; k <= 2*r; k++) for (l = 0; l <= 2*r; l++) {
                                if (in->arr[i-r+l][j-r+k] <= p) b++;
                        }
                        if (b > (max-1)/2) out->arr[i][j] = 0;
                        else out->arr[i][j] = 255;
                }
        }

? Rozdíl je jen pár procent, nicméně myslel jsem že bude na druhou stranu, nebo nebude žádný.

Odpovědi

[Dadam]

Mimochodem - že se to ocitne v záporných indexech vím, je to sice prasárna, nicméně počítám s tím.

[reqw5q23]

Zkus si vytknout 2*r a vyrazy v indexech poli si take predpocitej.

[heretic]

Neměly by bejt u toho druhýho příkladu indexy náhodou takhle ?

if (b > (max-1)/2) out->arr[i-r][j-r] = 0;
else out->arr[i-r][j-r] = 255;

[Dadam]

Ee, v in se prochází matice šířky 2r+1 kolem indexu [i][j], výsledek se zapíše do out indexu [i][j] (střed procházené matice). Je to mediánový filtr na černobílé bitmapě.

[pht]

Obvykle jsou procesorům milejší cykly s jedním koncem v nule. Podívejte se co z toho leze za assembler (gcc -S) a třeba to bude jasnější.

[Dadam]

Podívám. S assemblerem moc nepracuju, že existuje nějaká instrukce loop, která končí na nule, vím, ale ta rychlejší varianta na nule začíná, a vzhledem k tomu, že mám vypnuté optimalizace, počítám s tím, že mi to celé nepřeskládá.