操作系统概念LAB4—矩阵相乘—实验报告.docxVIP

LAB4 实验报告 实验目的：1.矩阵相乘 实验内容：给定两个矩阵 A 和 B，其中 A 是具有 M 行、K 列的矩阵，B 为 K 行、 N 列矩阵，A 和 B 的矩阵积为 C，C 为 M 行，N 列。矩阵 C 中第 i 行、第 j 列的元 素 Ci,j 就是矩阵 A 第 i 行每个元素和矩阵 B 第 j 列每个元素乘积的和。本实验的目 的是通过对矩阵乘法是多线程实现，让学生加深对线程的理解与调用。 实验步骤： 在 Linux 操作系统上，利用 Linux 处理多线程的优点，编写程序用多线程实现 矩阵乘法。实现 A、B 两个矩阵的乘法，并输出计算结果。 实验代码如下： #include stdio.h #include time.h #include pthread.h #include stdlib.h #include unistd.h #include memory.h #define RANGE 100 #define M 2 #define K 2 #define N 2 void gene_matrix(); void read_matrix(); int matrixA[M][N]={{1,1},{2,2}}; int matrixB[N][M]={{4,14},{5,25}}; int arr[M][M][N]; int res[M][N]={0}; void *func(void *arg); main() { clock_t start=clock(); int i,j,k,cnt,cnt2; pthread_t tids[N]; for(i=0;iN;i++) { if(pthread_create(tids[i],NULL,func,(void *)i)) { perror(pthread_create); exit(1); } pthread_join(tids[i],NULL); } for(i=0;iN;i++) pthread_join(tids[i],NULL); for(i=0;iM;i++) for(j=0;jM;j++) { for(k=0;kN;k++) res[i][j]+=arr[i][j][k]; printf(res[%d][%d]=%d\n,i,j,res[i][j]); } clock_t finish=clock(); printf(n,(long)(finish-start)/1E6); } void *func(void *arg) { int k=*(int *)arg; int i,j; printf(Enter func\n); for(i=0;iM;i++) { for(j=0;jN;j++) { arr[i][j][k]=matrixA[i][k]*matrixB[k][j]; printf(func=%d\n, arr[i][j][k]); } } pthread_exit(NULL); } void gene_matrix() { FILE *file1,*file2; if((file1=fopen(/home/timelong/Desktop/matrixA,wt))==NULL) { perror(fopen); exit(1); } if((file2=fopen(/home/timelong/Desktop/matrixB,wt))==NULL) { perror(fopen); exit(1); } int i,j; srand((unsigned)time(NULL)); for(i=0;iM;i++) { for(j=0;jN;j++) fprintf(file1,%-8d,rand()%RANGE); fprintf(file1,\n); } fclose(file1); for(i=0;iN;i++) { for(j=0;jM;j++) fprintf(file2,%-8d,rand()%RANGE); fprintf(file2,\n); } fclose(file2); } void read_matrix() { FILE *file1,*file2; if((file1=fopen(/home/timelong/Desktop/matrixA,rt))==NULL) { perror(fopen); exit(1); } if((file2=fopen(/home/timelong/Desktop/matrixB,rt))==NULL) { perror(fopen); exit(1); } int i,j; for(i=0;iM;i++) for(j=0;jN;j++) { fscanf(file1,%d,matrixA[