/****************************************矩陣基本運算庫****************************************
功能：
1、矩陣相乘?（直接相乘、分治法）
2、求解矩陣行列式值（部分選主元LU分解法）
3、特征值（QR分解）
4、特征向量（原點平移反冪法，只針對于非對稱實矩陣）
5、矩陣的逆（列主元單位陣變換法）
執行效率：一般，可使用Cache編碼技術優化
作者：李泉衡
創建時間：2013.6.19
完成測試時間：2013.6.20
注釋整理時間：2013.10.12
郵箱：lqh0220@gmail.com
說明：版權所有，僅供參考。
*********************************************************************************************/


#include?“stdafx.h“
#include?
#include?
#include?“MatrixMath.h“
?void?MatrixMulti（MATRIX_DATA_TYPE?**pOutMATRIX_DATA_TYPE?**Xlong?X_rowlong?X_colMATRIX_DATA_TYPE?**Ylong?Y_rowlong?Y_col）??//簡單直接相乘
{
	memset（pOut[0]0X_row*Y_col*sizeof（MATRIX_DATA_TYPE））;
	if（X_col!=Y_row）
		return;
	for（long?i=0;i		for（long?k=0;k			for（long?j=0;j				pOut[i][k]?+=X[i][j]*Y[j][k];
}
void?Matrix_Multi1（MATRIX_DATA_TYPE?**AMATRIX_DATA_TYPE?**BMATRIX_DATA_TYPE?**Clong?size）???//分治法1?只能是方陣
{
	long?rowcoln=0;
	long?sub_size=size/2;
	long?c=sub_size*sub_size;
	long?c0=sub_size*sizeof（MATRIX_DATA_TYPE）;
	MATRIX_DATA_TYPE?**p_sub=NULL;
	MATRIX_DATA_TYPE?***pp_sub=NULL;
	MATRIX_DATA_TYPE?*p_tmp=NULL;
	if（size!=1）
	{
		pp_sub=（MATRIX_DATA_TYPE?***）malloc（29*sizeof（MATRIX_DATA_TYPE））;
		p_sub=（MATRIX_DATA_TYPE?**）malloc（29*sub_size*sizeof（MATRIX_DATA_TYPE））;
		p_tmp=（MATRIX_DATA_TYPE?*）malloc（c*17*sizeof（MATRIX_DATA_TYPE））;
		for（int?i=0;i<29;i++）
		{
			pp_sub[i]=p_sub+i*sub_size;
			if（i>7&&i<25）
			{
				pp_sub[i][0]=p_tmp+（i-8）*c;
				for（n=1;n				{
					pp_sub[i][n]=pp_sub[i][0]+n*sub_size;
				}
			}
		}
		for（row=0;row		{
			?pp_sub[0][row]=A[row];//A11???????????????//拆分A矩陣
			?pp_sub[1][row]=pp_sub[0][row]+sub_size;//A12???????????
			?pp_sub[2][row]=A[sub_size+row];//A21
			?pp_sub[3][row]=pp_sub[2][row]+sub_size;//A22

			?pp_sub[4][row]=B[row];//B11???????????????//拆分B矩陣
			?pp_sub[5][row]=pp_sub[4][row]+sub_size;//B12???????????
			?pp_sub[6][row]=B[sub_size+row];//B21
			?pp_sub[7][row]=pp_sub[6][row]+sub_size;//B22

			?pp_sub[25][row]=C[row];//C11???????????????//拆分C矩陣
			?pp_sub[26][row]=pp_sub[25][row]+sub_size;//C12???????????
			?pp_sub[27][row]=C[sub_size+row];//C21
			?pp_sub[28][row]=pp_sub[27][row]+sub_size;//C22

			for（col=0;col			{
				pp_sub[8][row][col]=pp_sub[5][row][col]?-?pp_sub[7][row][col];????//B12-B22
				pp_sub[9][row][col]=pp_sub[0][row][col]?+?pp_sub[1][row][col];????//A11+A12
				pp_sub[10][row][col]=pp_sub[2][row][col]+?pp_sub[3][row][col];????//A21+A22
				pp_sub[11][row][col]=pp_sub[6][row][col]-?pp_sub[4][row][col];????//B21-B11
				pp_sub[12][row][col]=pp_sub[0][row][col]+?pp_sub[3][row][col];????//A11+A22
				pp_sub[13][row][col]=pp_sub[4][row][col]+?pp_sub[7][row][col];????//B11+B22
				pp_sub[14][row][col]=pp_sub[1][row][col]-?pp_sub[3][row][col];????//A12-A22
				pp_sub[15][row][col]=pp_sub[6][row][col]+?pp_sub[7][row][col];???//B21+B22

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1282??2013-10-12?16:42??MatrixMath.h

?????文件??????19711??2013-10-12?16:39??MatrixMath.cpp

-----------?---------??----------?-----??----

????????????????20993????????????????????2