algorithm_atomic.c
a24fbaf57e4clouis

algorithm_atomic.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;

			double temp_sum = 0;

			omp_set_num_threads(threads);

			for(int n=0; n < iterations; n++)
			{
			#pragma omp parallel for collapse(2)
			for(int i=1; i<length-1; i++)
			{
			for(int j=1; j<length-1; j++)
			{
			if(!(j == length/2 - 1)\|\|(j == length/2))
			{
			temp_sum = INPUT(i-1,j) + INPUT(i,j) + INPUT(i,j+1);

			if(!(i-1 == length/2 - 1)\|\|(i-1 == length/2))
			{
			#pragma omp atomic
			OUTPUT(i-1,j) += temp_sum;
			}
			if(!(i == length/2 - 1)\|\|(i == length/2))
			{
			#pragma omp atomic
			OUTPUT(i,j) += temp_sum;
			}
			if(!(i+1 == length/2 - 1)\|\|(i+1 == length/2))
			{
			#pragma omp atomic
			OUTPUT(i+1,j) += temp_sum;
			}
			}
			}
			}

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

			}