def test_callback_factory__numba():
args = x, y = symbols('x y')
expr = x + atan(y)
cb = _callback_factory(args, [expr], 'numpy', np.float64, 'C', use_numba=True)
n = 500
inp = np.empty((n, 2))
inp[:, 0] = np.linspace(0, 1, n)
inp[:, 1] = np.linspace(-10, 10, n)
assert np.allclose(cb(inp), inp[:, 0] + np.arctan(inp[:, 1]))
评论列表
文章目录