我正在使用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);
}
答案 0 :(得分:1)
你的支架不匹配。您可以在if(threadno==0)
if部分后立即关闭并行部分。
此外,您似乎在代码末尾有一个额外的大括号。
正如戴夫所说,你的陈述中也缺少一些支撑。这不是问题的原因,但会导致算法无法返回预期的结果。
在C
for(i = 0; i < n; i++)
function(i);
用于在for循环中执行单个操作
如果希望for循环包含多个操作,请使用大括号
for(i = 0; i < n; i++){
function(i);
function2(i);
function3(i);
}