0

根据 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)))

给出完全相同的结果。

4

1 回答 1

1

您需要restype像这样设置函数的:

libgsl.gsl_stats_correlation.restype = c_double

看看这个ctypes 教程(特别是这个部分)来了解如何正确指定参数和响应类型。

于 2017-10-12T01:40:43.840 回答