[Numpy-discussion] Arrays of poly1d objects, is this a bug?

[Stefan van der Walt]

On Fri, Jan 26, 2007 at 12:25:27PM -0700, Fernando Perez wrote:
> Hi all,
> 
> I'm puzzled by this behavior a colleague ran into:
> 
> In [38]: p1=N.poly1d([1.])
> 
> In [39]: a=N.array([p1],dtype='O')
> 
> In [40]: a
> Out[40]: array([], shape=(1, 0), dtype=object)
> 
> In [42]: print a
> []
> 
> In [43]: N.__version__
> Out[43]: '1.0.2.dev3512'
> 
> He saw it running r3520 as well.
> 
> This looks like a bug to me,  since it seems impossible to make arrays
> of poly1d objects:
> 
> In [44]: a=N.array([p1,p1],dtype='O')
> 
> In [45]: a
> Out[45]: array([], shape=(2, 0), dtype=object)
> 
> In [46]: print a
> []

I think the problem might be that an instance of poly1d can be viewed as an
array (e.g. N.asarray(p1)).

The following works:

x = N.empty((3,),object)
x[0] = p1

Cheers
Stéfan