#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <omp.h>

int main(int argc, char** argv)
{
	double* T;
	int Nx1=4096+4096;
	int seed=86456;
	int Total=0;
	int i, j, tnot;

	srand(seed);

	T=(double*) malloc(Nx1*sizeof(double));
	for(i=0;i<Nx1;i++){
		T[i]=rand();	
	}

#pragma omp parallel for default(shared) private(i,j) schedule(runtime) reduction(+:Total)
	for(i=0;i<Nx1;i++) {
		for(j=0;j<Nx1;j++)
			T[i] = T[i] + sqrt(fabs((cos(T[i])+sin(T[j]))))/T[j];
		Total = Total + T[i];
	}
	free(T);
	printf("total = %d\n", Total);
	
	return 0;
}