Re: [devel] atlas is broken again?

[Valery Pipin <pip@iszf.irk.ru>]

В сообщении от 5 апреля 2011 00:41:27 автор REAL написал:
> 05.04.2011 14:34, Valery Pipin пишет:
> >> Из сизифа?
> > 
> > Отсюда, http://www.enthought.com/
> 
> А зачем, если есть в сизифе?
> 
сизифовский медленный, и к тому же периодически ломается. В enthought сборке 
использована MKL и линейная алгебра реально паралеллится если есть куда.
> >> покажите rpm -qi libatlas
> > 
> > я писал что я вернулся е прежнюю, рабочую, версию
> 
> Тогда я вообще в растерянности, ибо ни провала в производительности,
> ни в падении ничего не заметил.
Возвращаемся к первому сообщению. Вот это с  libatlas-3.9.35:
IPython 0.10 -- 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.random import *

In [2]: from numpy.linalg import *

In [3]: m=randn(2000,2000)

In [4]: %time U,s,V=svd(m)
Ошибка сегментирования
Теперь пробуем для массивов меньших размеров:
In [1]: from numpy.linalg import *

In [2]: from numpy.random import *

In [3]: m=randn(100,100)

In [4]: %time U,s,V=svd(m)
CPU times: user 0.01 s, sys: 0.00 s, total: 0.01 s
Wall time: 0.01 s

In [6]: m=randn(200,200)

In [7]: %time U,s,V=svd(m)
---------------------------------------------------------------------------
LinAlgError                               Traceback (most recent call last)

/home/vv/<ipython console> in <module>()

/usr/lib/python2.6/site-packages/IPython/iplib.pyc in ipmagic(self, arg_s)
   1180         else:
   1181             magic_args = self.var_expand(magic_args,1)
-> 1182             return fn(magic_args)
   1183 
   1184     def ipalias(self,arg_s):

/usr/lib/python2.6/site-packages/IPython/Magic.pyc in magic_time(self, 
parameter_s)
   1969         else:
   1970             st = clk()
-> 1971             exec code in glob
   1972             end = clk()
   1973             out = None

/home/vv/<timed exec> in <module>()

/usr/lib64/python2.6/site-packages/numpy/linalg/linalg.pyc in svd(a, 
full_matrices, compute_uv)
   1318                                  work, lwork, iwork, 0)
   1319     if results['info'] > 0:
-> 1320         raise LinAlgError, 'SVD did not converge'
   1321     s = s.astype(_realType(result_t))
   1322     if compute_uv:

LinAlgError: SVD did not converge