矩阵相乘SSE优化算法

#include   
#include
#include
#include
#include
#define maxN 10
using namespace std;
void swap（float &afloat &b）//参数传递问题
{
	float t;
	t=a;
	a=b;
	b=t;
}
//矩阵相乘平凡算法
void matrixmulti（int nfloat **afloat **bfloat **c）
{
	int ijk;
	for（i=0;i	{
		for（j=0;j		{
			c[i][j]=0;
			for（k=0;k			{
				c[i][j]=a[i][k]*b[k][j];
			}
		}
	}
}
//矩阵相乘SSE优化算法
void sse_mul（int n float **a float **b float **c）{
	__m128 t1 t2 sum;
            //turning matrix into T（b）
	for （int i = 0; i < n; ++i） for （int j = 0; j < i; ++j） swap（b[i][j] b[j][i]）;
	for （int i = 0; i < n; ++i）{
		for （int j = 0; j < n; ++j）{
			c[i][j] = 0.0;
			sum = _mm_setzero_ps（）;
			for （int k = n - 4; k >= 0; k -= 4）{