bones-compiler 1.1.0

data/examples/element/example4.c

@@ -0,0 +1,49 @@
+//
+// 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...........element/example4.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+
+// This is 'example4', demonstrating two input arrays and an inner-loop as the computational body
+int main(void) {
+	int i,l;
+	float factor = 0;
+	
+	// Declare input/output arrays
+	float A[700];
+	float B[700];
+	float C[700];
+	
+	// Set the input data
+	for(i=0;i<700;i++) {
+		A[i] = i*2.3;
+		B[i] = i+6.0;
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:699|element ^ 0:699|element -> 0:699|element
+	for(i=0;i<700;i++) {
+		factor = 0.5;
+		for (l=0;l<3;l++) {
+			factor = factor + 0.2*factor;
+		}
+		C[i] = factor*A[i] + factor*B[i];
+	}
+	#pragma species endkernel example4
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/element/example5.c

@@ -0,0 +1,56 @@
+//
+// 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...........element/example5.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#define SIZE 2311
+
+// This is 'example5', demonstrating multiple inputs and outputs of different types
+int main(void) {
+	int i;
+	float result = 0;
+	
+	// Declare input/output arrays
+	float in1[SIZE];
+	int in2[SIZE];
+	char in3[SIZE];
+	float out1[SIZE];
+	float out2[SIZE];
+	
+	// Set the input data
+	for(i=0;i<SIZE;i++) {
+		in1[i] = -0.34*i;
+		in2[i] = i*3;
+		in3[i] = i%100;
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:SIZE-1|element ^ 0:SIZE-1|element ^ 0:SIZE-1|element -> 0:SIZE-1|element ^ 0:SIZE-1|element
+	for(i=0;i<SIZE;i++) {
+		if (in3[i] > 50) {
+			result = in2[i] / in1[i];
+		}
+		else {
+			result = in2[i] * in1[i];
+		}
+		out1[i] = result;
+		out2[i] = in3[i]/255.0;
+	}
+	#pragma species endkernel example5
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/element/example6.c

@@ -0,0 +1,46 @@
+//
+// 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...........element/example6.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#include <stdlib.h>
+
+// This is 'example6', demonstrating dynamically sized arrays, a dynamically sized kernel, and a classification including variables
+int main(void) {
+	int i;
+	int N = 2048*2048;
+	
+	// Declare input/output arrays
+	int *A = (int *)malloc(N*sizeof(int));
+	int *B = (int *)malloc(N*sizeof(int));
+	
+	// Set the input data
+	for(i=0;i<N;i++) {
+		A[i] = i;
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:N-1|element -> 0:N-1|element
+	for(i=0;i<N;i++) {
+		B[i] = A[i] + 3;
+	}
+	#pragma species endkernel example6
+	
+	// Clean-up and exit the function
+	free(A);
+	free(B);
+	fflush(stdout);
+	return 0;
+}
+

data/examples/element/example7.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...........element/example7.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+
+#include <stdio.h>
+#include <stdlib.h>
+
+// This is 'example7', demonstrating variable length arrays (C99 VLAs)
+int main(void) {
+	int a,b,c;
+	int dim_1A = 50;
+	int dim_1B = 40;
+	
+	// Declare input/output arrays
+	int in1[dim_1A][dim_1B][30];
+	int in2[dim_1A][dim_1B][30];
+	int out[dim_1A][dim_1B][30];
+	
+	// Set the input data
+	for(a=0;a<dim_1A;a++) {
+		for(b=0;b<dim_1B;b++) {
+			for(c=0;c<30;c++) {
+				in1[a][b][c] = a*b+c;
+				in2[a][b][c] = -(c/4)+6;
+			}
+		}
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:dim_1A-1,0:dim_1B-1,0:29|element ^ 0:dim_1A-1,0:dim_1B-1,0:29|element -> 0:dim_1A-1,0:dim_1B-1,0:29|element
+	for(a=0;a<dim_1A;a++) {
+		for(b=0;b<dim_1B;b++) {
+			for(c=0;c<30;c++) {
+				out[a][b][c] = in1[a][b][c] + in2[a][b][c];
+			}
+		}
+	}
+	#pragma species endkernel example7
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/element/example8.c

@@ -0,0 +1,45 @@
+//
+// 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...........element/example8.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+
+// This is 'example8', demonstrating a reading and writing from the same array
+int main(void) {
+	int i,j;
+	
+	// Declare input/output arrays
+	int A[100][16];
+	
+	// Set the input data
+	for(i=0;i<100;i++) {
+		for(j=0;j<16;j++) {
+			A[i][j] = i+j;
+		}
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:99,0:15|element -> 0:99,0:15|element
+	for(i=0;i<100;i++) {
+		for(j=0;j<16;j++) {
+			A[i][j] = 2*A[i][j];
+		}
+	}
+	#pragma species endkernel example8
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/element/example9.c

@@ -0,0 +1,48 @@
+//
+// 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...........element/example9.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#include <stdlib.h>
+
+// This is 'example9', demonstrating a for-loop that does not start at zero
+int main(void) {
+	int i;
+	int result;
+	int N = 2048*2048;
+	
+	// Declare input/output arrays
+	int *A = (int *)malloc(N*sizeof(int));
+	int *B = (int *)malloc(N*sizeof(int));
+	
+	// Set the input data
+	for(i=0;i<N;i++) {
+		A[i] = i;
+	}
+	
+	// Perform the computation
+	#pragma species kernel 2:N-1|element -> 2:N-1|element
+	for(i=2;i<N;i++) {
+		result = A[i] + 3;
+		B[i] = result * 3;
+	}
+	#pragma species endkernel example9
+	
+	// Clean-up and exit the function
+	free(A);
+	free(B);
+	fflush(stdout);
+	return 0;
+}
+

data/examples/neighbourhood/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...........neighbourhood/example1.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#define SIZE 60000
+#define NB 2
+
+// This is 'example1', demonstrating a basic 1D neighbourhood-based computation whose size is set by a define
+int main(void) {
+	int i,n;
+	float result = 0;
+	
+	// Declare input/output arrays
+	float A[SIZE];
+	float B[SIZE];
+	
+	// Set the input data
+	for(i=0;i<SIZE;i++) {
+		A[i] = i/2.0;
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:SIZE-1|neighbourhood(-NB:NB) -> 0:SIZE-1|element
+	for(i=0;i<SIZE;i++) {
+		if (i >= NB && i < SIZE-NB) {
+			result = 0;
+			for (n=-NB;n<=NB;n++) {
+				result = result + A[i+n];
+			}
+			B[i] = result / (NB*2+1);
+		}
+		else {
+			B[i] = A[i];
+		}
+	}
+	#pragma species endkernel example1
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/neighbourhood/example2.c

@@ -0,0 +1,55 @@
+//
+// 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...........neighbourhood/example2.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#define A 256
+#define B 512
+
+// This is 'example2', demonstrating a 2D array, a 2D neighbourhood and a for-loop-less notation of the neighbourhood accesses
+int main(void) {
+	int i,j;
+	
+	// Declare input/output arrays
+	float in[A][B];
+	float out[A][B];
+	
+	// Set the input data
+	for(i=0;i<A;i++) {
+		for(j=0;j<B;j++) {
+			in[i][j] = i+j;
+		}
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:255,0:511|neighbourhood(-1:1,-1:1) -> 0:255,0:511|element
+	for(i=0;i<A;i++) {
+		for(j=0;j<B;j++) {
+			if (i >= 1 && j >= 1 && i < (A-1) && j < (B-1)) {
+				out[i][j] = (in[i+1][j+1] + in[i+1][ j ] + in[i+1][j-1] +
+				             in[ i ][j+1] + in[ i ][ j ] + in[ i ][j-1] +
+				             in[i-1][j+1] + in[i-1][ j ] + in[i-1][j-1])/9.0;
+			}
+			else {
+				out[i][j] = in[i][j];
+			}
+		}
+	}
+	#pragma species endkernel example2
+	
+	// Clean-up and exit the function
+	fflush(stdout);
+	return 0;
+}
+

data/examples/neighbourhood/example3.c

@@ -0,0 +1,82 @@
+//
+// 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...........neighbourhood/example3.c
+// Author.............Cedric Nugteren
+// Last modified on...16-April-2012
+//
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#define A 1024
+#define B 1536
+
+// Forward declarations of helper functions for statically allocated 2D memory
+float ** alloc_2D(int size1, int size2);
+void free_2D(float ** array_2D);
+
+// This is 'example3', demonstrating a neighbourhood with only some values used (a cross) and a math.h square root function call
+int main(void) {
+	int i,j;
+	int sizea = A;
+	int sizeb = B;
+	
+	// Declare input/output arrays
+	float **in = alloc_2D(sizea,sizeb);
+	float **out = alloc_2D(sizea,sizeb);
+	
+	// Set the input data
+	for(i=0;i<sizea;i++) {
+		for(j=0;j<sizeb;j++) {
+			in[i][j] = i+j;
+		}
+	}
+	
+	// Perform the computation
+	#pragma species kernel 0:sizea-1,0:sizeb-1|neighbourhood(-1:1,-1:1) -> 0:sizea-1,0:sizeb-1|element
+	for(i=0;i<sizea;i++) {
+		for(j=0;j<sizeb;j++) {
+			if (i >= 1 && j >= 1 && i < (sizea-1) && j < (sizeb-1)) {
+				out[i][j] =                 in[i+1][ j ] + 
+				             in[ i ][j+1] + in[ i ][ j ] + in[ i ][j-1] +
+				                            in[i-1][ j ]                 ;
+			}
+			else {
+				out[i][j] = sqrt(in[i][j]);
+			}
+		}
+	}
+	#pragma species endkernel example3
+	
+	// Clean-up and exit the function
+	free_2D(in);
+	free_2D(out);
+	fflush(stdout);
+	return 0;
+}
+
+// Helper function to allocate a 2D-array
+float ** alloc_2D(int size1, int size2) {
+	int a;
+	float ** array_2D = (float **)malloc(size1*sizeof(float*));
+	float * array_1D = (float *)malloc(size1*size2*sizeof(float));
+	for (a=0; a<size1; a++) {
+		array_2D[a] = &array_1D[a*size2];
+	}
+	return array_2D;
+}
+
+// Helper function to free a 2D-array
+void free_2D(float ** array_2D) {
+	free(array_2D[0]);
+	free(array_2D);
+}
+