Question

首先，我重申一下这个问题：使用python docs和其他资源将Python嵌入到C / C ++应用程序中时。我发现你可以使用

导入一个模块

PyObject *pName = PyUnicode_FromString((char*)"scipy");
PyObject *pModule = PyImport_Import(pName);

但如果我尝试导入＆＃34; scipy.optimize＆＃34;

PyObject *pName = PyUnicode_FromString((char*)"scipy.optimize");
PyObject *pModule = PyImport_Import(pName);

然后程序无法初始化pModule。如果我尝试再次

PyObject *pName = PyUnicode_FromString((char*)"scipy");
PyObject *pModule = PyImport_Import(pName); 
pFunc = PyObject_GetAttrString(pModule, (char*)"optimize.nnls");

我将optimize放在函数名中，无法初始化pFunc。如何在子模块中导入或调用函数，即如何调用函数scipy.optimize.nnls？

接下来，我会列出我的代码，以防它有用：

/* Relelvant imports
    #include <Python.h>
    #include <numpy/arrayobject.h> 
*/

void nnls::update(const vec& x, const vec& y)
{
    mat B;
    vecToGslVec(x, gslx);
    generateX(X, gslx);
    A = gslMatToMat(X);
    B = A.transpose();
    long double *c_out;

    Py_Initialize();
    PyObject *pName = PyUnicode_FromString((char*)"scipy"); //Issue Here
    check(pName, "pName not initializes.");

    PyObject *pModule = PyImport_Import(pName);
    check(pModule, "pModule not initializes.");

    PyObject *pFunc, *pArgs, *pResult;
    PyArrayObject *pNpArray;

    npy_intp Adims[2]; Adims[0] = A.rows(); Adims[1] = A.cols();
    npy_intp bdim[1];  bdim[0] = y.size();

    if (pModule != NULL) {
        pFunc = PyObject_GetAttrString(pModule, (char*)"optimize.nnls"); //Issue Here
        check(pFunc, "pFunc not initializes.");

        /*program never advances past this point unless I remove this check, 
        in which case I get a segfault, because pFunc is not initialized.*/

        if (pFunc && PyCallable_Check(pFunc)) {
            pArgs = PyTuple_Pack(2,
                PyArray_SimpleNewFromData(2, Adims, NPY_FLOAT, B.data()),
                PyArray_SimpleNewFromData(1, bdim , NPY_FLOAT, \
                    const_cast<double*> (y.data()))
            );
            check(pArgs, "pArgs not initializes.");
            pResult = PyObject_CallObject(pFunc, pArgs);
            check(pResult, "pResult not initializes.");
            pNpArray = reinterpret_cast<PyArrayObject*>(pResult);
            log_info("not PyArray cast");
        }
        Py_DECREF(pFunc);
    }
    Py_DECREF(pModule);
    Py_DECREF(pArgs);
    c_out = reinterpret_cast<long double*>(PyArray_DATA(pNpArray));
    for(size_t i=0; i<order; i++){
        gsl_vector_set(gslc, i, c_out[i]);
    }
    if (pResult != NULL) Py_DECREF(pResult); Py_DECREF(pNpArray);
    Py_Finalize();
}

同样，我的问题在于

PyObject *pModule = PyImport_Import(pName);

和

pFunc = PyObject_GetAttrString(pModule, (char*)"optimize.nnls");

如果我设置

pName = PyUnicode_FromString((char*)"scipy.optimize");

或

pFunc = PyObject_GetAttrString(pModule, (char*)"optimize.nnls");

代码无法在这些点初始化这些变量。

最后一次重申我的问题如何导入模块scipy.optimize或在此上下文中调用函数optimize.nnls？谢谢。我希望它不会太混乱，如果它让我知道并且我会澄清。

Answer 1

与很少使用的内置__import__函数一样，当PyImport_Import名称为scipy.optimize时，导入系统会加载scipy.optimize模块，但会返回{ {1}}代替scipy。您需要将属性链向下追溯到scipy.optimize：

scipy.optimize.nnls

与往常一样，当您使用它们时，不要忘记PyObject *scipy = PyImport_Import(pName); PyObject *optimize = PyObject_GetAttrString(scipy, "optimize"); PyObject *nnls = PyObject_GetAttrString(optimize, "nnls");您拥有的任何引用。其中包括Py_DECREF等内容。

Answer 2

首先，我从不使用PyImport_Import()。我一直使用PyImport_ImportModule()，它使用const char * arg而不是PyObject *。这应该会在代码中保存一些步骤。

PyObject *pModule = PyImport_ImportModule("scipy.optimize");
PyObject *pFunc = PyObject_GetAttrString(pModule, "nnls");

这就是我用我编写的代码做的事情。 PyImport_ImportModule()处理包层次结构。 PyObject_getAttrString()不知道如何处理点。

我从未使用过PyImport_Import()，因此我不知道其行为可能会如何改变。

当在另一个应用程序中嵌入python时，如何导入或调用子模块中的函数（即scipy.optimize.nnls）？

2 个答案: