#include <stdio.h>
#include <omp.h>
#define NRA 1000
#define NCA 1000
#define NCB 7

void main()
{
	int tid,nthreads,i,j,k,chuck;
	double a[NRA][NCA];
	double b[NCA][NCB];
	double c[NRA][NCB];
	chuck=10;
	
	#pragma omp parallel shared(a,b,c,nthreads,chuck) private (tid, i,j,k)
	{
		tid=omp_get_thread_num();
		if (tid==0)
		{	
			nthreads=omp_get_num_threads();
			printf("Matris çarpımı için %d thread kullanılmaktadır\n",nthreads);
		}

	#pragma omp for schedule(static, chuck)
	for (i=0;i<NRA;i++)
		for(j=0;j<NCA;j++)
			a[i][j]=i+j;

	#pragma omp for schedule(static, chuck)
	for (i=0;i<NCA;i++)
		for(j=0;j<NCB;j++)
			b[i][j]=i*j;

	#pragma omp for schedule(static, chuck)
	for (i=0;i<NRA;i++)
		for(j=0;j<NCB;j++)
			c[i][j]=0;
		

	#pragma omp for schedule(static, chuck)
	for (i=0;i<NRA;i++)
		{
			for(j=0;j<NCB;j++)
				for (k=0;k<NCA;k++)
					c[i][j]+=a[i][k]*b[k][j];
		}
	}
	for (i=0;i<NRA;i++)
		{
			for(j=0;j<NCB;j++)
				printf("%f   ", c[i][j]);
			printf("\n");

		}
			
	
}