Cython int ** and int * types

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)

但是，当我尝试运行相关的设置文件时，我在代码的最后一行中得到一个引用maxmaxters的错误，该错误说:

不能将'int **'类型赋值给'int *'类型

你能告诉我怎么修理这个吗?

罗斯罗伊

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*类型。但是您随后将其作为int*的引用传递给cssor函数，使其成为int**。

正确的代码是:

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