在C中使用稀疏矩阵库文件
|
我正在从事一个大型项目,其中正在设计稀疏矩阵矢量应用程序,但我仍在努力理解代码。我从为应用程序建立基础开始,但是在执行程序时遇到了分段错误。我已经在MatrixRead函数中将此问题跟踪到此循环,并在下面附上代码。执行该程序时,我尝试在一些测试消息中进行编程,该程序似乎执行了所有循环,但最后返回分段错误。当然,这仅仅是猜测。任何帮助都是极好的。谢谢!
while (ret != EOF && row <= mat->rows)
{
if (row != curr_row) // Won\'t execute for first iteration
{
/* store this row */
MatrixSetRow(mat, curr_row, len, ind, val);
/* check if the previous row is zero */
i = 1;
while(row != curr_row + i)
{
mat->lens[curr_row+i-1] = 0;
mat->inds[curr_row+i-1] = 0;
mat->vals[curr_row+i-1] = 0;
i++;
}
curr_row = row;
/* reset row pointer */
len = 0;
}
ind[len] = col;
val[len] = value;
len++;
ret = fscanf(file, \"%lf %lf %lf\", &r1, &c1, &value);
col = (int) (c1);
row = (int) (r1);
}
/* Store the final row */
if (ret == EOF || row > mat->rows)
MatrixSetRow(mat, mat->rows, len, ind, val);
这是MatrixSetRow函数的代码:
/*--------------------------------------------------------------------------
* MatrixSetRow - Set a row in a matrix. Only local rows can be set.
* Once a row has been set, it should not be set again, or else the
* memory used by the existing row will not be recovered until
* the matrix is destroyed. \"row\" is in global coordinate numbering.
*--------------------------------------------------------------------------*/
void MatrixSetRow(Matrix *mat, int row, int len, int *ind, double *val)
{
row -= 1;
mat->lens[row] = len;
mat->inds[row] = (int *) MemAlloc(mat->mem, len*sizeof(int));
mat->vals[row] = (double *) MemAlloc(mat->mem, len*sizeof(double));
if (ind != NULL)
memcpy(mat->inds[row], ind, len*sizeof(int));
if (val != NULL)
memcpy(mat->vals[row], val, len*sizeof(double));
}
我还包括随之而来的Matrix.h文件的代码,其中定义了Matrix的成员:
#include <stdio.h>
#include \"Common.h\"
#include \"Mem.h\"
#ifndef _MATRIX_H
#define _MATRIX_H
typedef struct
{
int rows;
int columns;
Mem *mem;
int *lens;
int **inds;
double **vals;
}
Matrix;
没有找到相关结果
已邀请:
0 个回复