// 
// lock -- atomic -- critical 
//
#include <omp.h>
#include <stdio.h> 

const int n=1000000;
 
int main()
{
	long a=0;
	int i;
	double t0, t1;
	
	omp_lock_t lock;  
	omp_init_lock(&lock);   

// lock
	a = 0;
	t0=omp_get_wtime(); // 获取墙上时间，以秒为单位 
	#pragma omp parallel for num_threads(4) private(i)
		for(i=1; i<n; i++)
		{
			omp_set_lock(&lock);
			a = a + 1;
			omp_unset_lock(&lock);
		}
	omp_destroy_lock(&lock);	
	t1=omp_get_wtime(); 	
	printf("lock: a=%d, elapsed time=%.2e\n", a, t1-t0);
	
// atomic	
	a = 0;
	t0=omp_get_wtime(); // 获取墙上时间，以秒为单位 
	#pragma omp parallel for num_threads(4) private(i)
		for(i=1; i<n; i++)
		{
			#pragma omp atomic
				a = a + 1;
		}		
	t1=omp_get_wtime(); 	
	printf("atomic: a=%d, elapsed time=%.2e\n", a, t1-t0);

// atomic	
	a = 0;
	t0=omp_get_wtime(); // 获取墙上时间，以秒为单位 
	#pragma omp parallel for num_threads(4) private(i)
		for(i=1; i<n; i++)
		{
			#pragma omp critical
				a = a + 1;
		}		
	t1=omp_get_wtime(); 	
	printf("critical: a=%d, elapsed time=%.2e\n", a, t1-t0);

	return 0;
}