Yet another example of an unexpected behaviour:
>>> a=np.ma.array([], mask=0)
>>> b=np.ma.array([])
>>> np.ma.allequal(a,b)
True
>>> a.mean()
masked
>>> b.mean()
nan
But
>>>a
masked_array(data = [],
mask = [],
fill_value = 1e+20)
>>> b
masked_array(data = [],
f2py generated wrappers cannot be vectorized with numpy1.3.0 and Python2.6.2.
The reason is change to Python's getargs.c. Vectorize, or rather _get_nargs()
defined in lib/function_base.py tries to determine the number of arguments from
error message generated while the interpreter parses function i