algorithm.c
aa9b50c53d27ancarola

algorithm.c
View Options

			/*
			============================================================================
			Filename : algorithm.c
			Author : Your names go here
			SCIPER : Your SCIPER numbers

			============================================================================
			*/
			#include <math.h>

			#define INPUT(I,J) input[(I)*length+(J)]
			#define OUTPUT(I,J) output[(I)*length+(J)]


			void simulate(double input, double output, int threads, int length, int iterations)
			{
			double *temp;

			omp_set_num_threads(threads);

			for(int n=0; n < iterations; n++)
			{
			//int th_index;
			int tid;
			const int B = 4; // size of block
			const int jump = threads * B;
			const int N2 = length * (length-1);

			#pragma omp parallel shared(N2, B, jump)
			{
			tid = omp_get_thread_num();

			for(int off=length + tid * B; off < N2; off += jump)
			{
			for(int k=off; k<B + off; k++)
			{
			int i = k / length;
			int j = k % length;

			if ( (j == 0) \|\| (j == length - 1)
			\|\| (((i == length/2-1) \|\| (i== length/2)) && ((j == length/2-1) \|\| (j == length/2))) )
			continue;


			OUTPUT(i,j) = ( (INPUT(i-1,j-1) + INPUT(i-1,j) + INPUT(i-1,j+1)) +
			(INPUT(i,j-1) + INPUT(i,j) + INPUT(i,j+1)) +
			(INPUT(i+1,j-1) + INPUT(i+1,j) + INPUT(i+1,j+1)) ) / 9.0;
			}
			}
			}

			temp = input;
			input = output;
			output = temp;
			}

			}