我有一个用于Python的C扩展模块,它使用函数_xdr_read_xtc读取xtc轨迹。
该模块内置在 .so 库中没有问题,并且大部分时间运行良好。但是,有时我会收到"分段错误(核心转储)"。
static PyObject * _xdr_read_xtc(PyObject *self, PyObject *args)
{
int natoms;
XDRFILE *xd;
xd = (XDRFILE *) malloc(sizeof(XDRFILE));
if (xd == NULL){
return NULL;}
XDRFILE *dummy;
dummy = xd;
if (!PyArg_ParseTuple(args, "ii", &xd, &natoms)){
return NULL;
}
free(dummy);
int step = 0;
float time;
float prec;
matrix box;
rvec *x;
x = malloc(natoms * sizeof(*x));
if (x == NULL){
return NULL;}
// read frame
int status = read_xtc(xd, natoms, &step, &time, box, x, &prec);
if (status == 0 | status == 11){
npy_intp dims[2]= {natoms, 3};
PyArrayObject *matout = (PyArrayObject *) PyArray_SimpleNewFromData(2, dims, NPY_FLOAT, x);
PyArray_ENABLEFLAGS(matout, NPY_ARRAY_OWNDATA);
PyObject *Frame = Py_BuildValue("Oii", matout, status, step);
Py_DECREF(matout);
return Frame;
}
else{
free(x);
return NULL;
}
}
使用 Valgrind 进行调试时,我得到"进程终止,默认操作为信号 11 (SIGSEGV)。访问不在地址 0x195688988' 的映射区域内:
int status = read_xtc(xd, natoms, &step, &time, box, x, &prec);
代码有什么明显的问题吗?也许是无效指针?或者可能是内存问题?
谢谢!
您在xd
中分配内存以保存XDRFILE
。然后你xd
移动到dummy
,解析一个整数(文件句柄?)并放入xd
/dummy
。然后你通过释放dummy
来释放xd
。然后调用将访问释放内存的read_xtc(xd, ...
。如果您的内存分配器决定放弃该页面,您将获得一个SIGSEGV。
将free(dummy)
移动到实际上不再需要内存的位置。