Discussion :

[Arzar]

Alors (avec N=100 cette fois, j'avais mis à 64 par erreur plus haut)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
Test original :
3520.54 ms
Test avec float* Q = new float[N*N*N+16]+16;
1919.73 ms
Test avec le dernier code de screetch (alloca)
2018.73

Le static de static void compute... ne semble pas avoir d'effet particulier.

Donc au final, si je comprends bien, le problème venait du fait que le compilateur aligne en mémoire les deux adresses de P et Q (lors du new), ce qui fait qu'il ne peut plus y accéder simultanément après (mais alors pourquoi les aligner ?). Et dans le dernier code les calculs intermédiaires sont fait dans la pile (avec alloca) pour forcer le compilateur a comprendre qu'il n' y a pas d'aliasing, le tas et la pile étant bien distinct. C'est à peu près ça ?

[screetch]

Mathieu, le principal est de faire un bout de calcul, le mettre dans un buffer, copier le buffer, recommencer
je le fais sur une ligne, mais ca pourrait etre un nombre arbitraire, genre 128, qui ne coïncide pas avec les lignes

Arzar, le probleme c,est qu'en informatique beaucoup de choses reposent sur des puissances de 2, et il y a des chances que l'allocateur alloue au meme endroit (enfin, aligné), meme si tu lui demandes 1500
en faisant

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
	float* P = new float[N*N*N*2];
	float* Q = &P[N*N*N];

si N n'est pas une puissance de 2, ca marche
si N est une puissance de 2, kaboom, les perf meurent.

Sur la pile, tu as moins de chance d'etre aligné mais tu peux quand meme l'etre (ca serait la faute a pas de bol...)
l'idéal donc :

- avoir un buffer sur la pile (peu importe la taille, mais un truc qui permet de travailler en batch)
- désaligner de manière certaine le buffer et les paramètres d'entrée/sortie en sautant les premiers éléments (un ptit peu de math avec l'adresse des buffers...)
- faire le calcul et stocker dans le buffer
- lorsque le buffer est plein, on le copie dans le résultat

[loufoque]

Et en utilisant posix_memalign au lieu de new[] ?

[Matthieu Brucher]

Mathieu, le principal est de faire un bout de calcul, le mettre dans un buffer, copier le buffer, recommencer
je le fais sur une ligne, mais ca pourrait etre un nombre arbitraire, genre 128, qui ne coïncide pas avec les lignes

A voir, je note dans la liste des possibilités. Je vais déjà tester si le compilateur me permet d'avoir de bonnes perfs en mettant le stencil dans une sous-fonction et passer par structure.

edit: bon, il faut donner explicitement le stencil dans la boucle, sinon ça ne fonctionne pas