Matrix.h

#ifndef Matrix_h
#define Matrix_h
#ifndef MAX_SIZE   ///<矩阵最大维度
#define MAX_SIZE 50  
#endif

#ifndef DBL_EPSILON   ///<矩阵求逆0
#define DBL_EPSILON 0.001  
#endif

#include<stdio.h>
#include<string.h>
#include<math.h>
 
/**************************
 *@brief  Matrix  矩阵
 **************************/
struct Matrix
{
    int Row;  ///< 行
    int Col;  ///< 列
    double Val[MAX_SIZE*MAX_SIZE];
};

/**********************************************
 *@brief  Matrix_Add   矩阵加法 m3=m1+m2
 
 *@param  m1   I   输入矩阵m1
 *@param  m2   I   输入矩阵m2
 *@param  m3   O   输出矩阵m3
 **********************************************/
void Matrix_Add(const struct Matrix* m1,const struct Matrix *m2,struct Matrix *m3);

/**********************************************
 *@brief  Matrix_Subtract   矩阵减法 m3=m1-m2
 
 *@param  m1   I   输入矩阵m1
 *@param  m2   I   输入矩阵m2
 *@param  m3   O   输出矩阵m3
 **********************************************/
void Matrix_Subtract(const struct Matrix* m1,const struct Matrix *m2,struct Matrix *m3);  

/**********************************************
 *@brief  Matrix_Mutiply   矩阵乘法 m3=m1*m2
 
 *@param  m1   I   输入矩阵m1
 *@param  m2   I   输入矩阵m2
 *@param  m3   O   输出矩阵m3
 **********************************************/
void Matrix_Mutiply(const struct Matrix* m1,const struct Matrix *m2, struct Matrix *m3);

/**********************************************
 *@brief  Matrix_Transpose   矩阵转置 m2=m1_T
 
 *@param  m1   I   输入矩阵m1
 *@param  m2   O   输出矩阵m2
 **********************************************/
void Matrix_Transpose(const struct Matrix* m1, struct Matrix* m2);

/**********************************************
 *@brief  Matrix_Transpose   矩阵求逆 m2=inv(m1)
 
 *@param  m1   I   输入矩阵m1
 *@param  m2   O   输出矩阵m2
 **********************************************/
int Matrix_Inv(const struct Matrix* m1,struct Matrix* m2);


#endif