reduction.c

并行计算的演示程序

#include <omp.h>
int main ()  
{
	int   i, n, chunk;
	float a[100], b[100], result;
	/* Some initializations */
	n = 100;
	chunk = 10;
	result = 0.0;
	for (i=0; i < n; i++)  
	{  
		a[i] = i * 1.0; 
		b[i] = i * 2.0; 
 	}
	#pragma omp parallel for default(shared) private(i)
	for (i=0; i < n; i++)    
		result = result + (a[i] * b[i]);
	printf("Final result= %f\n",result);
}