我正在使用 20 个数组。并使用 openMP 将 1/4 任务按顺序分配给四个线程中的每个线程。然后将整个数组的结果存储到文件中。这里有什么问题?
在第一个数组中,我为每个元素分配 i*j 值,然后进行 20 x 20 的矩阵乘法。
#include<stdio.h>
#include<omp.h>
#include "head.h"
int sum=0;
int c[20][20];
//#include<conio.h>
int main(void) {
int A[20][20],B[20][20],C[20][20];
int i,j,e;
static sum=0;
FILE *fp;
unsigned long long a,b,c,d;
int threadno;
fp=fopen("m2.txt","w");
// clrscr();
printf("\n%d \n",h[20][20]);
#pragma omp parallel shared(a,b,c,d)
{
threadno=omp_get_thread_num();
if(threadno==0)
{
for (i=0;i<5;i++)
for (j=0;j<5;j++)
A[i][j]=i*j;
B[i][j]=i*j;
for (i=0;i<5;i++)
for (j=0;j<5;j++)
for (e=0;e<5;e++)
sum+=A[i][e]*B[e][j];
C[i][j]=sum;
//fprintf(fp,"%d \t",C[i][j]);
}
}
if(threadno==1)
{
for (i=5;i<10;i++)
for (j=5;j<10;j++)
A[i][j]=i*j;
B[i][j]=i*j;
for (i=5;i<10;i++)
for (j=5;j<10;j++) {
sum=0;
for (e=5;e<10;e++)
sum+=A[i][e]*B[e][j];
C[i][j]=sum;
//fprintf(fp,"%d \t",C[i][j]);
}
}
if(threadno==2)
{
for (i=10;i<15;i++)
for (j=10;j<15;j++)
A[i][j]=i*j;
B[i][j]=i*j;
for (i=10;i<15;i++)
for (j=10;j<15;j++) {
for (e=10;e<15;e++)
sum+=A[i][e]*B[e][j];
C[i][j]=sum;
//fprintf(fp,"%d \t",C[i][j]);
}
}
if(threadno==3)
{
for (i=15;i<20;i++)
for (j=15;j<20;j++)
A[i][j]=i*j;
B[i][j]=i*j;
for (i=15;i<20;i++)
for (j=15;j<20;j++) {
for (e=15;e<20;e++)
sum+=A[i][e]*B[e][j];
C[i][j]=sum;
//fprintf(fp,"%d \t",C[i][j]);
}
}
for (i=0;i<20;i++){
for (j=0;j<20;j++) {
fprintf(fp,"%d \t",C[i][j]);
}
}
}
fclose(fp);
}