在GSL中使用ctypes传递数组

Question

根据GSL文件，签名为：

double gsl_stats_correlation (const double data1[], 
                              const size_t stride1, 
                              const double data2[], 
                              const size_t stride2, 
                              const size_t n)

当我试着从PyPy调用它时：

from ctypes import CDLL, RTLD_GLOBAL
gslcblas = CDLL('libgslcblas.0.dylib',mode=RTLD_GLOBAL)
libgsl = CDLL('/usr/local/lib/libgsl.0.dylib')
from ctypes import c_double, c_size_t, pointer
a1 = (c_double * 5)(1, 2, 3, 4, 5)
a2 = (c_double * 5)(1, 2, 3, 6, 5)
print(libgsl.gsl_stats_correlation(a1, c_size_t(1), 
      a2, c_size_t(1), c_size_t(5)))

目前我的机器上的结果是1086463496，尽管它会随着运行而变化。这显然是不正确的。我做错什么了？请注意，将函数调用更改为：

libgsl.gsl_stats_correlation(pointer(a1), c_size_t(1), 
                             pointer(a2), c_size_t(1), c_size_t(5)))

给出了完全相同的结果。你知道吗

Answer 1

您需要如下设置函数的restype：

libgsl.gsl_stats_correlation.restype = c_double

看看thisctypes教程（特别是this section），了解如何正确指定参数和响应类型。你知道吗