nt=omp_get_num_threads()
#pragma omp parallel private(i,id) reduction(+:sum) shared(nt)
{
	id = omp_get_thread_num();
	for(i=id;i<SIZE;i+=nt)
	{
		sum += calculate(i);
	}
}