python - numpy roll 和掩码数组的奇怪行为

Question

我有以下代码（即它只是一个演示我的问题的虚拟代码）

#!/usr/bin/python
import numpy as np

#define array, -999 being a missing value
a=np.array([5.,-999.,7.,8.])
print "array:"
print a
print

a=np.ma.masked_values(a,-999.)
print "masked array:"
print a
print

#define new missing value
bad=1e30
a.set_fill_value(bad)

#roll array
a=np.roll(a,1,axis=0)

#plug-in new missing value 
a=a.filled()

print "array after rolling with new missing value (1e30):"
print a

在我使用 numpy 1.4.1 的 32 位机器上，我得到了预期的输出：

array:
[   5. -999.    7.    8.]

masked array:
[5.0 -- 7.0 8.0]

array after rolling with new missing value (1e30):
[  8.00000000e+00   5.00000000e+00   1.00000000e+30   7.00000000e+00]

在我的带有 numpy 1.6.1 的 64 位机器上，我得到了

[   5. -999.    7.    8.]

masked array:
[5.0 -- 7.0 8.0]

array after rolling with new missing value (1e30):
[  8.00000000e+00   5.00000000e+00   1.00000000e+20   7.00000000e+00]

所以，我对 1e30 的掩码数组的填充值在滚动数组后更改为 1e20。我做错了什么，或者这是一个已知的错误？

拉斐尔