Python C拡張モジュールのSIGSEGV

Question

_xdr_read_xtc関数を使ってxtc軌道を読み込むPythonのC拡張モジュールがあります。

モジュールは問題なく.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); 

コードに間違いがありますか？無効なポインタかもしれない？それとも記憶上の問題でしょうか？

ありがとうございます！

Answer 1

xdにメモリを割り当てて、XDRFILEを保持します。次にxdをdummyに移動し、整数（ファイルハンドル？）を解析し、xd/dummyに入れます。その後、無料でdummyを解放してxdを解放します。その後、フリーメモリにアクセスするread_xtc(xd, ...に電話します。あなたのメモリアロケータがそのページを放棄することを決めた場合、あなたはSIGSEGVを取得します。

free(dummy)を実際にメモリが不要になる場所に移動します。