matrixmultiply.cpp

测量时

// MatrixMultiply.cpp : Defines the entry point for the console application.
//
//实矩阵相乘
#include "stdafx.h"
void matrixMultiply(double a[4][5],double b[5][3],double c[4][3],int m,int n,int nk)
{
	for(long i=0;i<=m-1;i++)
		for(long j=0;j<nk-1;j++)
		{
		double u=i*nk+j;
		c[i][j]=0.0;
		for (long l=0;l<=n-1;l++)
			c[i][j] +=a[i][l]*b[l][j];
		}
	return;
}
int main(int argc, char* argv[])
{
	static double a[4][5]={{1,3,-2,0,4},{-2,-1,5,-7,2},{0,8,4,1,-5},{3,3,2,-4,1}};
	static double c[4][3];
	static double	b[5][3]={{4.0,5,-1},{2,-2,6},{7,8,1},{0,3,-5},{9,8,-6}};
	matrixMultiply(a,b,c,4,5,3);
	for (long i=0;i<=3;i++)
	{
		for (long j=0;j<=2;j++)
			printf("%16e",c[i][j]);
        printf("\n");
	}

	printf("\n");
	return 0;
}