Question

我正在尝试用Cython包装一些C代码，但是我遇到了一个我不理解的错误，尽管进行了大量搜索，我似乎无法找到任何内容。这是我的c代码

void cssor(double *U, int m, int n, double omega, double tol, int maxiters, int *info){
double maxerr, temp, lcf, rcf;
int i, j, k;
lcf = 1.0 - omega;
rcf = 0.25 * omega;
for (k =0; k < maxiters ; k ++){
    maxerr = 0.0;
    for (j =1; j < n-1; j++) {
        for (i =1; i < m-1; i++) {
            temp = U[i*n+ j];
            U[i*n+j] = lcf * U[i*n+j] + rcf * (U[i*n+j-1] + U [i*n+j+1] + U [(i-1)*n + j] + U [(i+1)*n+j]);
            maxerr = fmax(fabs(U[i*n+j] - temp), maxerr);
        }
    }
    if(maxerr < tol){break;}
}
if (maxerr < tol) {*info =0;}
else{*info =1;}

}

我的.pyx文件是

    cdef extern from "cssor.h":
        void cssor(double *U, int m, int n, double omega, double tol, int maxiters, int *info)

    cpdef cyssor(double[:, ::1] U, double omega, double tol, int maxiters, int *info):
        cdef int n, m
        m = U.shape[0]
        n = U.shape[1]
        cssor(&U[0, 0], m, n, omega, tol, maxiters, &info)

但是，当我尝试运行相关的安装文件时，我在代码的最后一行中引用了一个错误，指的是：

无法分配类型＆＃39; int **＆＃39;输入＆＃39; int *＆＃39;

你能告诉我如何解决这个问题吗？

Roy Roth

Answer 1

问题来自：

cpdef cyssor(double[:, ::1] U, double omega, double tol, int maxiters, int *info):
    cdef int n, m
    m = U.shape[0]
    n = U.shape[1]
    cssor(&U[0, 0], m, n, omega, tol, maxiters, &info)

您将info声明为int*类型。但是，您可以将其作为cssor的引用传递给int*函数，使其成为int**。

正确的代码是：

cssor(&U[0, 0], m, n, omega, tol, maxiters, info)

Cython int **和int *类型

1 个答案: