Numba中的发电机参数

时间:2017-01-25 20:22:56

标签: python generator numba

此问题的后续跟进:function types in numba

我正在编写一个需要将生成器作为其参数之一的函数。粘贴在这里太复杂了,所以请考虑这个玩具示例:

def take_and_sum(gen):
    @numba.jit(nopython=False)
    def inner(n):
        s = 0
        for _ in range(n):
            s += next(gen)
        return s
    return inner

它返回生成器的第一个n元素的总和。用法示例:

@numba.njit()
def odd_numbers():
    n = 1
    while True:
        yield n
        n += 2

take_and_sum(odd_numbers())(3) # prints 9

这是因为我想用nopython=True进行编译,然后我不能将gen(一个pyobject)作为参数传递。不幸的是,nopython=True我收到错误:

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

即使我nopython编译了我的生成器。

令人困惑的是,对输入进行硬编码工作:

def take_and_sum():
    @numba.njit()
    def inner(n):
        gen = odd_numbers()
        s = 0.0
        for _ in range(n):
            s += next(gen)
        return s
    return inner

take_and_sum()(3)

我也尝试将我的发电机变成一个类:

@numba.jitclass({'n': numba.uint})
class Odd:
    def __init__(self):
        self.n = 1
    def next(self):
        n = self.n
        self.n += 2
        return n

同样,这在对象模式下工作,但在nopython模式下我得到了不可搜索的:

LoweringError: Failed at nopython (nopython mode backend)
Internal error:
NotImplementedError: instance.jitclass.Odd#4aa9758<n:uint64> as constant unsupported

1 个答案:

答案 0 :(得分:2)

我实际上无法解决您的问题,因为据我所知,这是不可能的。我只是突出了一些方面(对numba 0.30有效):

无法创建一个numba - jitclass生成器:

import numba

@numba.jitclass({'n': numba.uint})
class Odd:
    def __init__(self):
        self.n = 1

    def __iter__(self):
        return self

    def __next__(self):
        n = self.n
        self.n += 2
        return n

试试吧:

>>> next(Odd())
TypeError: 'Odd' object is not an iterator

当您删除时,numba.jitclass可以正常运行:

>>> next(Odd())
1

使用硬编码生成器的示例不相同。您的原始尝试创建一个生成器对象将其传递给numba函数并修改生成器。您可能希望它更新生成器的状态

>>> t = odd_numbers()
>>> take_and_sum(t)(3)
9
>>> next(t)   # State has been updated, unfortunatly that requires nopython=False!
7

但这对于numba来说根本不可能。(

第二个例子不同,因为每次调用函数时都会创建生成器,因此函数外部没有需要更新的状态:

>>> take_and_sum()(3) # using your hardcoded version
9.0
>>> take_and_sum()(3) # no updated state so this returns the same:
9.0

完全可以改变它,但没有选择使用仲裁功能:

@numba.jitclass({'n': numba.uint})
class Odd:
    def __init__(self):
        self.n = 1

    def calculate(self, n):
        s = 0.0
        for _ in range(n):
            s += self.n
            self.n += 2
        return s

>>> x = Odd()
>>> x.calculate(3)
9.0
>>> x.calculate(3)
27.0

我知道这不是你想要的,但至少它在某种程度上有效: - )