[C++ & asm] najhitrejša inicializacija 2D matrike

Mene zanima kako najhitreje postavit vse elemente v matriki na 0?

Zaenkrat mam trotl ziher varianto:

for(int x=0;x<iNodes;x++)
{
    for(int y=0;y<iNodes;y++)
    {
        conntaken[x][y]=0;
    }
}

Bi šlo z memset?

Kaj pa recimo sse? A se tam 128 bitov enako hitro prenese kot če prenašaš 32 bitne podatke iz registra v mem?

Za dynamic memory ne moreš nič rečt kake vrednosti bojo notri po defaultu, tako da ko nekaj alociraš vsebina ni nujno 0, ne vem od kje ti taka ideja BigWhale.

Če imaš matriko kot statičen C++ array potem lako poskusiš nekaj takega:
int conntaken[10][10] = {0};
in pustiš prevajalniku da se spomni čimhitrejšega načina da to izvede. Ampak to gre le ob initializaciji, ko definiraš matriko. Zelo verjetno bi pa res bilo hitreje če uporabiš malo assemblerja. Mislim da se celo dela tako, da nastaviš en 128biten register na 0 in ga potem kopiraš v RAM. Morda celo obstaja 128 bitni STOS ukaz. Mislim da se ti splača pogledat, sem nekje bral da se to obnese (ali pa to morda velja le za kopiranje podatkov, ne vem točno).

Aha sem pogledal mal po manualu in ugotovil da gre takole:
PXOR xmm1,xmm1
MOVDQA mem,xmm1

PS: To je javendar bistvena razlika med malloc() in calloc()

int conntaken[Y][X];
    
int *p;
p = &conntaken[0][0];

for(i = 0; i < Y * X; i++)
      *(p + i) = 0;

Na koncu se ti to najverjetneje ne bo poznal (razen ce res skos sam inicializiras matrike 90% cajta) in se ne splaca tok pacat po kodi za zlo mal benefita... al lepo vse memsetas (ne bos verjel kaj dober kompajler zna delat s tem) al pa naredis se en byte (al pa int ce hoces met alignan) "identity" in ga das na 1. Pa ga upostevas pri vseh kalkulacijah... in bos se tam ful profitiral :) Tko smo mel mi neki cajta v enginu sam mislm, da zdej se tega nimamo vec :)

Hvala vsem za nasvete!

Predvsem mi je zanimiva tista od Gundolfa z MOVNTDQ—Store Double Quadword Using Non-Temporal Hint.


@Fury
Vem da se mi zadeva zelo velikokrat izvede... zdej kolk % časa to vzame ne vem.. sem mi pa splača vsako zadevo preštudirat.

Tale zbirka šablon bi znala biti zanimiva:
http://www.pixelglow.com/macstl/

edit: typo

Mene pa zanima - kakšne vrednosti shranjuješ v matriko? Koliko bitne?

No prov. Povedal nisi, kar sem te vprašal, zato bom vsaj jaz poskušal odgovoriti na tisto, kar si vprašal ti.

Čeprav, približno sem že zgoraj!

Najhitrejši način za reinicializacijo matrike je naslednji. Sploh je ne reinicializiraš vsakič! Samo vsake toliko.

Zdej pa denimo, da shranjuješ v matriko 10 bitne vrednosti. Ob prvem loopu postavljaš notri vrednosti med 0 in 1023. V drugem loopu med 1024 in 2047 ... itd?

Ker tudi tam v N-tem loopu, kjer testiraš, vzameš vse vrednosti pod N*1024 za 0! Med (N+1)*1024 pa za število med 0 in 1023.

Kapiš?

No, če kapiš, potem boš pri 32 bitnih možnih vrednostih matrike redimenzioniral samo vsakih MILIJON loopov!

Kapiš?

Ja, koda, ki to enostavno dela, je pa seveda MMX (ali GWXH ...) in je ne bom (ne smem) napisal. Dokaj lahko jo pa (eno) napišeš sam.

Se razumemo?

So to slučajno razpršene matrike?

> So to slučajno razpršene matrike?

Za veliko večino teh zadev to vse lepo prime. So pa možne matrike in primeri (vsaj extra designed) kjer se ne splača.

No sej. Po moje da ma in da se mu bo kar splačalo.

Dej uporab memset k se ti ne bo nc poznal.. s kaksnimi stvarmi se vi sekirate :)

> Ja, ne rabi initializacije,

Super, ane? Sej ne bi verjel da je možno, če ne bi vidu?

> pri premikanju, and-anju itd. mora pa še vedno vse elemente premaknit, and-at, ...

Lahko se zgodi da rabiš, lahko se zgodi da ne rabiš (vsega) tega.

Lahko se zgodi, da je potem še bolj pripravno.

> Torej boš še vedno reinicializiral na vsake par milijardokrat!

Me happy.