18 a16 (float *x
, float *y
, int *index
, int n
)
21 #pragma omp parallel for shared(x, y, index, n)
22 for (i
= 0; i
< n
; i
++)
25 x
[index
[i
]] += work1 (i
);
36 for (i
= 0; i
< 10000; i
++)
41 for (i
= 0; i
< 1000; i
++)
43 a16 (x
, y
, index
, 10000);
44 for (i
= 0; i
< 10; i
++)
45 printf ("x[%d] = %f, y[%d] = %f\n", i
, x
[i
], i
, y
[i
]);