1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
#include <stdio.h>
#include <string.h>
#define ARRAY_LEN 4096
double x[ ARRAY_LEN ];
double y[ ARRAY_LEN ];
int main( int argc, char **argv )
{
const double alpha = 2.1;
for ( int i = 0; i < ARRAY_LEN; ++i )
{
y[ i ] = 0.0;
x[ i ] = 1.0;
}
#pragma omp simd
for ( int i = 0; i < ARRAY_LEN; ++i )
{
y[ i ] += alpha * x[ i ];
}
#if 0
_Pragma( "omp simd" )
for ( int i = 0; i < ARRAY_LEN; ++i )
{
x[ i ] += alpha * y[ i ];
}
#endif
return 0;
}
|
