下列代码是如何读入一个稀疏矩阵的
struct sparseStruct
{
int numRow;
int numCol;
int numVal;
int *colPtr;
int *rowIdx;
double *value;
};
int *readMatrix(char *fname, sparseStruct *mat, int &numEmptyCol)
//从文件读入一个稀疏矩阵到mat中，返回一些空向量的IDs
{
char filename[FILENAME_LENGTH];
sprintf(filename, "%s%s", fname, "_dim");
std::ifstream dimFile(filename);//包含行数、列数、值的个数的文档
sprintf(filename, "%s%s", fname, "_row_ccs");
std::ifstream rowPtrFile(filename);//包含行ID的文档
sprintf(filename, "%s%s", fname, "_col_ccs");
std::ifstream colIdxFile(filename);//不明白是不是列ID的文档，这个文档的内容不太清楚。
sprintf(filename, "%s%s", filename, "_nz");
std::ifstream valFile(filename);//包含非零值的文档
cout << " Reading dimension file ... " << endl;
dimFile >> mat->numRow >> mat->numCol >> mat->numVal;//读入维度文件信息
dimFile.close();
mat->colPtr = new int[mat->numCol+1];
mat->rowIdx = new int[mat->numVal];
mat->value = new double[mat->numVal];
memoryUsed += (mat->numCol + 1 + mat->numVal) * sizeof(int) + mat->numVal * sizeof(double);//存储稀疏矩阵的空间值
int pre = -1;//不太清楚这一段是怎么存储空向量的
int *tempEmptyColId = new int[mat->numCol], *emptyColId;
numEmptyCol = 0;
cout << " Reading column pointer file ... " << endl;
for (int i = 0; i < mat->numCol+1; i++){
colIdxFile >> (mat->colPtr)[i];
if ((mat->colPtr)[i] == pre){
tempEmptyColId[numEmptyCol] = i - 1;
numEmptyCol++;
}
pre = (mat->colPtr)[i];
}
colIdxFile.close();
emptyColId = new int[numEmptyCol+1];
for(int i = 0; i < numEmptyCol; i++)
emptyColId[i] = tempEmptyColId[i];
emptyColId[numEmptyCol] = mat->numCol;
delete [] tempEmptyColId;
cout << " Reading row index file ... " << endl;读入行的ID
for (int i = 0; i < mat->numVal; i++)
rowPtrFile >> (mat->rowIdx)[i];
rowPtrFile.close();
cout << " Reading non-zero value file ... " << endl;读入非零的值
for (int i = 0; i < mat->numVal; i++)
valFile >> (mat->value)[i];
valFile.close();
cout << endl;
if (numEmptyCol == 0)
cout << " !!! No empty col found !!!" << endl;
else {
cout << " !!! " << numEmptyCol << " empty col(s) found !!!" << endl;
for(int i = 0; i < numEmptyCol; i++)
cout << emptyColId[i] << " ";
cout << endl;
}
return emptyColId;
}
sparseStruct *mat 是一个结构变量，存储从file中读取的稀疏矩阵，在Reading column pointer file 时遇到一些问题，比如一个稀疏矩阵
0 2 4 0 1
3 2 0 0 0
0 0 4 0 0
2 3 0 0 0
维度文件fname_dim的内容为4 5 8
行文件fname_row_ccs的内容应该是行的ID，因为从上述代码中可以看到它有8个值。
但是列文件fname_col_ccs的内容应该是什么？这个代码是怎么读入一个稀疏矩阵的？
[解决办法]
就是用：
dimFile >> mat->numRow >> mat->numCol >> mat->numVal;//读入维度文件信息
这种流的方式读的。。。
------解决方案--------------------
如果fname_row_ccs的内容是行的ID
如果fname_col_ccs的内容就应该是列的ID
