+__global__ void genere_diagos_rectangle(uint4 * d_diagos, int h, int l, int q){
+ int inci = h/q;
+ int incj = l/q;
+ int iM,jM, iN, jN ;
+ int idxDiago = 0 ;
+ // boucle double pour les positions du point NO de la diagonale
+ for ( iM = 0; iM < q-1; iM++){
+ for (jM = 0 ; jM < q-1 ; jM++){
+ //boucle double pour les positions du point SE de la diagonale
+ for (iN = iM+1; iN < q; iN++){
+ for (jN = jM+1; jN < q; jN++){
+ d_diagos[idxDiago].x = iM*inci;
+ d_diagos[idxDiago].y = jM*incj;
+ d_diagos[idxDiago].z = iN*inci;
+ d_diagos[idxDiago].w = jN*incj;
+ idxDiago++;
+ }
+ }
+ }
+ }
+}