python - numpy.reciprocal 在重复调用时返回不同的值

Question

<分区>

Answer 1

这不仅仅是互惠；使用 where 参数时会出现此问题。我已经能够使用 numpy 的主分支(np.__version__ 是 '1.15.0.dev0+c093997')重现该问题，其功能如 abs, sign, add, subtract 等

如果您仔细阅读 numpy“ufuncs”的文档字符串并正确解释它们，您会发现该行为不是错误。以下是 numpy.reciprocal 文档字符串中的相关描述:

out : ndarray, None, or tuple of ndarray and None, optional
    A location into which the result is stored. If provided, it must have
    a shape that the inputs broadcast to. If not provided or `None`,
    a freshly-allocated array is returned. A tuple (possible only as a
    keyword argument) must have length equal to the number of outputs.
where : array_like, optional
    Values of True indicate to calculate the ufunc at that position, values
    of False indicate to leave the value in the output alone.

请特别注意:

• 哪里说“值 的 False 表示单独保留输出中的值。”
• out 表示“如果未提供或无， 返回一个新分配的数组。”

您没有提供 out 参数，因此通过调用 reciprocal 分配了一个新数组。这个数组的内容没有被初始化；该数组包含分配的内存中发生的所有内容。当您使用 where 参数时，只有输出中 where 为 True 的那些位置才被赋值。 where 为 False 的位置未被触及，因此它们保存数组分配时存在的任何随机内容。对于浮点输出，输出中的随机值可能是 0.0、4.85105226e-309 或任何其他随机值。

要按预期方式使用 where 参数，您还应该提供自己的 out 参数，并在输出 where where 中使用您想要的值进行初始化 为假。在您的情况下，您应该传入一个零数组:

In [84]: ssh_sum
Out[84]: 
array([[0., 2., 1., 0., 0., 0.],
       [0., 0., 1., 2., 0., 0.],
       [0., 0., 0., 1., 0., 2.]])

In [85]: out = np.zeros_like(ssh_sum)

In [86]: np.reciprocal(ssh_sum, where=ssh_sum > 0.0, out=out)
Out[86]: 
array([[0. , 0.5, 1. , 0. , 0. , 0. ],
       [0. , 0. , 1. , 0.5, 0. , 0. ],
       [0. , 0. , 0. , 1. , 0. , 0.5]])

In [87]: out
Out[87]: 
array([[0. , 0.5, 1. , 0. , 0. , 0. ],
       [0. , 0. , 1. , 0.5, 0. , 0. ],
       [0. , 0. , 0. , 1. , 0. , 0.5]])