]> git.somenet.org - pub/jan/parprog.git/blob - openmp/prefix/hillis.c
add parallel computation to hillis
[pub/jan/parprog.git] / openmp / prefix / hillis.c
1 /*
2  * O(nlog n) work algorithm (Hillis-Steele)
3  */
4 #include "hillis.h"
5
6 void algorithm (numtype x[], unsigned long size, unsigned int ops[]) {
7         numtype tmp[size];
8         unsigned long k;
9         unsigned long i;
10
11         for(k=2; k <= size; k <<=1){
12                 #pragma omp parallel for shared(x, size, ops, k) private(i)
13                 for(i = (k-1); i < size; i+=k){
14                         //printf ("x[%2li] = x[%2li] + x[%2li]; // {i:%li, k:%li}\n", i, i-k, i, i, k);
15                         x[i] = x[i-(k/2)] + x[i];
16                 }
17                 //arrayState(x, size);
18                 //barrier;
19         }
20
21         /* Copy out array to in */
22 //      for (i = 0; i < size; i++){
23 //              x[i] = tmp[i+1];
24 //      }
25 }