bones-compiler 1.1.0

data/examples/benchmarks/seidel-2d.c

@@ -0,0 +1,51 @@
+//
+// 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/seidel-2d.c
+// Author.............Cedric Nugteren
+// Last modified on...05-April-2012
+//
+
+#include "common.h"
+
+// This is 'seidel-2d', a 2D Seidel stencil computation
+int main(void) {
+	int i,j,t;
+	
+	// Declare arrays on the stack
+	float A[N][N];
+	
+	// Set the input data
+	for (i=0; i<N; i++) {
+		for (j=0; j<N; j++) {
+			A[i][j] = ((float) i*i*(j+2) + 2) / N;
+		}
+	}
+	
+	// Perform the computation
+	for (t=0; t<TSTEPS-1; t++) {
+		for (i=1; i<=N-2; i++) {
+			for (j=1; j<=N-2; j++) {
+				A[i][j] = (A[i-1][j-1] + A[i-1][ j ] + A[i-1][j+1]
+				         + A[ i ][j-1] + A[ i ][ j ] + A[ i ][j+1]
+				         + A[i+1][j-1] + A[i+1][ j ] + A[i+1][j+1])/9.0;
+			}
+		}
+	}
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}

data/examples/benchmarks/symm.c

@@ -0,0 +1,74 @@
+//
+// 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/symm.c
+// Author.............Cedric Nugteren
+// Last modified on...23-May-2012
+//
+
+#include "common.h"
+
+// This is 'symm', a symmetric matrix multiplication kernel
+int main(void) {
+	int i,j,k;
+	float acc[1];
+	float bij;
+	
+	// Declare arrays on the stack
+	float A[NJ][NJ];
+	float B[NI][NJ];
+	float C[NI][NJ];
+	
+	// Set the constants
+	float alpha = 32412;
+	float beta = 2123;
+	
+	// Set the input data
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NJ; j++) {
+			C[i][j] = ((float) i*j) / NI;
+			B[i][j] = ((float) i*j) / NI;
+		}
+	}
+	for (i=0; i<NJ; i++) {
+		for (j=0; j<NJ; j++) {
+			A[i][j] = ((float) i*j) / NI;
+		}
+	}
+	
+	// Perform the computation (C := alpha*A*B + beta*C, with A symmetric)
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NJ; j++) {
+			acc[0] = 0;
+			bij = B[i][j];
+			#pragma species kernel 0:j-2,j:j|element ^ 0:j-2,i:i|element -> 0:j-2,j:j|element
+			for (k=0; k<j-1; k++) {
+				C[k][j] += alpha*A[k][i] * bij;
+			}
+			#pragma species endkernel symm-part1
+			#pragma species kernel 0:j-2,j:j|element ^ 0:j-2,i:i|element -> 0:0|shared
+			for (k=0; k<j-1; k++) {
+				acc[0] += B[k][j] * A[k][i];
+			}
+			#pragma species endkernel symm-part2
+			C[i][j] = beta*C[i][j] + alpha*A[i][i]*bij + alpha*acc[0];
+		}
+	}
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/benchmarks/syr2k.c

@@ -0,0 +1,65 @@
+//
+// 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/syr2k.c
+// Author.............Cedric Nugteren
+// Last modified on...03-April-2012
+//
+
+#include "common.h"
+
+// This is 'syr2k', an algorithm for symmetric rank-2k operations
+int main(void) {
+	int i,j,k;
+	
+	// Declare arrays on the stack
+	float A[NI][NJ];
+	float B[NI][NJ];
+	float C[NI][NI];
+	
+	// Set the constants
+	int alpha = 32412;
+	int beta = 2123;
+	
+	// Set the input data
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NJ; j++) {
+			A[i][j] = ((float) i*j) / NI;
+			B[i][j] = ((float) i*j) / NI;
+		}
+	}
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NI; j++) {
+			C[i][j] = ((float) i*j) / NI;
+		}
+	}
+	
+	// Perform the computation (C := alpha*A*B' + alpha*B*A' + beta*C)
+	#pragma species kernel 0:NI-1,0:NI-1|element ^ 0:NI-1,0:NJ-1|chunk(0:0,0:NJ-1) ^ 0:NI-1,0:NJ-1|chunk(0:0,0:NJ-1) -> 0:NI-1,0:NI-1|element
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NI; j++) {
+			C[i][j] *= beta;
+			for (k=0; k<NJ; k++) {
+				C[i][j] += alpha * A[i][k] * B[j][k];
+				C[i][j] += alpha * B[i][k] * A[j][k];
+			}
+		}
+	}
+	#pragma species endkernel syr2k
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}

data/examples/benchmarks/syrk.c

@@ -0,0 +1,62 @@
+//
+// 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/syrk.c
+// Author.............Cedric Nugteren
+// Last modified on...08-May-2012
+//
+
+#include "common.h"
+
+// This is 'syrk', an algorithm for symmetric rank-k operations
+int main(void) {
+	int i,j,k;
+	
+	// Declare arrays on the stack
+	float A[NI][NJ];
+	float C[NI][NI];
+	
+	// Set the constants
+	float alpha = 32412;
+	float beta = 2123;
+	
+	// Set the input data
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NJ; j++) {
+			A[i][j] = ((float) i*j) / NI;
+		}
+	}
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NI; j++) {
+			C[i][j] = ((float) i*j) / NI;
+		}
+	}
+	
+	// Perform the computation (C := alpha*A*A' + beta*C)
+	#pragma species kernel C[0:NI-1,0:NI-1]|element ^ A[0:NI-1,0:NJ-1]|chunk(0:0,0:NJ-1) ^ A[0:NI-1,0:NJ-1]|chunk(0:0,0:NJ-1) -> C[0:NI-1,0:NI-1]|element
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NI; j++) {
+			C[i][j] *= beta;
+			for (k=0; k<NJ; k++) {
+				C[i][j] += alpha * A[i][k] * A[j][k];
+			}
+		}
+	}
+	#pragma species endkernel syrk
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}

data/examples/benchmarks/trisolv.c

@@ -0,0 +1,57 @@
+//
+// 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/trisolv.c
+// Author.............Cedric Nugteren
+// Last modified on...03-Jul-2012
+//
+
+#include "common.h"
+
+// This is 'trisolv', a solver for linear systems with triangular matrices
+int main(void) {
+	int i,j;
+	float A_i_i;
+	
+	// Declare arrays on the stack
+	float A[NX][NX];
+	float c[NX];
+	float x[NX];
+	
+	// Set the input data
+	for (i=0; i<NX; i++) {
+		c[i] = ((float) i) / NX;
+		x[i] = ((float) i) / NX;
+		for (j=0; j<NX; j++) {
+			A[i][j] = ((float) i*j) / NX;
+		}
+	}
+	
+	// Perform the computation
+	for (i=0; i<NX; i++) {
+		x[i] = c[i];
+		A_i_i = A[i][i];
+		//#pragma species kernel i:i,0:i-1|element -> i:i|shared
+		for (j=0; j<= i-1; j++) {
+			x[i] = x[i] - A[i][j] * x[j];
+			x[i] = x[i] / A_i_i;
+		}
+		//#pragma species endkernel trisolv
+	}
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}

data/examples/benchmarks/trmm.c

@@ -0,0 +1,57 @@
+//
+// 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/trmm.c
+// Author.............Cedric Nugteren
+// Last modified on...26-Jun-2012
+//
+
+#include "common.h"
+
+// This is 'trmm', a triangular matrix multiplication kernel
+int main(void) {
+	int i,j,k;
+	
+	// Declare arrays on the stack
+	float A[NI][NI];
+	float B[NI][NI];
+	
+	// Set the constants
+	int alpha = 32412;
+	
+	// Set the input data
+	for (i=0; i<NI; i++) {
+		for (j=0; j<NI; j++) {
+			A[i][j] = ((float) i*j) / NI;
+			B[i][j] = ((float) i*j) / NI;
+		}
+	}
+	
+	// Perform the computation (B := alpha*A'*B, with A triangular)
+	for (i=1; i<NI; i++) {
+		for (j=0; j<NI; j++) {
+			#pragma species kernel i:i,0:i-1|element ^ j:j,0:i-1|element -> i:i,j:j|shared
+			for (k=0; k<i; k++) {
+				B[i][j] += alpha * A[i][k] * B[j][k];
+			}
+			#pragma species endkernel trmm
+		}
+	}
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/chunk/example1.c

@@ -0,0 +1,54 @@
+//
+// This file is part of the Bones source-to-source compiler examples. This C-code
+// example is meant to illustrate the use of Bones. For more information on Bones
+// use the contact information below.
+//
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........chunk/example1.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+
+// This is 'example1', a basic chunk to element example using a 2D tile
+int main(void) {
+	int i,j;
+	int i2,j2;
+	int result = 0;
+	
+	// Declare input/output arrays
+	int A[560][32];
+	int B[56][16];
+	
+	// Set the input data
+	for(i=0;i<560;i++) {
+		for(j=0;j<32;j++) {
+			A[i][j] = i+j;
+		}
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:559,0:31|chunk(0:9,0:1) -> 0:55,0:15|element
+	for(i=0;i<56;i++) {
+		for(j=0;j<16;j++) {
+			result = 0;
+			for (i2=0;i2<10;i2++) {
+				for (j2=0;j2<2;j2++) {
+					result = result + A[i*10+i2][j*2+j2];
+				}
+			}
+			B[i][j] = result;
+		}
+	}
+	#pragma species endkernel example1
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/chunk/example2.c

@@ -0,0 +1,44 @@
+//
+// This file is part of the Bones source-to-source compiler examples. This C-code
+// example is meant to illustrate the use of Bones. For more information on Bones
+// use the contact information below.
+//
+// == More information on Bones
+// Contact............Cedric Nugteren <c.nugteren@tue.nl>
+// Web address........http://parse.ele.tue.nl/bones/
+//
+// == File information
+// Filename...........chunk/example2.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#define SIZE 2048
+#define HALFSIZE (SIZE/2)
+
+// This is 'example2', demonstrating a chunk-example without an inner-loop, everything is unrolled manually
+int main(void) {
+	int i;
+	
+	// Declare input/output arrays
+	float A[SIZE];
+	float B[HALFSIZE];
+	
+	// Set the input data
+	for(i=0;i<SIZE;i++) {
+		A[i] = i%6+i;
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:SIZE-1|chunk(0:1) -> 0:HALFSIZE-1|element
+	for(i=0;i<HALFSIZE;i++) {
+		B[i] = A[i*2] + A[i*2+1];
+	}
+	#pragma species endkernel example2
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+