diff options
Diffstat (limited to 'benchmarks/polybench-syn/stencils/adi.c')
| -rw-r--r-- | benchmarks/polybench-syn/stencils/adi.c | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/benchmarks/polybench-syn/stencils/adi.c b/benchmarks/polybench-syn/stencils/adi.c index cc54e11..5dd406b 100644 --- a/benchmarks/polybench-syn/stencils/adi.c +++ b/benchmarks/polybench-syn/stencils/adi.c @@ -70,13 +70,14 @@ void kernel_adi(int tsteps, int n, d = -mul2 / 2; e = 1+mul2; f = d; + int ZERO = 0; for (t=1; t<=tsteps; t++) { for (i=1; i<n-1; i++) { - v[0][i] = 1; - p[i][0] = 0; - q[i][0] = v[0][i]; + v[ZERO][i] = 1; + p[i][ZERO] = 0; + q[i][ZERO] = v[ZERO][i]; for (j=1; j<n-1; j++) { p[i][j] = -c / (a*p[i][j-1]+b); q[i][j] = (-d*u[j][i-1]+(1+2*d)*u[j][i] - f*u[j][i+1]-a*q[i][j-1])/(a*p[i][j-1]+b); @@ -89,9 +90,9 @@ void kernel_adi(int tsteps, int n, } for (i=1; i<n-1; i++) { - u[i][0] = 1; - p[i][0] = 0; - q[i][0] = u[i][0]; + u[i][ZERO] = 1; + p[i][ZERO] = 0; + q[i][ZERO] = u[i][ZERO]; for (j=1; j<n-1; j++) { p[i][j] = -f / (d*p[i][j-1]+e); q[i][j] = (-a*v[i-1][j]+(1+2*a)*v[i][j] - c*v[i+1][j]-d*q[i][j-1])/(d*p[i][j-1]+e); |