bones-compiler 1.1.0 → 1.3.1

data/test/examples/benchmarks/PolyBench/adi_species.c

@@ -0,0 +1,89 @@
+//
+// This file is part of the Bones source-to-source compiler examples. The C-code
+// is largely identical in terms of functionality and variable naming to the code
+// found in PolyBench/C version 3.2. For more information on PolyBench/C or Bones
+// please use the contact information below.
+//
+// == More information on PolyBench/C
+// Contact............Louis-Noel Pouchet <pouchet@cse.ohio-state.edu>
+// Web address........http://polybench.sourceforge.net/
+// 
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........benchmark/adi.c
+// Author.............Cedric Nugteren
+// Last modified on...29-May-2012
+//
+
+#include "common.h"
+
+// This is 'adi', an alternating direction implicit solver
+int main(void) {
+	int t,i,j,i1,i2;
+	
+	// Declare arrays on the stack
+	float X[N][N];
+	float A[N][N];
+	float B[N][N];
+	
+	// Set the input data
+	for (i=0; i<N; i++) {
+		for (j=0; j<N; j++) {
+			X[i][j] = ((float) i*(j+1) + 1) / N;
+			A[i][j] = ((float) i*(j+2) + 2) / N;
+			B[i][j] = ((float) i*(j+3) + 3) / N;
+		}
+	}
+	
+	// Perform the computation
+	#pragma scop
+	{
+		for (t = 0; t < TSTEPS; t++) {
+			for (i1 = 0; i1 < N; i1++) {
+				for (i2 = 1; i2 < N; i2++) {
+					X[i1][i2] = X[i1][i2] - X[i1][i2 - 1] * A[i1][i2] / B[i1][i2 - 1];
+					B[i1][i2] = B[i1][i2] - A[i1][i2] * A[i1][i2] / B[i1][i2 - 1];
+				}
+			}
+			#pragma species kernel X[0:N-1,N-1:N-1]|element ^ B[0:N-1,N-1:N-1]|element -> X[0:N-1,N-1:N-1]|element
+			for (i1 = 0; i1 < N; i1++) {
+				X[i1][N - 1] = X[i1][N - 1] / B[i1][N - 1];
+			}
+			#pragma species endkernel adi_k3
+			for (i1 = 0; i1 < N; i1++) {
+				for (i2 = 0; i2 < N - 2; i2++) {
+					X[i1][N - i2 - 2] = (X[i1][N - 2 - i2] - X[i1][N - 2 - i2 - 1] * A[i1][N - i2 - 3]) / B[i1][N - 3 - i2];
+				}
+			}
+			for (i1 = 1; i1 < N; i1++) {
+				#pragma species kernel X[i1-1:i1,0:N-1]|chunk(i1-1:i1,0:0) ^ A[i1:i1,0:N-1]|element ^ B[i1-1:i1,0:N-1]|chunk(i1-1:i1,0:0) -> X[i1:i1,0:N-1]|element ^ B[i1:i1,0:N-1]|element
+				for (i2 = 0; i2 < N; i2++) {
+					X[i1][i2] = X[i1][i2] - X[i1 - 1][i2] * A[i1][i2] / B[i1 - 1][i2];
+					B[i1][i2] = B[i1][i2] - A[i1][i2] * A[i1][i2] / B[i1 - 1][i2];
+				}
+				#pragma species endkernel adi_k14
+			}
+			#pragma species kernel X[N-1:N-1,0:N-1]|element ^ B[N-1:N-1,0:N-1]|element -> X[N-1:N-1,0:N-1]|element
+			for (i2 = 0; i2 < N; i2++) {
+				X[N - 1][i2] = X[N - 1][i2] / B[N - 1][i2];
+			}
+			#pragma species endkernel adi_k6
+			for (i1 = 0; i1 < N - 2; i1++) {
+				#pragma species kernel X[N-i1-3:N-i1-2,0:N-1]|chunk(N-i1-3:N-i1-2,0:0) ^ A[N-i1-3:N-i1-3,0:N-1]|element ^ B[N-i1-2:N-i1-2,0:N-1]|element -> X[N-i1-2:N-i1-2,0:N-1]|element
+				for (i2 = 0; i2 < N; i2++) {
+					X[N - 2 - i1][i2] = (X[N - 2 - i1][i2] - X[N - i1 - 3][i2] * A[N - 3 - i1][i2]) / B[N - 2 - i1][i2];
+				}
+				#pragma species endkernel adi_k18
+			}
+		}
+	}
+	#pragma endscop
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/test/examples/benchmarks/PolyBench/atax_species.c

@@ -0,0 +1,69 @@
+//
+// This file is part of the Bones source-to-source compiler examples. The C-code
+// is largely identical in terms of functionality and variable naming to the code
+// found in PolyBench/C version 3.2. For more information on PolyBench/C or Bones
+// please use the contact information below.
+//
+// == More information on PolyBench/C
+// Contact............Louis-Noel Pouchet <pouchet@cse.ohio-state.edu>
+// Web address........http://polybench.sourceforge.net/
+// 
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........benchmark/atax.c
+// Author.............Cedric Nugteren
+// Last modified on...03-April-2012
+//
+
+#include "common.h"
+
+// This is 'atax', a matrix transpose and vector multiplication kernel
+int main(void) {
+	int i,j;
+	
+	// Declare arrays on the stack
+	float A[NX][NY];
+	float x[NY];
+	float y[NX];
+	float tmp[NX];
+	
+	// Set the input data
+	for (i=0; i<NY; i++) {
+		x[i] = i*3.14159;
+	}
+	for (i=0; i<NX; i++) {
+		for (j=0; j<NY; j++) {
+			A[i][j] = ((float) i*(j+1)) / NX;
+		}
+	}
+	
+	// Perform the computation (y := A'Ax)
+	#pragma scop
+	{
+		#pragma species kernel A[0:NX-1,0:NY-1]|chunk(0:0,0:NY-1) ^ x[0:NY-1]|full -> tmp[0:NX-1]|element
+		for (i = 0; i < NX; i++) {
+			tmp[i] = 0;
+			for (j = 0; j < NY; j++) {
+				tmp[i] = tmp[i] + A[i][j] * x[j];
+			}
+		}
+		#pragma species endkernel atax_k1
+		#pragma species kernel A[0:NX-1,0:NY-1]|chunk(0:NX-1,0:0) ^ tmp[0:NX-1]|full -> y[0:NY-1]|element
+		for (j = 0; j < NY; j++) {
+			y[j] = 0;
+			for (i = 0; i < NX; i++) {
+				y[j] = y[j] + A[i][j] * tmp[i];
+			}
+		}
+		#pragma species endkernel atax_k2
+	}
+	#pragma endscop
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/test/examples/benchmarks/PolyBench/bicg_species.c

@@ -0,0 +1,71 @@
+//
+// This file is part of the Bones source-to-source compiler examples. The C-code
+// is largely identical in terms of functionality and variable naming to the code
+// found in PolyBench/C version 3.2. For more information on PolyBench/C or Bones
+// please use the contact information below.
+//
+// == More information on PolyBench/C
+// Contact............Louis-Noel Pouchet <pouchet@cse.ohio-state.edu>
+// Web address........http://polybench.sourceforge.net/
+// 
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........benchmark/bicg.c
+// Author.............Cedric Nugteren
+// Last modified on...03-April-2012
+//
+
+#include "common.h"
+
+// This is 'bicg', a biconjugate gradients method kernel
+int main(void) {
+	int i,j;
+	
+	// Declare arrays on the stack
+	float A[NX][NY];
+	float p[NY];
+	float q[NX];
+	float r[NX];
+	float s[NY];
+	
+	// Set the input data
+	for (i=0; i<NY; i++) {
+		p[i] = i*3.14159;
+	}
+	for (i=0; i<NX; i++) {
+		r[i] = i*3.14159;
+		for (j=0; j<NY; j++) {
+			A[i][j] = ((float) i*(j+1)) / NX;
+		}
+	}
+	
+	// Perform the computation
+	#pragma scop
+	{
+		#pragma species kernel r[0:NX-1]|full ^ A[0:NX-1,0:NY-1]|chunk(0:NX-1,0:0) -> s[0:NY-1]|element
+		for (j = 0; j < NY; j++) {
+			s[j] = 0;
+			for (i = 0; i < NX; i++) {
+				s[j] = s[j] + r[i] * A[i][j];
+			}
+		}
+		#pragma species endkernel bicg_k1
+		#pragma species kernel A[0:NX-1,0:NY-1]|chunk(0:0,0:NY-1) ^ p[0:NY-1]|full -> q[0:NX-1]|element
+		for (i = 0; i < NX; i++) {
+			q[i] = 0;
+			for (j = 0; j < NY; j++) {
+				q[i] = q[i] + A[i][j] * p[j];
+			}
+		}
+		#pragma species endkernel bicg_k2
+	}
+	#pragma endscop
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/test/examples/benchmarks/PolyBench/cholesky_species.c

@@ -0,0 +1,68 @@
+//
+// This file is part of the Bones source-to-source compiler examples. The C-code
+// is largely identical in terms of functionality and variable naming to the code
+// found in PolyBench/C version 3.2. For more information on PolyBench/C or Bones
+// please use the contact information below.
+//
+// == More information on PolyBench/C
+// Contact............Louis-Noel Pouchet <pouchet@cse.ohio-state.edu>
+// Web address........http://polybench.sourceforge.net/
+// 
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........benchmark/cholesky.c
+// Author.............Cedric Nugteren
+// Last modified on...03-Jul-2012
+//
+
+#include "common.h"
+
+// This is 'cholesky', a cholesky decomposition kernel
+int main(void) {
+	int i,j,k;
+	float x[1];
+	float p_i;
+	
+	// Declare arrays on the stack
+	float A[N][N];
+	float p[N];
+	
+	// Set the input data
+	for (i=0; i<N; i++) {
+		for (j=0; j<N; j++) {
+			A[i][j] = i*2.3 + 1;
+		}
+	}
+	
+	// Perform the computation
+	#pragma scop
+	{
+		for (i = 0; i < N; i++) {
+			x[0] = A[i][i];
+			#pragma species kernel x[0:0]|full ^ A[i:i,0:i-1]|element -> x[0:0]|shared
+			for (j = 0; j <= i - 1; j++) {
+				x[0] = x[0] - A[i][j] * A[i][j];
+			}
+			#pragma species endkernel cholesky_k3
+			p[i] = 1.0 / sqrt(x[0]);
+			p_i = p[i];
+			#pragma species kernel A[i:i,N-1:i+1]|element ^ A[N-1:i+1,0:i-1]|chunk(0:0,0:i-1) ^ A[i:i,0:i-1]|full -> x[0:0]|shared ^ A[N-1:i+1,i:i]|element
+			for (j = i + 1; j < N; j++) {
+				x[0] = A[i][j];
+				for (k = 0; k <= i - 1; k++) {
+					x[0] = x[0] - A[j][k] * A[i][k];
+				}
+				A[j][i] = x[0] * p_i;
+			}
+			#pragma species endkernel cholesky_k6
+		}
+	}
+	#pragma endscop
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}

data/test/examples/benchmarks/PolyBench/correlation_species.c

@@ -0,0 +1,97 @@
+//
+// This file is part of the Bones source-to-source compiler examples. The C-code
+// is largely identical in terms of functionality and variable naming to the code
+// found in PolyBench/C version 3.2. For more information on PolyBench/C or Bones
+// please use the contact information below.
+//
+// == More information on PolyBench/C
+// Contact............Louis-Noel Pouchet <pouchet@cse.ohio-state.edu>
+// Web address........http://polybench.sourceforge.net/
+// 
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........benchmark/correlation.c
+// Author.............Cedric Nugteren
+// Last modified on...07-Feb-2013
+//
+
+#include "common.h"
+
+// This is 'correlation', a correlation computation algorithm
+int main(void) {
+	int i,j,j1,j2;
+	float meanj;
+	
+	// Declare arrays on the stack
+	float data[N][M];
+	float mean[M];
+	float stddev[M];
+	float symmat[M][M];
+	
+	// Set the constants
+	float float_n = 1.2;
+	float eps = 0.1;
+	
+	// Set the input data
+	for (i=0; i<N; i++) {
+		for (j=0; j<M; j++) {
+			data[i][j] = ((float) i*j) / M;
+		}
+	}
+	
+	// Perform the computation
+	#pragma scop
+	{
+		#pragma species kernel data[0:N-1,0:M-1]|chunk(0:N-1,0:0) -> mean[0:M-1]|element
+		for (j = 0; j < M; j++) {
+			mean[j] = 0.0;
+			for (i = 0; i < N; i++) {
+				mean[j] += data[i][j];
+			}
+			mean[j] /= float_n;
+		}
+		#pragma species endkernel correlation_k1
+		#pragma species kernel mean[0:M-1]|element ^ data[0:N-1,0:M-1]|chunk(0:N-1,0:0) -> stddev[0:M-1]|element
+		for (j = 0; j < M; j++) {
+			stddev[j] = 0.0;
+			meanj = mean[j];
+			for (i = 0; i < N; i++) {
+				stddev[j] += (data[i][j] - meanj) * (data[i][j] - meanj);
+			}
+			stddev[j] /= float_n;
+			stddev[j] = sqrt(stddev[j]);
+			stddev[j] = stddev[j] <= eps ? 1.0 : stddev[j];
+		}
+		#pragma species endkernel correlation_k2
+		#pragma species kernel data[0:N-1,0:M-1]|element ^ mean[0:M-1]|element ^ stddev[0:M-1]|element -> data[0:N-1,0:M-1]|element
+		for (i = 0; i < N; i++) {
+			for (j = 0; j < M; j++) {
+				data[i][j] -= mean[j];
+				data[i][j] /= sqrt(float_n) * stddev[j];
+			}
+		}
+		#pragma species endkernel correlation_k3
+	}
+	#pragma endscop
+	
+	// Calculate the MxM correlation matrix
+	for (j1=0; j1<M-1; j1++) {
+		symmat[j1][j1] = 1.0;
+		for (j2=j1+1; j2<M; j2++) {
+			symmat[j1][j2] = 0.0;
+			for (i = 0; i<N; i++) {
+				symmat[j1][j2] += (data[i][j1] * data[i][j2]);
+			}
+			symmat[j2][j1] = symmat[j1][j2];
+		}
+	}
+	symmat[M-1][M-1] = 1.0;
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/test/examples/benchmarks/PolyBench/covariance_species.c

@@ -0,0 +1,78 @@
+//
+// This file is part of the Bones source-to-source compiler examples. The C-code
+// is largely identical in terms of functionality and variable naming to the code
+// found in PolyBench/C version 3.2. For more information on PolyBench/C or Bones
+// please use the contact information below.
+//
+// == More information on PolyBench/C
+// Contact............Louis-Noel Pouchet <pouchet@cse.ohio-state.edu>
+// Web address........http://polybench.sourceforge.net/
+// 
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........benchmark/covariance.c
+// Author.............Cedric Nugteren
+// Last modified on...26-Jun-2012
+//
+
+#include "common.h"
+
+// This is 'covariance', a covariance computation algorithm
+int main(void) {
+	int i,j,j1,j2;
+	
+	// Declare arrays on the stack
+	float data[N][M];
+	float mean[M];
+	float symmat[M][M];
+	
+	// Set the constants
+	float float_n = 1.2;
+	
+	// Set the input data
+	for (i=0; i<N; i++) {
+		for (j=0; j<M; j++) {
+			data[i][j] = ((float) i*j) / M;
+		}
+	}
+	
+	// Perform the computation
+	#pragma scop
+	{
+		#pragma species kernel data[0:N-1,0:M-1]|chunk(0:N-1,0:0) -> mean[0:M-1]|element
+		for (j = 0; j < M; j++) {
+			mean[j] = 0.0;
+			for (i = 0; i < N; i++) {
+				mean[j] += data[i][j];
+			}
+			mean[j] /= float_n;
+		}
+		#pragma species endkernel covariance_k1
+		#pragma species kernel data[0:N-1,0:M-1]|element ^ mean[0:M-1]|element -> data[0:N-1,0:M-1]|element
+		for (i = 0; i < N; i++) {
+			for (j = 0; j < M; j++) {
+				data[i][j] -= mean[j];
+			}
+		}
+		#pragma species endkernel covariance_k2
+	}
+	#pragma endscop
+	
+	// Calculate the MxM covariance matrix
+	for (j1=0; j1<M; j1++) {
+		for (j2=j1; j2<M; j2++) {
+			symmat[j1][j2] = 0.0;
+			for (i=0; i<N; i++) {
+				symmat[j1][j2] += data[i][j1] * data[i][j2];
+			}
+			symmat[j2][j1] = symmat[j1][j2];
+		}
+	}
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}