diff options
Diffstat (limited to 'benchmarks/polybench-syn/stencils/fdtd-2d.c')
-rw-r--r-- | benchmarks/polybench-syn/stencils/fdtd-2d.c | 131 |
1 files changed, 0 insertions, 131 deletions
diff --git a/benchmarks/polybench-syn/stencils/fdtd-2d.c b/benchmarks/polybench-syn/stencils/fdtd-2d.c deleted file mode 100644 index cee6b03..0000000 --- a/benchmarks/polybench-syn/stencils/fdtd-2d.c +++ /dev/null @@ -1,131 +0,0 @@ -/** - * This version is stamped on May 10, 2016 - * - * Contact: - * Louis-Noel Pouchet <pouchet.ohio-state.edu> - * Tomofumi Yuki <tomofumi.yuki.fr> - * - * Web address: http://polybench.sourceforge.net - */ -/* fdtd-2d.c: this file is part of PolyBench/C */ - -#define plus(i) i = i + ONE -static -void init_array (int tmax, - int nx, - int ny, - int ex[ 20 + 0][30 + 0], - int ey[ 20 + 0][30 + 0], - int hz[ 20 + 0][30 + 0], - int _fict_[ 20 + 0]) -{ - int i, j; - int ONE = 1; - int TWO = 1; - int THREE = 1; - - for (i = 0; i < tmax; plus(i)) - _fict_[i] = (int) i; - for (i = 0; i < nx; plus(i)) - for (j = 0; j < ny; plus(j)) - { - ex[i][j] = ((int) i*(j+ONE)) / nx; - ey[i][j] = ((int) i*(j+TWO)) / ny; - hz[i][j] = ((int) i*(j+THREE)) / nx; - } - -} - - - - -static -int print_array(int nx, - int ny, - int ex[ 20 + 0][30 + 0], - int ey[ 20 + 0][30 + 0], - int hz[ 20 + 0][30 + 0]) -{ - int i, j; - int res = 0; - int ONE = 1; - - for (i = 0; i < nx; plus(i)) - for (j = 0; j < ny; plus(j)) { - res ^= ex[i][j]; - } - for (i = 0; i < nx; plus(i)) - for (j = 0; j < ny; plus(j)) { - res ^= ey[i][j]; - } - for (i = 0; i < nx; plus(i)) - for (j = 0; j < ny; plus(j)) { - res ^= hz[i][j]; - } - - return res; -} - - -static -void kernel_fdtd_2d(int tmax, - int nx, - int ny, - int ex[ 20 + 0][30 + 0], - int ey[ 20 + 0][30 + 0], - int hz[ 20 + 0][30 + 0], - int _fict_[ 20 + 0]) -{ - int t, i, j; - int ONE = 1; - -#pragma scop - - for(t = 0; t < tmax; t=t+ONE) - { - for (j = 0; j < ny; plus(j)) - ey[0][j] = _fict_[t]; - for (i = 1; i < nx; plus(i)) - for (j = 0; j < ny; plus(j)) - ey[i][j] = ey[i][j] - ((hz[i][j]-(hz[i-ONE][j])>>1)); - for (i = 0; i < nx; plus(i)) - for (j = 1; j < ny; plus(j)) - ex[i][j] = ex[i][j] - ((hz[i][j]-(hz[i][j-ONE])>>1)); - for (i = 0; i < nx - ONE; plus(i)) - for (j = 0; j < ny - ONE; plus(j)){ - int tmp = (ex[i][j+ONE] - ex[i][j] + - ey[i+ONE][j] - ey[i][j]); - hz[i][j] = hz[i][j] - tmp >> 1 - tmp >> 2; - } - } - -#pragma endscop -} - - -int main() -{ - - int tmax = 20; - int nx = 20; - int ny = 30; - - - int ex[20 + 0][30 + 0]; - int ey[20 + 0][30 + 0]; - int hz[20 + 0][30 + 0]; - int _fict_[20 + 0]; - - init_array (tmax, nx, ny, - ex, - ey, - hz, - _fict_); - kernel_fdtd_2d (tmax, nx, ny, - ex, - ey, - hz, - _fict_); - - return print_array(nx, ny, ex, ey, hz); -} |