[SciPy-user] Speed differences in sqrt calculation: what is good to know

David Cournapeau david@ar.media.kyoto-u.ac...
Thu Jun 12 22:23:00 CDT 2008

Lubos Vrbka wrote:
> hi,
> slightly off-topic, but still regarding the speed differences... what's 
> the difference between the following 2 functions?
> numpy.linalg.solve
> scipy.linalg.solve
> the first is almost 2 times faster on the same data (4096 repetitions of 
> a A*a = b problem with A = 2x2 float32 matrix) - typically 0.30 vs. 0.57 
> sec for the whole operation...

I think you will find a lot of those differences with very small data. I 
would not even be surprised if it was faster to do the inversion in 
python. Using numpy and scipy for those kind of data just do not make 
much sense. I am sure you will have exactly the same problem in matlab, 
albeit in a less significant way, because recent matlab engines have JIT 
(and a much simpler language, which certainly helps here).

Difference in functionalities between numpy and scipy are much more 
bothering to me (sqrt(-1) giving a different answer for numpy and scipy, 
for example).



More information about the SciPy-user mailing list