为什么guvectorize在numba中返回奇怪的输出

Question

例如，我无法理解以下函数的输出：

@nb.guvectorize(["void(float64, float64, float64, float64)"], "(),()->(),()")
def add_subtract(a, b,  res1, res2):
res1 = a + b
res2 = a - b

当我称之为

add_subtract(np.array([3.0,5.0]), np.array([6.0,7.0]))

奇怪的是，它会显示：

(array([6.95307160e-310, 1.23643049e-311]), array([0., 0.]))

此函数似乎正在尝试返回0。为什么函数不返回[[9，-3]，[12，-2]]？

Answer 1

输出数组在未初始化的情况下传递给函数。如果为res1和res2局部变量指定一个新值，则将用新值替换数组。这不会更改传递给函数的数组，这些数组仍处于未初始化状态。

相反，您需要替换它们的内容。以下示例有效：

@nb.guvectorize("void(float64[:], float64[:], float64[:], float64[:])",
"(),()->(),()")
def add_subtract(a, b, res1, res2):
res1[:] = a + b
res2[:] = a - b