采用压缩稀疏行格式(CSR)存储核心代码段
1 for (i=0; i<M; i++) 2 { 3 t = 0.0; 4 for (j=row_start[i]; j<row_start[i+1]; j++) { 5 t += val[j] * x[col_idx[j]]; 6 } 7 y[i]=t; //y[i]存储积的第i行 8 }