[Numpy-discussion] Re: indexing problem

Gary Ruben gruben at bigpond.net.au
Mon Feb 13 20:50:02 CST 2006


Hi David,
So, I think what you had done would be OK provided you removed the 
x**0.5 case to avoid the problem Tim raised and checked that the 
exponent is an integer, not just a scalar.
Does anyone see a problem with this approach.
Gary R.

David M. Cooke wrote:
> Gary Ruben <gruben at bigpond.net.au> writes:
> 
>> Tim Hochberg wrote:
>> <snip>
>>> However, I'm not convinced this is a good idea for numpy. This would
>>> introduce a discontinuity in a**b that could cause problems in some
>>> cases. If, for instance, one were running an iterative solver of
>>> some sort (something I've been known to do), and b was a free
>>> variable, it could get stuck at b = 2 since things would go
>>> nonmonotonic there.
>> I don't quite understand the problem here. Tim says Python special
>> cases integer powers but then talks about the problem when b is a
>> floating type. I think special casing x**2 and maybe even x**3 when
>> the power is an integer is still a good idea.
> 
> Well, what I had done with Numeric did special case x**0, x**1,
> x**(-1), x**0.5, x**2, x**3, x**4, and x**5, and only when the
> exponent was a scalar (so x**y where y was an array wouldn't be). I
> think this is very useful, as I don't want to microoptimize my code to
> x*x instead of x**2. The reason for just scalar exponents was so
> choosing how to do the power was lifted out of the inner loop. With
> that, x**2 was as fast as x*x.
> 




More information about the Numpy-discussion mailing list