numba 中的函数类型

Question

目前在 numba 中处理高阶函数的最佳方法是什么？

def secant_method_curried (f):
    def inner (x_minus1, x_0, consecutive_tolerance):
        x_new = x_0
        x_old = x_minus1
        x_oldest = None
        while abs(x_new - x_old) > consecutive_tolerance:
            x_oldest = x_old
            x_old = x_new
            x_new = x_old - f(x_old)*((x_old-x_oldest)/(f(x_old)-f(x_oldest)))
        return x_new
    return numba.jit(nopython=False)(inner)

问题是没有办法告诉 numba f 是 doube(double)，所以上面的代码中断 nopython=True:

TypingError: Failed at nopython (nopython frontend)
Untyped global name 'f'

之前的版本好像有FunctionType，但是得到了removed/renamed：http://numba.pydata.org/numba-doc/0.8/types.html#functions

On this page，他们提到了一个叫做 numba.addressof() 的东西，这似乎有点帮助，但又可以追溯到 4 年前。

Answer 1

经过一些实验，我可以重现您的错误。在这种情况下，jit 传递给您的函数就足够了 secant_method_curried:

>>> from numba import njit

>>> def func(x):  # an example function
...     return x

>>> p = secant_method_curried(njit(func))  # jitted the function

>>> p(1,2,3)
2.0

你也可以在传入njit(func)或jit(func)时声明签名。

在 documentation 中也有一个很好的使用 numba 的闭包示例，并且还提到：

[...] you should JIT-compile that function if it is called from another jitted function.

numba 中的函数类型

function types in numba

python

numba