线性代数 1：基本的数据类型和数据结构选择

线性代数中涉及到了矩阵和行列式。其中行列式的行列值是相等的但矩阵就不一定满足这个条件为了方便设计，可以将行列式看成是一种特殊的矩阵。既二者公用一个结构体类型。

/*基本数据类型*/ typedef struct data{ lbian hang; /*行数*/ lbian lie; /*列数 */ datatype **top; /*指向行列式或者矩阵的首个数据的指针*/ }data; 其中数据类型的定义是： typedef unsigned char byte; typedef char int8; typedef unsigned short int uint16; typedef short int int16; typedef long int int32; typedef unsigned long int uint32; typedef double datatype; /*矩阵元素类型*/ typedef uint32 lbian; /*矩阵长宽类型*/ typedef uint32 cnttype; /*计数类型*/

结构体中的hang，lie分别代表这个矩阵的行数和列数，当然对于行列式这两个值是相等的。
top指向一块内存，这块内存存放这矩阵的数据值。

例如定义一个结构体变量data D;
并在以后为D的成员赋值：D.hang = 4 ;D.lie = 3;
通过调用app_mem （）函数为D.top成员申请内存：

void apply_mem(data *D ) { lbian i; datatype **top; top = (datatype **)calloc(D->hang+1, sizeof (datatype *)); if(!top) printf("内存错误!/n"); for(i = 1; i <= D->hang;i ++){ top[i] = (datatype *)calloc(D->lie+1,sizeof(datatype)); if(!top[i]) printf("内存错误/n"); } D ->top = top; }

因为在矩阵中对元素位置的定位的说法是n行m列，其中m和n都是大于1的数，所以应当从top[1]开始使用故多申请一块内存：
（也可以将top--；这样top[1]就是申请的第一块内存，就可以节约一块内存地址空间了）

见下图：

使用完了以后还得调用free_mem（）函数将top所指向的内存释放：

void free_mem(data *D) { lbian i; if(D != NULL ){ for(i = 1; i <= D->hang; i++){ free((D->top)[i]); } free(D->top); D ->top = NULL; /*释放后将D的top成员置为NULL；*/ } }