[Numpy-discussion] Unexpected output usingnumpy.ndarray and__radd__

Mark Hoffmann Mark.Hoffmann at dk.manbw.com
Mon Dec 18 07:20:57 CST 2006


Excellent, thank you - it solved the problem!

/Mark 

-----Original Message-----
From: numpy-discussion-bounces at scipy.org [mailto:numpy-discussion-bounces at scipy.org] On Behalf Of Tim Hochberg
Sent: 18. december 2006 14:06
To: Discussion of Numerical Python
Subject: Re: [Numpy-discussion] Unexpected output usingnumpy.ndarray and__radd__

Mark Hoffmann wrote:
> I appreciate the answer and the solution suggestion. I see that it is possible to make a work around by subclassing from ndarray. Still, in the "print a+c" statement, I don't understand why a.__add__(c) doesn't return NotImplemented (because ndarray shouldn't recognize the Cyclehist class) and directly call c.__radd__(a) implemented in my Cyclehist class. I tried the exactly same programme using Python 2.4.1 and Scipy 0.3.2 (based on numeric/numarray) and the result of the "print a+c" didn't get nested as I expect.
>
> Regards,
> Mark
>   
I'm not sure what this is doing -- it looks kind of bizzare -- however, you can fix this case without resorting to subclassing to ndarray. Just toss an '__array_priority__ = 10' up at the top of the class definition and it will use your __methods__ in preference to the ndarrays. I don't have time to look into this further right now unfortunately.

-tim


> -----Original Message-----
> From: numpy-discussion-bounces at scipy.org 
> [mailto:numpy-discussion-bounces at scipy.org] On Behalf Of Stefan van 
> der Walt
> Sent: 18. december 2006 10:36
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] Unexpected output using numpy.ndarray 
> and__radd__
>
> Hi Mark
>
> On Mon, Dec 18, 2006 at 08:30:20AM +0100, Mark Hoffmann wrote:
>   
>> The following issue has puzzled me for a while. I want to add a 
>> numpy.ndarray and an instance of my own class. I define this 
>> operation by implementing the methods __add__ and __radd__. My programme (including output) looks like:
>>
>> #!/usr/local/bin/python
>>
>> import numpy
>>
>> class Cyclehist:
>>     def __init__(self,vals):
>>         self.valuearray = numpy.array(vals)
>>
>>     def __str__(self):
>>         return 'Cyclehist object: valuearray = '+str(self.valuearray)
>>
>>     def __add__(self,other):
>>         print "__add__ : ",self,other
>>         return self.valuearray + other
>>    
>>     def __radd__(self,other):
>>         print "__radd__ : ",self,other
>>         return other + self.valuearray
>>
>> c = Cyclehist([1.0,-21.2,3.2])
>> a = numpy.array([-1.0,2.2,-2.2])
>> print c + a
>> print a + c
>>     
>
> In the first instance, c.__add__(a) is called, which works fine.  In the second, a.__add__(c) is executed, which is your problem, since you rather want c.__radd__(a) to be executed.  A documentation snippets:
>
> """For instance, to evaluate the expression x-y, where y is an instance of a class that has an __rsub__() method, y.__rsub__(x) is called if x.__sub__(y) returns NotImplemented.
>
> Note: If the right operand's type is a subclass of the left operand's type and that subclass provides the reflected method for the operation, this method will be called before the left operand's non-reflected method. This behavior allows subclasses to override their ancestors' operations."""
>
> Since a.__add__ does not return NotImplemented, c.__radd__ is not called where you expect it to be.  I am not sure why broadcasting takes place here, maybe someone else on the list can elaborate.
>
> To solve your problem, you may want to look into subclassing ndarrays, as described at http://www.scipy.org/Subclasses.
>
> Cheers
> Stéfan
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at scipy.org
> http://projects.scipy.org/mailman/listinfo/numpy-discussion
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at scipy.org
> http://projects.scipy.org/mailman/listinfo/numpy-discussion
>
>
>   


_______________________________________________
Numpy-discussion mailing list
Numpy-discussion at scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion


More information about the Numpy-discussion mailing list