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 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
