[devel] atlas is broken again?

[Valery Pipin]

05.04.2011 22:15, REAL пишет:
> 06.04.2011 12:04, Valery Pipin пишет:
>> >>> m=randn(2000,2000)
>
> Странно, сейчас на этой строке не падает. Понятия не имею, что у Вас 
> за ситуация, она какая-то размытая.
>
Она никогда здесь не падала. Могу хоть 10000x10000 задать
In [1]: from numpy.linalg import *

In [2]: from numpy.random import *

In [3]: from numpy import dot

In [4]: m1=randn(5000,5000)

In [5]: %time m2=randn(5000,5000)
CPU times: user 1.90 s, sys: 0.33 s, total: 2.23 s
Wall time: 6.11 s

In [7]: %time m3=dot(m1,m2)
CPU times: user 104.00 s, sys: 0.67 s, total: 104.67 s
Wall time: 114.83 s

а 10000x10000 долго ждать в сизифе ;-)

А вот как EPD справляется:
Python 2.7.1 |EPD 7.0-2 (64-bit)| (r271:86832, Nov 29 2010, 13:51:37)
Type "copyright", "credits" or "license" for more information.

IPython 0.10.1 -- An enhanced Interactive Python.
?         -> Introduction and overview of IPython's features.
%quickref -> Quick reference.
help      -> Python's own help system.
object?   -> Details about 'object'. ?object also works, ?? prints more.

In [1]: from numpy.linalg import *

In [2]: from numpy.random import *

In [3]: from numpy import dot

In [4]: m1=randn(5000,5000)

In [5]: %time m2=randn(5000,5000)
CPU times: user 1.86 s, sys: 0.07 s, total: 1.93 s
Wall time: 1.93 s

In [7]: %time m3=dot(m1,m2)
CPU times: user 27.88 s, sys: 0.08 s, total: 27.95 s
Wall time: 16.12 s

Как говорится почувствуйте разницу ;-).