在C擴展段錯誤中創建numpy數組

Question

我只是想通過創建一個numpy數組開始，甚至還沒有開始編寫擴展名。 這是一個超級簡單的程序：

#include <stdio.h>
#include <iostream>
#include "Python.h"
#include "numpy/npy_common.h"
#include "numpy/ndarrayobject.h"
#include "numpy/arrayobject.h"

int main(int argc, char * argv[])
{
    int n = 2;
    int nd = 1;
    npy_intp size = {1};
    PyObject* alpha = PyArray_SimpleNew(nd, &size, NPY_DOUBLE);
    return 0;
}

該程序在PyArray_SimpleNew調用上出現段PyArray_SimpleNew ，我不明白為什么。 我正在嘗試遵循以前的一些問題（例如numpy數組C api和C數組到PyArray ）。 我究竟做錯了什么？

Answer 1

例如， PyArray_SimpleNew典型用法

int nd = 2;
npy_intp dims[] = {3,2};
PyObject *alpha = PyArray_SimpleNew(nd, dims, NPY_DOUBLE);

請注意， nd的值不得超過dims[]數組的元素數。

還：擴展必須調用 import_array()來設置C API的功能指針表。 例如在Cython：

import numpy as np
cimport numpy as np

np.import_array()  # so numpy's C API won't segfault

cdef make_array():
  cdef np.npy_intp element_count = 100
  return np.PyArray_SimpleNew(1, &element_count, np.NPY_DOUBLE)