csr的实现方法

采用压缩稀疏行格式（CSR）存储核心代码段

for (i=0; i<M; i++) 
{ 
    t = 0.0; 
    for (j=row_start[i]; j<row_start[i+1]; j++) { 
        t += val[j] * x[col_idx[j]]; 
    } 
    y[i]=t; //y[i]存储积的第i行
}