[Numpy-discussion] Second try: possible bug in assignment to complex array

[josef.pktd at gmail.com]

On Fri, Aug 10, 2012 at 11:41 AM, <josef.pktd at gmail.com> wrote:

>
>
> On Fri, Aug 10, 2012 at 10:00 AM, Travis Oliphant <travis at continuum.io>wrote:
>
>>
>> On Aug 10, 2012, at 5:37 AM, Paul Anton Letnes wrote:
>>
>> >
>> >
>> > On 10. aug. 2012, at 09:54, Mark Bakker wrote:
>> >
>> >> I am giving this a second try. Can anybody help me out?
>> >>
>> >> I think there is a problem with assigning a 1D complex array of length
>> one
>> >> to a position in another complex array.
>> >>
>> >> Example:
>> >>
>> >> a = ones(1,'D')
>> >> b = ones(1,'D')
>> >> a[0] = b
>> >>
>> ---------------------------------------------------------------------------
>> >> TypeError                                 Traceback (most recent call
>> last)
>> >> <ipython-input-37-0c4fc6d780e3> in <module>()
>> >> ----> 1 a[0] = b
>> >>
>> >> TypeError: can't convert complex to float
>> >>
>> >> This works correctly when a and b are real arrays:
>> >>
>> >> a = ones(1)
>> >> b = ones(1)
>> >> a[0] = b
>> >>
>> >> Bug or feature?
>> >
>> > The exact same thing happens on OS X 10.7.4, python 2.7.3, numpy 1.6.1.
>> >
>> > Looks like a bug to me - or at least very surprising behavior.
>>
>> This is definitely an inconsistency.    The error seems more correct
>> (though the error message needs improvement).
>>
>> Can someone try this on NumPy 1.5 and see if this inconsistency existed
>> there as well.
>>
>
> >>> np.__version__
> '1.5.1'
>
> >>> a = np.ones(1,'D')
> >>> b = np.ones(1,'D')
> >>> a[0] = b
> Traceback (most recent call last):
>   File "<stdin>", line 1, in <module>
>
> TypeError: can't convert complex to float
> >>> a = np.ones(1)
> >>> b = np.ones(1)
> >>> a[0] = b
>

and

>>> a = np.ones(1,'D')
>>> b = 2*np.ones(1)
>>> a[0] = b
>>> a
array([ 2.+0.j])
>>> c = 3*np.ones(1, int)
>>> a[0] = c
>>> a
array([ 3.+0.j])



>
> Josef
>
>
>> Thanks,
>>
>> -Travis
>>
>> >
>> > Paul
>> > _______________________________________________
>> > NumPy-Discussion mailing list
>> > NumPy-Discussion at scipy.org
>> > http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>
>> _______________________________________________
>> NumPy-Discussion mailing list
>> NumPy-Discussion at scipy.org
>> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20120810/05588327/attachment.html>