稀疏矩阵到三角形表结构转换算法实现

void list_table(LINKED_LIST *list,int *prow,int *pcol,int *old_newr,int *old_newc,TRIANGULAR_TABLE *table)
/////////将稀疏矩阵的双链表list格式转换成三角形表table。prow/pcol是新行/列号到旧行/列号的映射表；old_newr/old_newc是旧行/列号到新行/列号的映射表
{
int n;////矩阵阶数
int i,ii,jj,kr,kc,x,y,nn;
int si,sj,sm,sim,st,sst,ss[500];
n=(*list).n;
(*table).n=n
(*table).nza=(*list).nza;
for(i=1;i<n+1;i++)   (*table).roco[i]=i;/////////对角元素优先
      nza=n+1;//////已建立n个非零元素
for(i=1;i<n+1;i++)
{
   (*table).urp[i]=nza;
   ii=prow[i];////取原行号
   kr=(*list).rp[ii];
   sm=0;
   while(kr!=0){//扫描第i/ii行，将新列号保存在ss[]
          y=(*list).col[kr];
          y=old_newc[y];/////////将旧列号转换为新列号
          if(y>i){////////U行元素
                   ss[sm]=y;
                   sm++;//////U行非零元个数
                  }
          kr=(*list).right[kr];
         }
if(sm>0){
         for(si=0;si<sm-1;si++){////////将U行非零元列号记入rocol[]
               sim=si;
               sst=ss[si];