作者：未知 来源：月光软件站 加入时间：2005-2-28　月光软件站

1）使用D3DXMatrixMultiply进行矩阵相乘
2）使用优化的4元数算法模拟矩阵相乘（相乘的结果与1相同）
3）最笨的办法计算矩阵相乘（要算64个乘法、48个加法、以及若干赋值）
4）使用SSE汇编计算矩阵相乘（要算64个乘法、48个加法、以及若干赋值）
5）使用xmmintrin指令集计算矩阵相乘

结果表明：
1）最快，估计不仅指令进行了优化，算法也进行了优化
2）是1的4～5倍，与3相比，证明算法起了作用
3）是1的10倍左右
4）与3差不多，我ft，这是为什么？我明明用的是汇编啊？怎么速度一点也没快！！！
5）与4差不多。

经验教训：
1）不能直接声明float[4] f这样的浮点数组，由于对其问题，在汇编语言中运行不了，应该这样声明__declspec(align(16)) float[4] f; 或者这样声明
union sse4 {
    __m128 m;
    float f[4];
};
2）那么矩阵就是这样声明的了：
union sse16 {
    __m128 m[4];
    float f[4][4];
};

相关软件：