利用稀疏矩阵的程序.docVIP

#define R 1 #define G 2 #define VS 3 #define CS 4 #define VCCS 5 #define VCVS 6 #define CCCS 7 #define CCVS 8 #define OPAMP 9 typedef struct tagBranch{ int number; //支路号 char type; //类型号 int nfrom; //始节点 int nto; //终节点 float value; //参数值 int ncfrom; //控制始节点/控制终节点；G/R支路的“设置为电压定义支路”标记 int ncto; //控制终节点；G/R支路的“控制用电压定义支路”标记 } BRANCH; BRANCH Branch0; *BRANCH; Branch = new BRANCH[BranchNumber+1]; int BranchNumber; int NodeNumber; int vBranchNumber; int NN; double val[NN]; int row[NN]; int col[NN]; int up[NN]; int down[NN]; int left[NN]; int right[NN]; int rp[NN]; int cp[NN]; int del[NN]; int nza; int n; typedef struct tagLINKED_LIST{ int nza; //非零元素个数 int n; //矩阵阶数 int *row; //元素行号 int *col; //元素列号 int *up; //元素的上邻元素号 int *down; //元素的下邻元素号 int *left; //元素的左邻元素号 int *right; //元素的右邻元素号 int *rp; //行链指针 int *cp; //列链指针 char *del; //元素删除标志；1=删除 } LINKED_LIST; typedef struct tagTRIANGULAR_TABLE{ int nza; //非零元素个数 int n; //矩阵阶数 double *val; //矩阵及右端向量/中间解向量元素值 int *roco; //对角元素行列号，U元素列号，L元素行号 int *urp; //U阵行指针 int *lcp; //L阵列指针 int lun; //由符号LU分解确定的需修正的元素个数 int fen; //由符号前消确定的需修正的向量元素个数 int *lup; //由符号LU分解确定的需修正的元素号 int *fep; //前消需修正的元素号 } TRIANGULAR_TABLE; typedef struct tagTRIANGULAR_TABLE_C{ int nza; int n; int lun; int fen; double *vre; double *vim; int *roco; int *urp; int *lcp; int *lup; int *fep; } tagTRIANGULAR_TABLE_C; #include XISHU.h //创建双重链接表 void new_list(int rank_max,int nonzero,LINKED_LIST *list) { row=new int[nonzero]; col=new int[nonzero]; up=new int[nonzero]; down=new int[nonzero]; left=new int[nonzero]; right=new int[nonzero]; rp=new int[rank_max]; cp=new int[rank_max]; del=new int [nonzero]; nza=0; //创建时尚无元素 int i; for(i=0;irank_max;i++){ //初始化 rp[i]=0; cp[i]=0; } for(i=0;inonzero+1;i++){ row[i]=0; col[i]=0; up[i]=0; down[i]=0; left[i]=0; right[i]=0; del[i]=0; } } //-------------------