[Numpy-discussion] Square All Elements Of A Matrix

[andrea_gavana at tin.it]

Hello David,

>Could you check the speed of
>
>td = exp(2*log(dx))
>
>If it's comparable to dx**2, then your pow() library function isn't
>any good :-( -- it's not special-casing integer powers.

I get:

start = time.time(); td = Numeric.exp(2*Numeric.log(dx)); print time.time()
- start
0.530999898911

Does this make any sense to anyone? I am a bit puzzled...

Thanks however for your kind help/suggestions.

Andrea.