matrix_mpyc.c

自己编写的任意矩阵相乘的DSP程序

#define m 10
#define n 30
#define k 10

void matrix_mpyc(short A[m][n],short B[n][k],int C[m][k],int max[3])
{
    int i,w,j,sum;
    short a,b;

    for(i=0;i<m;i++)
	   for(w=0;w<k;w++)
	   {  sum=0;
          for(j=0;j<n;j++)
	      {
		     a=A[i][j];
			 b=B[j][w];
	         sum+=(int)a*(int)b;
	      }
		  C[i][w]=sum;
          if(sum>max[0]) {max[0]=sum;max[1]=i+1;max[2]=w+1;}
		}

}